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Abstract 


Sorting algorithms have attracted a great deal of attention and study, as they 
have numerous applications to Mathematics, Computer Science and related 
fields. In this thesis, we first deal with the mathematical analysis of the Quick¬ 
sort algorithm and its variants. Specifically, we study the time complexity of 
the algorithm and we provide a complete demonstration of the variance of the 
number of comparisons required, a known result but one whose detailed proof 
is not easy to read out of the literature. We also examine variants of Quicksort, 
where multiple pivots are chosen for the partitioning of the array. 

The rest of this work is dedicated to the analysis of finding the true order by 
further pairwise comparisons when a partial order compatible with the true 
order is given in advance. We discuss a number of cases where the partially 
ordered sets arise at random. To this end, we employ results from Graph and 
Information Theory. Finally, we obtain an alternative bound on the number of 
linear extensions when the partially ordered set arises from a random graph, 
and discuss the possible application of Shellsort in merging chains. 
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Chapter 1 


Preface 


The first Chapter serves as an introduction to this work, where in a simple 
manner, useful notions of algorithmic analysis in general are presented, along 
with definitions that will be used throughout the thesis. We present an intro¬ 
duction to Quicksort and this Chapter ends with an outline and a summary of 
the main contributions of this thesis. 


1.1 Preliminaries 

An algorithm is a valuable tool for solving computational problems. It is a 
well defined procedure that takes some value or a set of values as input and is 
guaranteed to produce an answer in finite time. See e.g. [15]. However the 
time taken may be impractically long. 

As a useful introduction, we present a simple and intuitive algorithm known as 
Euclid’s algorithm which is still used today to determine the greatest common 
divisor of two integers. Its definition follows [45]: 

Definition 1.1.1. Given two integers A > C greater than unity, we want to find 
their greatest common divisor, g.c.d.(A, C). 


1 



2 


1. IfC divides A, then the algorithm terminates with C as the greatest common 
divisor. 

2. If A mod C is equal to unity, the numbers are either prime or relatively 
prime and the algorithm terminates. Othenvise set A •<— C, C •<— A mod C 
and return to step 1. 

By A mod C, we denote the remainder of the division of A by C, namely 


A mod C A - I (7 • - 


which is between 0 and (7 — 1. Here the floor function [xj of a real number 
x is the largest integer less than or equal to x. We observe that this algorithm 
operates recursively by successive divisions, until obtaining remainder equal 
to 0 or to 1. Since the remainder strictly reduces at each stage, the process is 
finite and eventually terminates. 

Two main issues in relation to an algorithm are its running time or time com¬ 
plexity, which is the amount of time necessary to solve the problem, and its 
space complexity, which is the amount of memory needed for the execution 
of the algorithm in a computer. In case of Euclid’s algorithm, 3 locations of 
memory are required for storing the integer numbers A, C and A mod C. In 
this thesis, we will mainly be concerned with time complexity. 

The aim is usually to relate the time complexity to some measure of the size of 
the instance of the problem we are considering. Very often we are particularly 
concerned with what happens when the size n is large - tending to infinity. 
The notations O, H, 0 and o, u> are useful in this context. 
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1.2 Definitions 

In this section, we present primary definitions to the analysis presented in the 
thesis. These definitions come from [15], [26]. 

Definition 1.2.1. Let f{n ) and g{n) be two functions with n £ N. We say that 
f(n) = 0(g(n )) if and only if there exists a constant c > 0 and n 0 £ N, such that 

I/HI < c- | g(n)\, Vn > n 0 . 

Definition 1.2.2. Also, we say that f(n) = f7 (< 7 ( 71 )) if and only if there exists a 
constant c > 0 and n 0 e N, such that \f(n)\ > c ■ \g(n)\, Vn > n 0 . 

Definition 1.2.3. Furthermore, we say that f(n) = Q[g(n)) if and only if there 
exist positive constants c 1} c 2 andn 0 £ N, such that cr\g(n)\ < \f(n)\ < c 2 -\g(n)\, 
Vn > no- Equivalently, we can state that if f(n) = 0(g(n)) and f(n) = M(g(n)), 
then f(n) = Q(g(n)). 

We will also use the notations o, u. They provide the same kind of limiting 
bounds with the respective upper case notations. The difference is that for 
two functions f(n) and g(n), the upper case notation holds when it does exist 
some positive constant c. Whereas, the respective lower case notation is true 
for every positive constant c [15]. In other words, o is stronger statement than 
O, since f(n) = o(g(n)) implies that / is dominated by g. Equivalently, this 
can be stated as 


fin) = o(g(n)) 


lim 

n—» 00 


fjn) 

g(n) 


0, 


provided that g(n) is non-zero. 
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The relation f(n) = u(g(n)) implies that / dominates g, i.e. 


f(n) = u>(g{n)) 


lim 

n—>■ oo 


fH 

g(n) 


= oo. 


The relation f(n) ~ g(n) denotes the fact that f(n) and g(n) are asymptotically 
equivalent, i.e. 


lim 

n—» oo 


f(n) 

d(n) 


1. 


Further, best, worst and average case performance denote the resource us¬ 
age, e.g. amount of memory in computer or running time, of a given algo¬ 
rithm at least, at most and on average, respectively. In other words, these 
terms describe the behaviour of an algorithm under optimal circumstances 
(e.g. best-case scenario), worst circumstances and on average [15]. In this 
work, the study will be concentrated on the average and worst case perfor¬ 
mance of Quicksort and its variants. 


In our analysis, we will frequently come across with harmonic numbers, whose 
definition we now present. 

Definition 1.2.4. The sum 


:= V 


1 

i k 


11 1 
l fc + 2 k + + n* 


is defined to be the generalised n t h harmonic number of order k. When k — 1, the 
sum denotes the n th harmonic number, which we simply write H n . We define also 
H 0 := 0. 

There are numerous interesting properties of harmonic numbers, which are 
not yet fully investigated and understood. Harmonic series have links with 
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Stirling numbers [26] and arise frequently in the analysis of algorithms. For n 
large, it is well-known that [44], 


H n = log e (n) +7 + ±- J l- 


+ 


120 n 4 


+ 0 




where 7 = 0.57721... is the Euler-Mascheroni constant. We will use this often, 
especially in the form H n = log e (n) + 7 + o(l). 

Note that throughout this thesis, we shall adopt the convention of writing 
explicitly the base of logarithms. For example, the natural logarithm of n is 
denoted by log e (n), instead of ln(n). Also, the end of a proof will be denoted 
by the symbol ■ . 


1.3 Introduction to Quicksort 

Sorting an array of items is clearly a fundamental problem, directly linked 
to efficient searching with numerous applications. The problem is that given 
an array of keys, we want to rearrange these in non-decreasing order. Note 
that the order may be numerical, alphabetical or any other transitive relation 
defined on the keys [46]. In this work, the analysis deals with numerical order, 
where the keys are decimal numbers and we particularly focus on Quicksort 
algorithm and variants of it. Quicksort was invented by C. A. R. Hoare [29, 
30], Here is the detailed definition. 

Definition 1.3.1. 

The steps taken by the Quicksort algorithm are: 


1. Choose an element from the array, called pivot. 
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2. Rearrange the array by comparing every element to the pivot, so all elements 
smaller than or equal to the pivot come before the pivot and all elements 
greater than or equal to the pivot come after the pivot. 

3. Recursively apply steps 1 and 2 to the subarray of the elements smaller than 
or equal to the pivot and to the subarray of the elements greater than or 
equal to the pivot. 

Note that the original problem is divided into smaller ones, with (initially) 
two subarrays, the keys smaller than the pivot, and those bigger than it. Then 
recursively these are divided into smaller subarrays by further pivoting, until 
we get trivially sorted subarrays, which contain one or no elements. Given 
an array of n distinct keys A = (ai, a 2 ,..., a n } that we want to quick sort, 
with all the n\ permutations equally likely, the aim is to finding the unique 
permutation out of all the n\ possible, such that the keys are in increasing 
order. The essence of Quicksort is the partition operation, where by a series 
of pairwise comparisons, the pivot is brought to its final place, with smaller 
elements on its left and greater elements to the right. Elements equal to pivot 
can be on either or both sides. 

As we shall see, there are numerous partitioning schemes, and while the details 
of them are not central to this thesis, we should describe the basic ideas. A 
straightforward and natural way (see e.g. [46]) uses two pointers - a left 
pointer, initially at the left end of the array and a right pointer, initially at the 
right end of the array. We pick the leftmost element of the array as pivot and 
the right pointer scans from the right end of the array for a key less than the 
pivot. If it finds such a key, the pivot is swapped with that key. Then, the 
left pointer is increased by one and starts its scan, searching for a key greater 
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than the pivot: if such a key is found, again the pivot is exchanged with it. 
When the pointers are crossed, the pivot by repeated exchanges will “float” to 
its final position and the keys which are on its left are smaller and keys on its 
right are greater. The data movement of this scheme is quite large, since the 
pivot is swapped with the other elements. 

A different partitioning scheme, described in [30] is the following. Two point¬ 
ers i (the left pointer, initially 1) and j (the right pointer, initially n) are set 
and a key is arbitrarily chosen as pivot. The left pointer goes to the right un¬ 
til a key is found which is greater than the pivot. If one is found, its scan is 
stopped and the right pointer scans to the left until a key less than the pivot is 
found. If such a key is found, the right pointer stops and those two keys are 
exchanged. After the exchange, both pointers are stepped down one position 
and the lower one starts its scan. When pointers are crossed, i.e. when i > j, 
the final exchange places the pivot in its final position, completing the parti¬ 
tioning. The number of comparisons required to partition an array of n keys is 
at least n — 1 and the expected number of exchanges is f + Jp 

A third partitioning routine, called Lomuto’s partition, is mentioned in [6] - 
this involves exactly n — 1 comparisons, which is clearly best possible, but the 
downside is the increased number of exchanges. The expected number of key 
exchanges of this scheme is [48]. 

We now consider the worst case and best case, analysis of Quicksort. Suppose 
we want to sort the following array, (oi < a 2 < ... < a n } and we are very 
unlucky and our initial choice of pivot is the largest element a n . Then of 
course we only divide and conquer in a rather trivial sense: every element is 
below the pivot, and it has taken us n — 1 comparisons with a n to get here. 
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Suppose we now try again and are unlucky again, choosing a n _i as pivot this 
time. Again the algorithm performs n — 2 comparisons and we are left with 
everything less than a„_i. If we keep being unlucky in our choices of pivot, 
and keep choosing the largest element of what is left, after i recursive calls the 
running time of the algorithm will be equal to (n — 1) + (n — 2) + ... + (n — i ) 
comparisons, so the overall number of comparisons made is 

. . n ■ (n — 1) 

1 + 2 + • • • + (n — 1) =---. 

Thus Quicksort needs quadratic time to sort already sorted or reverse-sorted 
arrays if the choice of pivots is unfortunate. 

If instead we always made good choices, choosing each pivot to be roughly in 
the middle of the array we are considering at present, then in the first round 
we make n — 1 comparisons, then in the two subarrays of size about n/2 we 
make about n/2 comparisons, then in each of the four subarrays of size about 
n/4 we make n /4 comparisons, and so on. So we make about n comparisons 
in total in each round. The number of rounds will be roughly log 2 (n) as we are 
splitting the arrays into roughly equally-sized subarrays at each stage, and it 
will take log 2 (n) recursions of this to get down to trivially sorted arrays. 

Thus, in this good case we will need 0(nlog 2 (n)) comparisons. This is of 
course a rather informal argument, but does illustrate that the time complexity 
can be much smaller than the quadratic run-time in the worst case. This is 
already raising the question of what the typical time complexity will be: we 
address this in the next Chapter. 
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We briefly discuss the space complexity of the algorithm. There are n memory 
locations occupied by the keys. Moreover, the algorithm, due to its recursive 
nature, needs additional space for the storage of subarrays. The subarrays’ 
boundaries are saved on to a stack, which is a data structure providing tem¬ 
porary storage. At the end of the partition routine, the pivot is placed in its 
final position between two subarrays (one of them possibly empty). Recur¬ 
sively, the algorithm is applied to the smaller subarray and the other one is 
pushed on to stack. Since, in best and average case of Quicksort, we have 
0(log 2 (n)) recursive calls, the required stack space is 0(log 2 (n)) locations in 
memory. However, in worst case the stack may require 0(n ) locations, if the 
algorithm is applied to the larger subarray and the smaller one is saved to the 
stack [63]. 

This discussion makes it clear that the pivot selection plays a vital role in 
the performance of the algorithm. Many authors have proposed various tech¬ 
niques to remedy this situation and to avoid worst case behaviour, see [15], 
[30], [46], [62], [63] and [68]. These include the random shuffling of the 
array prior to initialisation of the algorithm, choosing as pivot the median of 
the array, or the median of a random sample of keys. 

Scowen in his paper [62], suggested choosing as pivot the middle element of 
the array: his variant is dubbed “Quickersort”. Using this rule for the choice 
of partitioning element, the aim is the splitting of the array into two halves 
of equal size. Thus, in case where the array is nearly sorted, quadratic time 
is avoided but if the chosen pivot is the minimum or maximum key, the al¬ 
gorithm’s running time attains its worst case and this variant does not offer 
any more than choosing the pivot randomly. Singleton [68] suggested a better 
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estimate of the median, by selecting as pivot the median of leftmost, rightmost 
and middle keys of the input array. Hoare [30] suggested the pivot may be 
chosen as the median of a random sample from the keys to be sorted, but he 
didn’t analyse this approach. 

One point is that Quicksort is not always very fast at sorting small arrays. 
Knuth [46] presented and analysed a partitioning scheme, which takes n + 
1 instead of n — 1 comparisons and the sorting of small subarrays (usually 
from about 9 to 15 elements) is implemented using insertion sort, since the 
recursive structure of Quicksort is better suited to large arrays. Insertion sort 
is a simple sorting algorithm, which gradually ‘constructs’ a sorted array from 
left to right, in the following manner. The first two elements are compared and 
exchanged, in case that are not in order. Then, the third element is compared 
with the element on its left. If it is greater, it is left at its initial location, 
otherwise is compared with the first element and accordingly is inserted to its 
position in the sorted array of 3 elements. This process is iteratively applied to 
the remaining elements, until the array is sorted. See as well in Cormen et al. 
[15], for a description of the algorithm. 


1.4 Outline and contributions of thesis 

This thesis consists of seven Chapters and one Appendix. After the first, intro¬ 
ductory Chapter, the rest of the thesis is organised as follows: 

In Chapter 2, we consider the first and second moments of the number of 
comparisons made when pivots are chosen randomly. The result for the mean 
is known and easy: the result for the variance is known, but less easy to find a 
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full proof of in the literature. We supply one. We briefly discuss the skewness 
of the number of comparisons and we study the asymptotic behaviour of the 
algorithm. 

In Chapter 3, we analyse the idea of choosing the pivot as a centered statistic 
of a random sample of the keys to be sorted and we obtain the average number 
of comparisons required by these variants, showing that the running time can 
be greatly improved. Moreover, we present relevant definitions of entropy. 

Not much of this is original, but some details about why various differential 
equations that arise in the analysis have the solutions they do (i.e. details 
about roots of indicial polynomials) are not in literature. 

In Chapter 4, we analyse extensions of Quicksort, where multiple pivots are 
used for the partitioning of the array. The main contributions in this Chapter 
are in sections 4.1 and 4.2. The results in the former section were published 
in the paper [34], where the expected costs related to the time complexity and 
the second moment of the number of comparisons are computed. The latter 
section contains the analysis of the generalisation of the algorithm. We study 
the general recurrence model, giving the expected cost of the variant, provided 
that the cost during partitioning is linear, with respect to the number of keys. 

We also present the application of Vandermonde matrices for the computation 
of the constants involved to the cost of these variants. 

In Chapter 5, various cases of partially ordered sets are discussed and the 
number of comparisons needed for the complete sorting is studied. The ‘information- 
theoretic lower bound’ is always u(n) in these cases and we show that the time 
needed for the sorting of partial orders is 0(nlog 2 (n)). The main contribution 
of this Chapter, is the derivation of the asymptotic number of comparisons 
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needed, for the sorting of various partially ordered sets. The basic ideas used 
here are due to, amongst others, Cardinal et al. [13], Fredman [25], Kahn and 
Kim [40], Kislitsyn [41], but the working out of the detailed consequences for 
these partial orders seems to be new. 

In Chapter 6, we consider random graph orders, where the ‘information- 
theoretic lower bound’ is of the same order of magnitude as the number of 
keys being sorted. We derive a new bound on the number of linear extensions 
using entropy arguments, though it is not at present competitive with an older 
bound in the literature [4]. 

In Chapter 7, we conclude the thesis, presenting future research directions. 
At this final Chapter, we derive another bound of the number of comparisons 
required to sort a random interval order and we discuss the merging of linearly 
ordered sets. 

In Appendix A, we present the Maple calculations, regarding the derivation 
of the variance of the number of comparisons of dual pivot Quicksort, analysed 
in subsection 4.1.1. 



Chapter 2 


Random selection of pivot 


In this Chapter, the mathematical analysis of Quicksort is presented, under 
the assumption that the pivots are uniformly selected at random. Specifically, 
the major expected costs regarding the time complexity of the algorithm and 
the second moment are computed. The derivation of the average costs is uni¬ 
fied under a general recurrence relation, demonstrating the amenability of 
the algorithm to a complete mathematical treatment. We also address the 
asymptotic analysis of the algorithm and we close this Chapter considering the 
presence of equal keys. 


2.1 Expected number of comparisons 

This discussion of lucky and unlucky choices of pivot suggests the idea of se¬ 
lecting the pivot at random, as randomisation often helps to improve running 
time in algorithms with bad worst-case, but good average-case complexity 
[69]. For example, we could choose the pivots randomly for a discrete uni¬ 
form distribution on the array we are looking at each stage. Recall that the 
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uniform distribution on a finite set assigns equal probability to each element 
of it. 

Definition 2.1.1. C n is the random variable giving the number of comparisons 
in Quicksort of n distinct elements when all the n\ permutations of the keys are 
equiprobable. 

It is clear that for n — 0 or n — 1, C 0 — Ci = 0 as there is nothing to sort. 
These are the initial or “seed” values of the recurrence relation for the number 
of comparisons, given in the following Lemma. 

Lemma 2.1.2. The random number of comparisons C n for the sorting of an array 
consisting ofn > 2 keys, is given by 

C n = Cjj n —1 + C*_ Un +71—1, 

where U n follows the uniform distribution over the set { 1 , 2 ,..., n} and C*_ v is 
identically distributed to C Un ~i and independent of it conditional on U n . 

Proof. The choice of U n as pivot, and comparing the other n — 1 elements 
with it, splits the array into two subarrays. There is one subarray of all U n — 1 
elements smaller than the pivot and another one of all n — U n elements greater 
than the pivot. Obviously these two subarrays are disjoint. Then recursively 
two pivots are randomly selected from the two subarrays, until the array is 
sorted, and so we get the equation. ■ 
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This allows us to find that the expected complexity of Quicksort applied to n 
keys is: 


E(C„) = E(C k „-i + C* n _ Un + n - 1) 

= E(C I/ ,_ 1 )+E(c;_ [/ J+ I i-1. 


Using conditional expectation and noting that U n = k has probability 1 /n, we 
get, writing a k for E (C k ), that 


^ 1 

a n = y - ' —(a k -1 + a n -k ) + n 


k =1 


n 


2 n_1 

a n \ T n 


n 


k =0 


1 . 


We have to solve this recurrence relation, in order to obtain a closed form for 
the expected number of comparisons. The following result is well-known (e.g. 
see in [15], [63]): 

Theorem 2.1.3. The expected number a n of comparisons for Quicksort with uni¬ 
form selection of pivots is a n = 2 (n + 1 )H n — 4 n. 


Proof. We multiply both sides of the formula for a n by n, getting 

n— 1 

na n = 2 a k + n(n — 1) 
k =o 

and similarly, multiplying by n — 1, 

n —2 

(■n - l)a„_i = 2 ^ a fc + (n - 2)(n - 1). 

k =0 
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Subtracting (n — l)a n _i from na n in order to eliminate the sum - see [46], we 
obtain 


na r 


- (n - 1 )a n _i = 2a n _i + 2(n - 1) 


na n = (n + l)a n _i + 2(n - 1) 

1 2(77 1 ) 

71+1 77 77(77 + 1) 


“Unfolding” the recurrence we get 


77 + I 


2 Xi ~ T ~~1 ~ 2 XI (—--- 


= 4 Hn + 


— 2 H n + 


77+1 

4 

77+1 


4 - 2 + 


4 . 


Finally, 


a n = 2(77 + !)//„ - 4t7. 


We now show a slick way of solving the recurrence about the expected number 
of comparisons using generating functions. This approach is also noted in 
various places, e.g. [44], [55], [63], We again start from 


n— 1 


Q"n. 


n 


J2 a j + n- 1 . 


i=o 
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We multiply through by n to clear fractions, getting 


n— 1 

na n = n{n — 1) + 2 dj 
3=0 

CO CO n— 1 

- 1 = n ( n - + 2 ^ ^ ajX n ~ l 

n= 1 n=l j=0 

co co n —1 

, na n x n ~ l = x ^ n(n — l)^ -2 + 2 ^ ^ djX n ~ l . 

n=l n =2 n=l j=0 


E 

n=l 

co 

E 


nd„x n = 


Letting /(x) = XT=o a n^ n , 


CO CO 


/'(x) = x^n{n- l)x n 2 + 2^2 

j= 0 n=J+l 


a.,a; n \ 


n=2 


co 

changing round the order of summation. To evaluate ^ n(n — l)a; n ~ 2 , sim- 

n =2 
co 

ply note this is the 2nd derivative of ^a; n . Since the latter sum is equal to 

n =0 

1/(1 — x), for |x| < 1 , its second derivative is easily checked to be 2(1 — a;) -3 . 
Multiplying the sum by x now gives 


f'{x) = 2x(l — x) 3 + 2 ^ dj ^ x n 1 

j =o n=j+1 


We evaluate the last double-sum. The inner sum is of course x j /(I — x) being 
a geometric series. Thus we get 

2 °° 

f(x) = 2a;(l — x)~ 3 + -—— ^ djX j 

~ X 3=0 

= 2x(l - a;) -3 + —— f(x). 

1 — x 
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This is now a fairly standard kind of differential equation. Multiplying both 
sides by (1 — x ) 2 , we see 

(1 — x) 2 f(x) = 2x{l — x ) _1 + 2(1 — x)f(x) 

==►(1 - x) 2 f'(x) - 2(1 - x)f(x) = —^-2 

1 — x 

^ (( 1 - x) 2 f(x))' - —^ - 2 
=>(1 — x) 2 f(x) = —2 log e (l — x) — 2x + c. 


Setting x = 0, we get /(0) =0 on the left-hand side, and on the right-hand 
side we get c, so c = 0. Therefore 


f(x) 


—2^1og e (l - x) + x^j 
(1 — x) 2 


Expanding out log e (l — x) as a series, — x — x 2 /2 — x 3 /3 — ..., and similarly 
writing 1/(1 — x) 2 as the derivative of 1/(1 — x), we obtain 

OO fc OO 

f(x) = + 

k=2 j =0 

Thus, looking at coefficients of x n on both sides, we get on the left-hand side 
a n . On the right-hand side, we get the coefficient for each x j in the first series 
(which is 1 / j) times the term for the x n ~ j in the other, namely (n — j + 1). So 



19 


we get 


a n — 


n 7,1 n ,1 n 

2 E n-k ± l =2E n + l_ 2l:i 


k =2 


k =2 


k=2 


1 


2(n + i) Et:- 1 - 2 (^- 1 ) 


k 

'k=l 

= 2 (n + i)H n — An. 


Some texts give a different argument for this, as follows. 

Theorem 2.1.4 (Mitzenmacher and Upfal [54]). Suppose that a pivot is cho¬ 
sen independently and uniformly at random from an array of n keys, in which 
Quicksort is applied. Then, for any input, the expected number of comparisons 
made by randomised Quicksort is 2nlog e (n) + 0(n). 

Proof. Let {xi,x 2 ,.. . ,x n } be the input values and the output, after the ter¬ 
mination of Quicksort, (i.e. keys in increasing order) be a permutation of 
the initial array { 2 / 1 , 2 / 2 , • • • ,2 In}- For * < j, let X tJ be a {0, l}-valued random 
variable, that takes the value 1, if y % and y 3 are compared over the course of 
algorithm and 0, otherwise. Then, the total number of comparisons X satisfies 

n— 1 n 

x = E E 

i =1 j=i +1 

( n —1 n \ n— 1 n 

E E x ‘i LEE E ( x «)- 

i= 1 j=i +1 / i= 1 j=i -\-1 

Since X,j e {0,1}, then E(X,j) = P(X^- = 1). This event occurs when y, and 
Dj are compared. Clearly, this happens when y, or y 3 is the first number chosen 
as pivot from the set E = {y», y»+i, • • • ,yj_i,yj}. (Because otherwise if some 
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element between them is chosen, y k say, y t and y 3 are compared to y k and y t is 
put below y k and y 3 above it with the result that y { and y 3 are never compared). 


Since the pivot is chosen uniformly at random, the probability the one of these 

1 

two elements is the first of the j — i +1 elements chosen is equal to-h 

j-i + 1 

1 2 

-=-. Substituting k = j — i + 1, we obtain 

j — i + 1 j — i + 1 


f n— 1 n 


E(X) = E EE; 


, 2=1 j=i +1 


j-i + 1 


n— 1 n— 2+1 

EE 


i =1 k =2 


2 

k 


n n-\-l—k 

E E 


k =2 i=l 


2 

k' 


The change of sum is justified by writing out the possible rows for fixed i 
and columns for a fixed k, then changing round from summing rows first to 
summing columns first. This is, as 2 /k does not depend on i, equal to 


n ~ n 0 n 

^(n + l-A;)- = (n + l)^--2^1. 


Thus, 


1 

2 (n + 1) V - - 2n + 2 = 2 (n + 1 )(H n - 1) - 2n + 2 
k 

k =2 

= 2 (n + l)i? n — 4n, 


proving the claim. ■ 

Also of some interest is the mean number of partitioning stages E (P n ) of the 
algorithm applied to n keys as input. For the case where we simply use Quick¬ 
sort for all the sorting, it is obvious that we will have P 0 = 0 and for n > 1, the 
number of partitioning stages P n obeys the following recurrence conditional 
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on j being the rank of the pivot 


Pn — 1 + Pj -1 + Pn-j- 

Observe that P\ — l and P 2 = 2. Taking expectations, and noting that the 
pivot is uniformly chosen at random and that the two sums E(P J -_ 1 ) and 
E( P n -j ) are equal, we see 

2 n ^ 

E(P„) = 1 + - VE(P i )- 

n 

3=0 

Multiplying both sides by n and differencing the recurrence relation, as we did 
for the derivation of the expected number of comparisons, we have 


nE(P„) - (n - l)E(P n _i) 
__ E(P n ) 
n + 1 

=► E(P n ) 


1 + 2E(P n _!) 


1 | E(P n _i) 

n(n + 1) n 



Finally, E(P n ) = n. 


2.2 Expected number of exchanges 

Here we consider the number of exchanges or swaps performed by the algo¬ 
rithm, which is mentioned by Hoare [30] as a relevant quantity. We assume 
that each swap has a fixed cost and as in the previous section, we assume that 
the keys are distinct and that all n\ permutations are equally likely to be the 
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input: this in particular implies that the pivot is chosen uniformly at random 
from the array. 

We should specify the partitioning procedure. Assume that we have to sort n 
distinct keys, where their locations in the array are numbered from left to right 
by 1,2 ,,n. Set two pointers i •<— 1 and j <— n — 1 and select the element at 
location n as a pivot. First, compare the element at location 1 with the pivot. 
If this key is less than the pivot, increase i by one until an element greater 
than the pivot is found. If an element greater than the pivot is found, stop and 
compare the element at location n — 1 with the pivot. If this key is greater than 
the pivot, then decrease j by one and compare the next element to the pivot. 
If an element less than the pivot is found, then the j pointer stops its scan and 
the keys that the two pointers refer are exchanged. 

Increase i by one, decrease j by one and in the same manner continue the 
scanning of the array until i > j. At the end of the partitioning operation, 
the pivot is placed in its final position k, where 1 < k < n, and Quicksort is 
recursively invoked to sort the subarray of k — 1 keys less than the pivot and 
the subarray of n — k keys greater than the pivot [29, 30]. 

Note that the probability of a key being greater than the pivot is 


n — k 

n — 1 


The number of keys which are greater than pivot, and were moved during 
partition is 


n — 1 
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Therefore, considering also that pivots are uniformly chosen and noting that 
we have to count the final swap with the pivot at the end of partition operation, 
we obtain 


y- (n — k)(k — 1) n 2 

2— n(n — 1 ) 6 3 ' 

k =1 v ' 

Let S n be the total number of exchanges, when the algorithm is applied to an 
array of n distinct keys. We have that So = Si = 0 and for n > 2, the following 
recurrence holds 

S n = “Number of exchanges during partition routine” + S k _i + S n _ k . 

Since the pivot is chosen uniformly at random, the recurrence for the expected 
number of exchanges is 

9 1 n 

E(Sn) = 7 + 77 + - E + HSn-k)) 

6 3 n ' V / 

fc=i 

E(S' n ) = — + - H— ^E(S'/ C ). 

6 i n 

k= o 

This recurrence relation is similar to the recurrences about the mean number 
of comparisons and will be solved by the same way Subtracting (n— l)E(S' n _i) 
from riE(S n ), the recurrence becomes 

2r> —i— 3 

nE(S n ) - {n - l)E(5' n _ 1 ) = —+ 2E(5' n _ 1 ) 

6 

^ E(S' n ) _ 2 n + 3 E(S' n _i) 

n+1 6 n(n + 1) n 
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Telescoping, the last relation yields 


E (S n ) 

n + 1 


y, 2j + 3 1 

% 6 iU +1) + 3 


E 


30 + 1 ) 


n 


+E 


1 

2j(j + 1) 



1 

3 

1 

3 




Et 


, 3=3 


+ 



1 ^ 1 
n + 1 ) + 3' 


1 

3 


Tidying up, the average number of exchanges in course of the algorithm is 


E (S n ) 


(■n + 1 )H n n 5 
3 9 _ 18‘ 


Its asymptotic value is 

nlog e (n) 

3 ' 

It follows that asymptotically the mean number of exchanges is about 1/6 of 
the average number of comparisons. 

In a variant of the algorithm analysed in [46] and [63], which we briefly men¬ 
tioned in the introduction, partitioning of n keys takes n + 1 comparisons and 
subfiles of m or fewer elements are sorted using insertion sort. Then the av¬ 
erage number of comparisons, partitioning stages and exchanges respectively, 
are 


E(C' n ) — (n + 1)(2 H n+ i — 2H m+ 2 + 1), 

77-I-I 

E(P B ) = 2 • - 1 and 

m + 2 

E(-S'n) = ( n + 1) (^H n+ 1 — —H m+ 2 + — —) + —, for n > m. 

\3 3 6 m + 2/ 2 
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For m = 0, we obtain the average quantities when there is no switch to inser¬ 
tion sort. Note that in this case the expected costs are 

E(C n ) = 2 (n + 1 )H n - 2 n, 

E(P„) = n and 
E(5„) = 


2.3 Variance 

We now similarly discuss the variance of the number of comparisons in Quick¬ 
sort. Although the result has been known for many years - see [46], exercise 
6.2.2-8 for an outline - there is not really a full version of all the details writ¬ 
ten down conveniently that we are aware of, so we have provided such an 
account - this summary has been put on the arXiv, [33]. The sources [43], 
[55] and [63] were useful in putting the argument together. Again, generating 
functions will be used. The result is: 

Theorem 2.3.1. The variance of the number of comparisons of Quicksort on n 
keys, with a pivot chosen uniformly at random is 

Var (C n ) = 7n 2 - 4(n + 1 ) 2 H^ - 2 (n + T)H n + 13n. 

We start with a recurrence for the generating function of C n , namely f n (z ) = 

OO 

^P(Q r) = k)z k . We will use this to reduce the proof of the Theorem to 

k =0 

proving a certain recurrence formula for the expression 
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Theorem 2.3.2. In Random Quicksort of n keys, the generating functions f, 
satisfy 

~n—1 n 

fn{z) = 

' 1 
3 =1 

Proof. We have, using the equation 


C n — C Un _ i + C*_ Un + n — 1 


that 

n 

P(C' n = k) = J2 V(Cn = k\U n = m)- 

m= 1 

n k—(n— 1) 

= E E 

z —' z —' n 

m= 1 j=l 

noting that C ' m _1 and C„,_ m are conditionally independent subject to the pivot. 

Thus 

n fc-(n-l) 

P(a = k)z k = - E E nan-i=J> j P(C'„_ m = A:-(n-l)-j > fc -^- 1 ^^- 1 

n z —' z —' 

ra=l j=l 

Multiplying by and summing over k, so as to get the generating function f n 
of C n on the left, we obtain 


fn{ z ) = ~ ZEE P(C m _ 1 = j)z j P(C n - m = k - (n - 1) - 

fc=l m=l j = l 

n _l n fe-(n-l) n-l+J 

= yE E E P(C'„-m = A:-(n-l)-j> fc - (n - 1) - J ' 

m=l jr = l fc=l 

~n-l n 

= 1 - E fm-l(z)f n -m(z), (2.1) 

m— 1 



27 


as required. ■ 

This of course will give a recursion for the variance, using the well-known 
formula for variance in terms of the generating function f x (z ): 

Var(A') = /"(1) + f' x ( 1) - (A-(l)) 2 . 


We use this formula together with Eq. (2.1). The first order derivative of f n (z) 
is 


( n _ 1 Vn-2 w ^n-1 w 

/,'(-) = —- — - Y fj-l(z)fn-j(z) + Y fj-l( Z )fn-j(z) 

1=1 3=1 


v n—1 


+ Y fj-l( Z )fn-j( Z )- 


3=1 


From standard properties of generating functions, it holds that 


/n(l) — E(C' n ). 


Differentiating again we obtain 


m = (,i - 1)( ” ~ 2) "" ± + (ra - 

l=i 1=1 

(n — l)^ n 2 ^ „ / \ £i / \ ('n — l)^ n 2 ^ , . . , . 

+ ---X! fj-l( Z )fn-j( Z ) + --- Y fj-l( Z )fn-j(z) 


1=1 


1=1 


^n— 1 


-n— 1 


+ V £ + VS /I-.WA-j M 

1=1 1=1 

(n — l)z n_2 — " n_1 


n 


Y fi-lWn-fc) + ^ f'j-lf'n-M) 


1=1 


1=1 


..n—1 


+ X! fj~l( Z )fn-j( Z )- 


1=1 
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Setting z = 1 [55], 

2 n 2 n 

fn ( 1 ) = (n - l)(w - 2) + ~(n - !) M]-\ + ~( n ~ X ) Mn ~i 

j =i j=i 

1 n g " 

+ n XXX-iW + + “ 

3 = 1 j=l 

where are /j_ 1 (l), /^(l), i.e. the mean number of comparisons 

to sort an array of (j — 1) and (n — j) elements respectively. Setting B n = ^^2, 
we obtain 


/i„ = 


n - 1\ 2(n - 1) A . r 2 ^ 1 ^ „ r „ r 

+- X M I ~ 1 + ~ X B i~ l + ~ X 

j= 1 i =1 j=i 


n 


using the symmetry of the sums. What this argument has shown for us is the 
following - compare [43] where it is also shown that this recurrence has to be 
solved, though no details of how to solve it are given. 

Lemma 2.3.3. In order to prove Theorem 2 . 3 . 1 , it is sufficient to show that this 
recurrence is given hy 


B n = 2(n + l) 2 H 2 n - (8 n + 2 ){n + 1 )H n + ?7(23?? 2 + 1?) - 2 (n + 1 

Proof. Combining the equations B n = and Var(A") = f x ( 1) + f x ( 1) — 
(/y( 1)) 2 , the result follows. ■ 

It will be convenient first to develop some theory on various sums involving 
harmonic numbers. Often we used Maple to verify these relations initially, 
but we provide complete proofs here. As a brief reminder, Mj denotes the 
expected number of comparisons needed to sort an array of j keys. Recall 
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that 

Mj — 2(j + l)Hj — Aj. 

Thus, 

n n n n 

E m j-i = E( 2 ^j-‘ - 4 « -1)) = - 4 E<^ - *)• (2 - 2) 

j= i i =1 i =1 i =1 

For the computation of the first sum of Eq. (2.2), a Lemma follows 

Lemma 2.3.4. For n e N 


E^-' = 

j=i 


n(n + l)iF„+i 


n(n + 5) 
4 


Proof. The sum can be written as 


3 =1 


E — 2H-3(l + — j + ... + n/ l+ - + ...+ 


n — 1 


tt /1 0 \ (i 1 + 2 1 + 2 + ... + n — 1 

— -ffn-i(l + 2 + ... + n) — ( 1 H-—-h • • • 4--- 

2 n — 1 


n(n + 1) 
2 

n(n + 1) 


n—1 


fl»-i - E 

i=i 


n=i * 


H n _ 1 - - 
2 2 


1 /n(n + 1) 


1 . 


The last equation can be easily seen to be equivalent with the statement of the 
Lemma. ■ 


Thus we can find out about the sum of the Mj s, that it holds for n e N 
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Corollary 2.3.5. 


_ n e;, n 2 j _ n 

Y Mj -1 = n{n + l)H n+1 -——- 

3 =1 


Proof. Using Lemma 2.3.4 and Eq. (2.2), the proof is immediate. ■ 


n 

Now, we will compute the term E Mj_\M n _ r 

j =i 


the following proof. 

Lemma 2.3.6. For n e N, it holds that 


We shall use three Lemmas for 


El 

3 =i 


4 E - j + l)H n _j 

3 =1 



l)-ldn+l + 



!)■ 


Proof. To do this, we will again use the formula obtained previously for M,, 
We have 


l=i 


!=i 


X] = Y ( (VHj -1 - 4 J + 4 ) (2(n - J + 1) H n —j - 4n + 4j) 


.(" - j + I.)", , - SnYjHj-i + *Yi' 2n 


3 ~ 1 


1=1 


1=1 


l=i 


- 8 E 3 ( n ~ 3 + 1 )#n-j + 16n E 4 _ 16 E 


l=i 

n 


3 = 1 1=1 

n 


E(- _ d + 1) H n —j - 16n 2 + 16 Ed- 

l=i i=i 


n 

We need to work out the value of E 

l=i 



Lemma 2.3.7. For n e N holds 


£iW. 

3 =1 


6 n(n + 1)(2 n + l)H n+ i — n(n + l)(4n + 23) 
36 ' 


Proof. Using the same reasoning as in Lemma 2.3.4, 


E^- I = 2» + s*(i + i) + ... + »»(i + i + ... + 4 1 

3 =1 


— H n — i(l +2 + ... + n ) 

12 + 2 2 l 2 + 2 2 + + (n - l) 2 

- ( 1 H-5-+-- 

2 n — 1 

n(n+l)(2n + l) E7 vVELi* 2 


6 




j=i 


J 


n(n + l)(2n + 1) rr 1 . i , 9 . 

-4-^- -H n _i -4n 3 + 3n 2 - n - 6 , 

6 36 v 


completing the proof. 

n 

We also need to compute ^ j(n — j + l)H n _ r A Lemma follows 

j'=i 


Lemma 2.3.8. For n e N 

n 

^2 jin - 3 + 1 )H n -j = 

3 =1 


6nF[ n+ i(n 2 + 3n + 2) - 5n 3 - 27n 2 - 22 n 


36 
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Proof. We can write j — n + 1 — (n — j + 1). Then, substituting k — n — j + 1 


y - j + = y (( n+_ (■ n ~ j + x )) (■ n _ •? + l ) H - 


n-3 


3 =1 


j=l 


= n 


= n 


1 ) XX" _ •? + “ y ( n _ ■? + 1 ) 2//; 


n-j 


i=i 


i=i 


i) j2 kH k-i- y* 2 ^ 


k- 1- 


fc=l 


fc=l 


These sums can be computed using Lemmas 2.3.4 and 2.3.7. ■ 

n 

In the same manner we shall compute y^(n — j + l)i7 n _j. Changing variables, 


3 =i 


the expression becomes y kH k _ i. Using the previous results, we have 


fc=i 


y, = 4JjT h ( R - j + - 8n 


i=l 


i=i 


n(n + l)if re+ i n(n + 5) 


6n(n + l)(2n + l)i7 n+1 — n(n + l)(4n + 23) 
36 


+ 16n y.) 

i=i 


, / 6nH n+1 (n 2 + 3ra + 2) - 5n 3 - 27n 2 - 22n 
’V 36 

n(n + l)i7 n+ i n(n + 5) 




3 =1 


16n 2 + 16 y j 


i=i 


4 y jHj_i(n -j + l)H n -j - 4n(n + l)(n - l)i7„+i 

3 =1 

4 2 176n 3 176n 

+ -n(n 2 - + -jj- - IT’ 


finishing the proof of Lemma 2.3.6. 
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After some tedious calculations, the recurrence relation becomes 


4 V -j + 1 )H n -j 

j= i ,2 ^ d , -9n 2 + 5n + 4 

- n - + n2^ Bi ~ 1 + -2- 

3 =1 

2 . 9 , 44 . 9 

- g(n - l)^n+i + y - !)■ 


Subtracting n£> n from (n + l)-B n+ i, 


(n + 1 )B n+ i - nB n 


( n+1 n 

^ jHj_i(n - j + 2)H n+1 _j - Y - j + l)#n- : 

i=i i=i 

_ _ , .. —9(n + l) 2 + 5 (n + 1) + 4 -9n 2 + 5n + 4 

+ 25 n + (n + 1)— 4 --2- n --- 

— ( n + l)-((n + l) 2 — l)-ff n +2 + g n ( n2 — ^)H n+ \ 

44 44 

+ ( n + i)y(( n +1) 2 ~ x ) ~ n ^ n2 ~ x ) 


9 


4 ( Y i H 3-^ n ~ j + 2 ) H n+l-j - Y ~ j + l ) H n~; 

S'=1 3 =1 

27n 2 + 17n 2 


+ 2f?„ 


-n(n + l)(n + 2 )H n+2 


2 2 44n(n +1) 

+ -ni/ n +i (n — 1) 4 --• 


We obtain 


(n + l)B n+ i - nB n 


4 ( Y, i^j-i + Yj 
S= i j=i 


4 * 2B n — 2 n{n + + —n(n + 11 ). 
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We have to work out the following sum 


j Hj -1 H n+ 1- j . 
7=1 

We note that 


'y ] jHj-\H n+ i—j — 2/7 1 i + 3//2 //„-2 + ... + (n — 1)77 ,,-2 772 + nH n -iHi 

3 =1 

= (2.3) 

7=1 


Sedgewick [63], presents and proves the following result: 

Lemma 2.3.9. 


£ tfA+i-i = (« + 2 )(H 2 n+1 - H^) - 2(n + 1)(H W , - 1), 


Proof. By the definition of harmonic numbers, we have 

n n 1 

y HiH n+ i~i =y Hi y - 

i= 1 i =1 j=l 

and the above equation becomes 

n . n+1—j n 

y - y ^< = y - ((n + 2 - j)H n+1 _j - (n + 1 - j)) , (2.4) 

i=i *—1 7=1 J 


^ Hj = (n + l)H n — n. 
7=i 


using the identity [44], 
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Eq. (2.4) can be written as 


I i TT I L 

(n + 2 ) ^2 n+1 ~ J ~ Hn + l ~j ~ ( n + l ) Hn + n 

3 = 1 J 3=1 

n rr 

= (n + 2) ^ v±lzi _ (( n +1 )H n - n) - ( n + 1 )H n + n 
3=1 J 

n TT 

— (n + 2) y] -:—- — 2(n + l)(if n+1 — 1). 

U J 

It can be easily verified that 


E 


3=1 




j 


n 


E 


H n ~j 

j 


n 


+E 


1 

j(n + l-j) 


n-l 

_ n—j 
3=1 J 


+ 2 


//» 

n + 1 


Making repeated use of Eq. (2.5), we obtain the identity 


E 

3=1 


tfn+i-i 


2 E 

fc=i 


k + 1 


We have then 


2 yJ^ = 2 y 

Z^ i. _i_ i z^ 


n+1 


fc=l 


A; + 1 


k =2 
»l+l fc 


tffc-1 

k 


n+1 

2 E 

fc=i 


+ ; -i 


n+1 y- tt n+1 

2 E + 2 E 


fc=l 
n+1 n+1 


fc=l 

r( 2 ) 


A; 2 


2 E E+ 2H 'nk = 2 E E+ 

^=1 .7=1 j =1 k=j 


k= 1 j= 
n+1 n+1 


2 EE^- 2 ". 


/c=l j=k 


( 2 ) 
n+1' 


(2.5) 
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The order of summation was interchanged. We can sum on all j and for k = j, 
we must count this term twice. We obtain 




fc=l 


//;, 

k+1 



_ o / t ( 2 ) — h 2 — H k 

z -"n+1 — rL n +1 -"n+l- 


r( 2 ) 


( 2 . 6 ) 


Finally 


E = (« + 2)(^ +1 - - 2(n + l)(ff„ +1 - 1). 1 

2=1 


The following Corollary is a direct consequence of Eq. (2.5) and (2.6). 
Corollary 2.3.10. For n e N, it holds 




n+1 


rr( 2 ) _ 

-"n+1 — 




3 = 1 


J + 1‘ 


Proof. 


ff„ 2 +i - = K - + 2^_ 


and by iteration, 


/7 2 — H ^ — 2 

n n +1 -"n+l — z 


n 

E 


H, 


-J + l 


We will use the above Lemma and Corollary in our analysis. We have that 




n 


i =1 


n+1 —2 i Hi ( H n —i + 

2=1 


E //,+, 




Hi 


n + l-ij I +-f-' n + 1 - i' 

2=1 2=1 


The second sum substituting j = n + 1 — i becomes 


n rr n TT 

E n i _ \ -"n+l-j 

n + 1 — i 

i= 1 j=l 


J 
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As we have seen it is equal to 


H n+ i_j_" H. 


.7 = 1 J 


Hence, by Lemma 2.3.9, 


£ HjHn-i = (n + 2)(77 2 + , - 77®,) - 2(n + l)(Jf„ +1 - 1) - (77 2 + , - 77®,) 
i=i 

= («+ 1)((J7 2 +, - ffSi) - 2(77„ +I - I))- 


Using the above result and Eq. (2.3), we have 


£7/7,-,^.-, = (" + 2 ) (( 7 r- - 77 «) + MUiU 


Having worked out all the expressions involved in the following relation 


l)-Bn+i - nB n = 4f y ZjHj-i + ^.///; |//„ j-i 

\j=i j= 1 


+ 2B n — 2 77(77 + l)iL n _|_i + — 71(71 + 11). 


This becomes 


(77 + l)£ n +i - nB n 


77(77 + l)iT n+1 77(77 + 5) / 77 + 2 


~ + l 2 


((77 2 +1 - 77®,) - 2(77,,+, - 1)) 


+ 2 B n — 277(77 + l)iL n+ i + — 77(77 + 11) 

= 2(n + l)(n + 2)(77 2 +1 - 77®,) - 4(n + l)(n + 2)(77„ + , - 1) - ~ 11 + 2B„. 
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Dividing both sides by (n + 1 ){n + 2) and unwinding the recurrence, 


B„ Bn 




71+1 1 


i= 1 


i= 1 


(i ~ l)(i ~ 2) 

+ Hi + 1 ) 


Hence 


Bn 


n + 1 


2(n + l)(tf 2 - i^ 2) ) + 4n - 4nW n - 4^(n + l)H n - 2n 
^-w'i + 2 3 


t=i 


2i i + 1 


= 2(n + l)(tf 2 - H^) - H n (8n + 2) + ^ - 3 + 


= 2 (n + 1)(H 2 - HW) - H n (8n + 2) + 


2 n + 1 

23n 2 + 17n 
2(n + l) 


Finally, multiplying by (n + 1) we obtain 


= 2(n + If {Hi - Hff - H n (n + l)(8n + 2) + 


23n 2 + 17n 


Consequently, by Lemma 2.3.3 the variance of the number of comparisons of 
randomised Quicksort is 

7n 2 - 4 (n + 1) 2 W, (2) - 2 (n + l)tf n + 13n, 

completing the proof of Theorem 2.3.1. 


2.4 “Divide and Conquer” recurrences 


We have computed the mean and variance of the number of comparisons made 
by Quicksort that mainly contribute to its time complexity. Because of the sim- 
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pie structure of the algorithm (dividing into smaller subproblems) we can in 
fact approach many other related problems in the same spirit. Let F(n) de¬ 
note the expected value of some random variable associated with randomised 
Quicksort and T(n) be the average value of the “toll function”, which is the 
needed cost to divide the problem into two simpler subproblems. Then F(n) 
is equal to the contribution T(n), plus the measures required sort the resulting 
subarrays of (i — 1) and (n — i ) elements, where the pivot i can be any key of 
the array with equal probability. 

Thus, the recurrence relation is 

1 n 

F(n) = T(n ) H— £(f (* - 1) + F(n - 0) 

n i =1 
n 

= T(n) + - £ F(i — 1). 

i= 1 

This is the general type of recurrences arising in the analysis of Quicksort, 
which can be manipulated using the difference method or by generating func¬ 
tions. Since an empty array or an one having a unique key is trivially solved, 
the initial values of the recurrence is 


F(0) = F{ 1) = 0. 


The first method leads to the elimination of the sum, by subtracting (n — 
1 )F(n — 1) from nF(n ) - see [46]. The recurrence becomes 


nF(n) — (n — 1 )F(n — 1) = nT(n) — (n — 1 )T(n — 1) + 2 F(n — 1) 
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and dividing by n(n + 1) we have 

F(n) nT(n ) — (n — 1 )T(n — 1) F(n — 1) 

n+1 n(n +1) n 


This recurrence can be immediately solved by “unfolding” its terms. The gen¬ 
eral solution is 


F(n) = (n+!)['£ ^ ~ fc ^ ~ 1} + m 


= n 



When the sorting of subarrays of m keys or less is done by insertion sort, the 
solution of the recurrence is 


F(n) 



jT(j) ~ {j ~ 1 )T(j - 1) 
j{j + !) 

jT(j) - (j - 1TO -1) 

j{j + 1) 


+ 

+ 


F{m + 1)\ 
m + 2 / 

T(m + 1)\ 
m + 2 / 


since n — 1 > m. 

Another classic approach, which is more transparent and elegant, is the ap¬ 
plication of generating functions. The recurrence is transformed to a differ¬ 
ential equation, which is then solved. The function is written in terms of 
series and the extracted coefficient is the solution. Multiplying by nx n and 
then summing with respect to n, in order to obtain the generating function 
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OO 

G(x) = ^ F(n)x n , we have 

n =0 


y ^nF(ri)x n 


J2nT(n): r" + 2^^F( 


i - l)x n . 


n=0 n=0 n=0 i=l 

The double sum is equal to 


oo n oo / \ 

n=0 2=1 n=l 

and the differential equation is 


xG'(x) = yy nT(n)x n + 

n =0 


2a:G(a;) 
1 — x 


Cancelling out x and multiplying by (1 — x) 2 , 


G"(a;)(l — x) 2 = (1 — x) 2 nT(n)x n 1 + 2(1 — x)G(x) 

n= 1 
oo 

(G(x)(l — x) 2 )' = (1 — x) 2 YnTWx^ 1 

n =1 

/ oo 

(1 — x) 2 y^wT(n)x n ~ 1 d x + C 

n =i 

/ OO 

(1 — x) 2 Y riT(n)x n - ] dx + C 

n =i 


G(x) = 


(l-xf 


where C is constant, which can be found using the initial condition G(0) = 0. 
The solution then is being written as power series and the coefficient sequence 
found is the expected sought cost. 

Now, one can obtain any expected cost of the algorithm, just by using these 
results. The “toll function” will be different for each case. Plugging in the 
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average value, the finding becomes a matter of simple operations. This type of 
analysis unifies the recurrences of Quicksort into a single one and provides an 
intuitive insight of the algorithm. 


2.5 Higher moments 

We have effectively calculated the first and second moments of C n in Quicksort. 
Existing literature does not seem to address much questions about skewness 
and kurtosis, which are sometimes held to be interesting features of a random 
variable. Here, we present an inconclusive discussion about the sign of the 
skewness. 

Using the probability generating function, we can obtain higher moments of 
the algorithm’s complexity. A Lemma follows 

Lemma 2.5.1. Let a random variable X having probability generating function: 

OO 

fx(z) = J2n x = n)z n . 

71=0 

For the k-th order derivative it holds that 

= E(X • (X - 1) • ... • (X - k + 1)). 

Z= 1 

Proof. Simply by computing successively the k-th order derivative of fx(z), 
we obtain 


d k fx{z) 

dz k 


d k fx{z ) 

dz k 


n ■ (n — 1) • ... • (n — k + 1)P(X = n) ■ z n k 

n=0 
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Setting z — 1, the proof follows directly. Note that the argument is similar to 
continuous random variables. ■ 


Using Maple, we obtained a recursive form for the general A;-th order deriva¬ 
tive of the generating function. 


Lemma 2.5.2. Let 


~n —1 n 

fn(z) = - -5Z/j-l(2)/n-j(*) 

3=1 


be the generating function of Quicksort’s complexity in sorting n keys. The k e N 
order derivative is given by 


d k fn(z) = 1. _ 
d z k n 



r(n) 

T(n — i) 



where the V function is defined for complex variable z f 0, —1, —2,... as 


r(*) := / p~ 'e 




dp 


and when z is a positive integer, then r(z) = (z — 1)!. 
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Proof. For k — 0, the result follows trivially. Assume that the statement of the 
Lemma holds for k = m. The (m + l)-th order derivative is 


d m+1 f n (z) = lA / m\ T(n) 
d z m+1 n \i ) T(n — i) 


(n — i — 1 )z 


n—i —2 


dz r 




3 =1 




m+1 


,-i-i d 


m—i +1 


d-m-i+i 


fj-l(z) ■ fn-j(z) 

3 =1 


I \ - / m 
n +' \/ — : 


r(n) 


1 / T(n - i + 1) 


\n — i\z 


i-i-i d 


m— 2+1 


d~m—i+i 


5^/j-l(-) ' fn-j(z) 

3=1 


£ 

t=0 v 7 


"A r(n) d 


m—i+1 


r(n - i) 


5^/i-iW ' +(+• 

j'=i 


Note that 

r(n) , _ r(n) 
r(ii-j + i) 1 r(n-<)' 


Therefore, 


d m+1 f n (z) 

d z m+1 



£H 

T (n — i) 




X 


qra—2+1 


d^m-i+i 


n X 

' /n-j(*)V 

3 = 1 ' 


The well-known identity, (e.g. see [26]) 



+ 



I'm. + 1\ 

v i r 


completes the proof. ■ 

We should point out that Lemma 2.5.2 is an immediate consequence of Leib¬ 
niz’s product rule. Next, we shall ask a Question about the sign of the skewness 
of the time complexity of the algorithm, as it is moderately difficult to solve the 
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recurrence involved, in order to compute the third moment. We already have 
seen that the possibility of worst-case performance of the algorithm is rather 
small and in the majority of cases the running time is close to the average time 
complexity which is O (n log 2 (n)). Intuitively, this suggests that the complexity 
is negatively skewed. We present the following Question: 

Question 2.5.3. Is the skewness S(C n ) of the random number of key comparisons 
of Quicksort for the sorting ofn> 3 keys negative? 

Note that the cases n — 1,2 are deterministic, since we always make 0 and 1 
comparisons for the sorting. This Question may have an affirmative answer, 
which can be possibly proven by an induction argument on the number of keys. 
However, great deal of attention must be exercised to the fact that the random 
number of comparisons required to sort the segment of keys less than the 
pivot and the segment of keys that are greater than the pivot are conditionally 
independent, subject to the choice of pivot. 


2.6 Asymptotic analysis 

After having examined the number of comparisons of Quicksort, in terms of 
average and worst case scenarios, and its variance, it is desirable also to study 
the concentration of the random variable C n about its mean. One might hope, 
by analogy with other probabilistic situations, that for large values of n the 
number of comparisons is highly likely to be very close to the mean. 

The analysis will be confined to the number of comparisons, because this is the 
most relevant measure for this thesis. Since our results will be asymptotic in 
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nature, we need to have some relevant ideas about convergence of sequences 
of random variables. The following definitions come from [8], [21]. 

Definition 2.6.1. 

(i) A sequence of random variables {Xi,X 2 ,.. .}is said to converge in distribution 
(or converge weakly) to a random variable X if and only if 

lim F n (x) = F(x) 

n—» oo 


at every point x where F is continuous. Here F n (x) and F(x) are respectively 
the cumulative distribution functions of random variables X n and X. We shall 
denote this type of convergence by X n —» X. 

(ii) A sequence of random variables {X x , X 2 ,...} is said to converge in probability 
to a random variable X if and only if\/e> 0 holds 

lim P(|X n - X\ > e) = 0. 

71—» OO 

'P 

We will denote this type of convergence by X n —> X. 

(iii) A sequence {Ai, X 2 ,. . .} of random variables is said to converge in L p -norm 
to a random variable X if and only if 

lim E(|X n - X\ p ) = 0. 

71—)■ OO 

Note that convergence in L p -norm, for p > 1, implies convergence in proba¬ 
bility, and it is also easy to see that convergence in probability implies con¬ 
vergence in distribution: see e.g. [8]. Both converse statements are false in 
general. 
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We also present the definition of martingale, which shall be employed in a 
later stage of our analysis. 

Definition 2.6.2. Let {Z\ , Z 2 ,.. .} be a sequence of random variables. We say Z n 
is a martingale if and only if 

(i) E(|Z n |) < oo, Vn e N. 

(ii) K(Z n \Z n _i, Z n _ 2 , • • •, Z\) = Z n _ i. 

2.6.1 Binary trees 

In this subsection, a central notion to the analysis of Quicksort, and in general 
to the analysis of algorithms is discussed. We begin with a definition. 

Definition 2.6.3. A graph is defined as an ordered pair of two sets (V, E). The 
set V corresponds to the set of vertices or nodes. The set E is the set of edges, 
which are pairs of distinct vertices. 

One kind of graph we concentrate on are trees. A definition of a tree (which 
suits us) is as follows [44]: 

Definition 2.6.4. Tree is a finite set A of nodes, such that: 

(i) There is a unique node called the root of the tree. 

(ii) The remaining nodes are partitioned into k e N disjoint sets A 1? A 2 , ... ,A k 
and each of these sets is a tree. Those trees are called the subtrees of the root. 

A particularly common tree is a binary tree. It is defined as a tree with the 
property that every node has at most 2 subtrees, i.e. each node - excluding 
the root - is adjacent to one parent, so to speak, and up to two offspring, 
namely left and right child nodes. Note here that nodes which do not have 
any child nodes are called external. Otherwise, they are called internal. The 
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size of a binary tree is the number of its nodes. The depth of a node is simply 
the number of edges from that node to the root in the (unique) shortest path 
between them and the height of a binary tree is the length from the deepest 
node to the root. 

An extended binary tree is a binary tree with the property that every internal 
node has exactly two offspring [44]. Let Dj be the depth of insertion of a key 
in a random binary tree of size j — 1. Since the root node is first inserted to an 
empty tree, it holds that D l := 0. The next inserted key is compared with the 
key at the root and if it is smaller, is placed to the left; otherwise is attached as 
a right subtree, thus D 2 := 1. 

Thus, we define the internal path length of an extended binary tree having n 
internal nodes, to be equal to the sum over all internal nodes, of their distances 
to the root. Let us denote this quantity by X n . We then have that 

n 

X n = Y J Dj- 

3 =1 

Similarly, the external path length is defined as the total number of edges in 
all the shortest paths from external nodes to the root node. The next Lemma 
gives a relationship between those two quantities. 

Lemma 2.6.5 (Knuth [44]). For an extended binary tree with n internal nodes 
it holds that 

Y n = X n + 2 n, 


where Y n denotes the external path length of the tree. 
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Proof. Suppose that we remove the two offspring of an internal node v, with 
its offspring being external nodes of the tree. We suppose that v is at distance 
h from the root. Then the external path length is reduced by 2(h+l), as its two 
offspring are removed. At the same time, the external path length is increased 
by h, because the vertex v has just become an external node. Thus, the net 
change is equal to —2(h+l)+h = — (h+2). For internal path length the change 
is a reduction by — h as v is no longer internal. Thus, overall, the change in 
Y n — X n is equal to — (h + 2) — (—h) = —2. The Lemma follows by induction. ■ 

Binary trees play a fundamental and crucial role in computing and in the anal¬ 
ysis of algorithms. They are widely used as data structures, because fast inser¬ 
tion, deletion and searching for a given record can be achieved. In Quicksort, 
letting nodes represent keys, the algorithm’s operation can be depicted as a 
binary tree. The root node stores the initial pivot element. Since the algorithm 
at each recursion splits the initial array into two subarrays and so on, we have 
an ordered binary tree. The left child of the root stores the pivot chosen to sort 
all keys less than the value of root and the right child node stores the pivot for 
sorting the elements greater than the root. 

The process continues until the algorithm divides the array into trivial subar¬ 
rays of order 0 or 1, which do not need any more sorting. These elements are 
stored as external or leaf nodes in this binary structure. It easily follows that 
for any given node storing a key k, its left subtree store keys less than k and 
similarly its right subtree contain keys greater than k. 

In variants of Quicksort, where many pivots are utilised to partitioning pro¬ 
cess, the generalisation of binary trees provides a framework for the analysis, 
though we do not develop this in detail here. In the next subsection, the lim- 
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iting distribution of the number of comparisons will be analysed, in terms of 
trees. 

2.6.2 Limiting behaviour 

We have previously seen that the operation of the algorithm can effectively 
be represented as a binary tree. Internal nodes store pivots selected at each 
recursion step of the algorithm: so the root vertex, for example, stores the first 
pivot with which all other elements are compared. We are interested in the 
total number of comparisons made. To understand this, we note that every 
vertex is compared with the first pivot which is at level 0 of the tree. The array 
is divided into two parts - those above the pivot, and those below it. (Either 
of these subarrays may be empty). If a subarray is not empty, a pivot is found 
in it and is attached to the root as a child on the left, for the elements smaller 
than the first pivot or on the right for the elements larger than the pivot. 

The process then continues recursively and each element at level k in the tree 
is compared with each of the k elements above it. Thus the total number of 
comparisons made is the sum of the depths of all nodes in the tree. This is 
equivalent to the internal path length of the extended binary tree. Thus, 

Y n = C n + 2 n, 

where Y n is the external path length of the tree. To see this, we simply use 
Lemma 2.6.5. This fact can be also found in [49]. 

Generally, assume that we have to sort an array of n distinct items with pivots 
uniformly chosen. All n\ orderings of keys are equally likely. This is equivalent 
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to carrying out n successive insertions [46], Initially, the root node is inserted. 
The second key to be inserted is compared with the key at the root. If it is less 
than that key, it is attached as its left subtree. Otherwise, it is inserted as the 
right subtree. This process continues recursively by a series of comparisons 
of keys, until all n keys have been inserted. Traversing the binary search tree 
in order, i.e. visiting the nodes of the left subtree, the root and the nodes of 
the right subtree, keys are printed in ascending order. Thus, Quicksort can be 
explicitly analysed in this way. 

Recall that each internal node corresponds to the pivot at a given recursion of 
the Quicksort process (e.g. the depth of a given node). Thus the first pivot 
corresponds to the root node, its descendants or child nodes are pivots chosen 
from the two subarrays to be sorted, etc. Eventually, after n insertions, we 
have built a binary search tree from top to bottom. We can use this approach 
to understand the following Theorem of Regnier [58]. 

Theorem 2.6.6 (Regnier [58]). Let random variables Y n and X n denote re¬ 
spectively the external and internal path length of binary search tree, built by n 
successive insertions of keys. Then the random variables 

_ Y n — 2(n + l)(H n+ i — 1) _ X n — 2 (n + 1 )H n + 4 n 
71+1 71+1 


form a martingale with null expectations. For their variances it holds 


Vaii z n) = 7 - 2 X- 2 XiM + 0 (t h 

3 n \n 


Proof. By induction on n, the base case being trivial. Suppose we have an 
(n — 2)-vertex tree and consider the insertion of the (n — l)-th key in the 
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random binary tree. Its depth of insertion is D n _i, so that a formerly external 
node becomes an internal and we see that its two (new) descendants are both 
at depth _D n _i + 1. Recall that D n is the random variable counting the depth 
of insertion of a key in a random binary tree of size n — 1, thus Di 0 and 
D -2 := 1. The following equation concerning conditional expectations holds: 

nE(D n \Di, ..., D n _ i) = ( n — l)E(D n _ 1 |Zi) 1 ,..., D n _ 2 ) — D n _i + 2(D n _i + 1). 

This recurrence yields 

nE(D n |T) 1 ,..., Z? n _i) — ( n — l)E(D ra _ 1 |T) 1 ,..., D n _ 2 ) = -D n - i + 2. 


Summing and using that the left-hand side is a telescopic sum, 


n— 1 


nE(£> n |£>i,..., D n _i) — + 2) — T) 


n—1 • 


2=1 


The last equation is justified by Lemma 2.6.5. Also, we have 


71 1 

, D n _ i) = E(lj ) ,_ 1 + 2 + D n \Di ,..., D n _ i) =-Tn.-i + 2. 

n 


Thus, taking expectations and using E(E(C/|V")) = E(C7), 


E(y„) = —Efy.^) + 2 ^ 5L1 = 3AA + 

n n + 1 n n + 1 


This recurrence has solution 


E(y n ) = 2(n+l)(H n+1 -l). 
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For Z n , we deduce that 

E(Z n | D u D n ) — E ( Yn ~^ n ) d D \ 

\ n + 1 / 

2 | y ra _! E(y ra ) 

n + 1 n n + 1 

— ^n-l- 

Therefore, Z n form a martingale. Further, note that Z n is a linear transforma¬ 
tion of the internal path length X n , so we get that 

Var (Z n ) = . * • Var(X n ). 

yn + l) 2 

Previously, we saw that the number of comparisons is just the internal path 
length of a binary search tree. As a reminder the variance of the number of 
comparisons is equal to 

Var(C n ) = 7 n 2 - 4(n + 1) 2 #, ( , 2) - 2(n + l)if n + 13n. 


Hence, 


Var(Z n ) = 7 


71+1 


lg (2) ^ 

n 77+1 (77+ l) 2 ' 


For the purpose of obtaining the asymptotics of the variance, an important 
family of functions, which is called polygamma functions are discussed. The 
digamma function is 

fp( z ) = T-!oge r (-) 
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and for complex variable z ^ 0, — 1 , —2, ... can be written as [1], 


^(z) = -7 + Y1 


3 + 1 j + z 


(2.7) 


In general, \/k e N, the set 


A k-\-l 

= ___ iog e r (^) ? 


with ^°\z) = ip(z), forms the family of polygamma functions. Differentiating 
Eq. (2.7), 




(j + A' 2 


( 2 . 8 ) 


By Eq. (2.7), it easily follows that 


H n = ip(n + 1) + 7. 


Further, using the fact that £(2) = lim^,*, = — [1], where £(s) = 

1 ® 

Y^jLi — is Riemann zeta function for 9fe(s) > 1, we obtain 


H m + = £ _ v,a» ( „ +1) + 2 W» + i) + 7) 

n +16 n +1 


(2.9) 


Eq. (2.9) is asymptotically equivalent to 


21og e (w) 
6 n 
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thus, the asymptotic variance is 


Var(Z„) = 7 _ + 0 (T 

3 n \n 


= 7 


2vr 2 


O 


loge(n) 


n 


Remark 2.6.7. We note that there is a typo in the expression for the asymptotic 
variance in Regnier’s paper [58], which is given as 


7-^ + 0( i|. 


The correct formula for the asymptotic variance is stated in Fill and Janson [22], 
The key point about martingales is that they converge. 

Theorem 2.6.8 (Feller [21]). Let Z n be a martingale, and suppose further that 
there is a constant C such that E (Z 2 ) < C for all n. Then there is a random 
variable Z to which Z n converges, with probability 1. Further, E (Z n ) = E (Z) for 
all n. 

We showed that 

Z n A z. 

It is important to emphasise that the random variable Z to which we get con¬ 
vergence is not normally distributed. We saw that the total number of com¬ 
parisons to sort an array of n > 2 keys, when the pivot is a uniform random 
variable on (1,2is equal to the number of comparisons to sort the 
subarray of U n — 1 keys below pivot plus the number of comparisons to sort 
the subarray of n — U n elements above pivot plus n — 1 comparisons done to 
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partition the array. Therefore, 


X n — X Un _i + x^_ v + n - 1, 


where the random variables Xu n -i and X*_ Un are identically distributed and 
independent conditional on U n . 

Consider the random variables 

v , X n -E(X n ) 

* n 

n 


The previous equation can be rewritten in the following form 


Y n = 


x u n -1 + X n-u n + n-l- E(X n ) 


n 


By a simple manipulation, it follows that [22], [59] 


Yn = Y Un _i ■ l Xl + Y' ■ —^ + C n (U n ), 


n 


n 


where 

77 — 1 1 

c„(j) = — + - (E(A-,-_.) + nx'n-j) - E(A'„)) . 


The random variable U n /n converges to a uniformly distributed variable 5 on 
[0,1]. A Lemma follows 


Lemma 2.6.9. Let U n be a uniformly distributed random variable on { 1,2, , n}. 
Then 


U n 



n 
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where 5 is uniformly distributed on [0,1]. 


Proof. The moment generating function of U n is given by 


EE 1 EE i p s(n+ 1 ) _ s 

= E p w> = = - ■ E'' = - ' • 

k =1 k =1 


For the random variable U n /n, it is 


i EE 1 „s(ra+l)/n _ pS/n 

Mu./M = Mu„(s/n) = - ■ E e* 1 ’ 7 ” = - • 


n ^ n e 

k =1 


*/n _ i 


The random variable 5 has moment generating function 


-l 


e s — 1 


M*(s) = / e st dt = 

Jo 5 


Now the moment generating function of U n /n is an approximation to the av¬ 
erage value of e sx over the interval [0,1] and so, as n tends to infinity, we can 
replace it by its integral 


e sx dx = 


e s -1 


and now all that is required has been proved. 


For the function 


77 — 1 1 

C„U) = — + - ■ (E(AV,) + E(x;_j) - E(X„)) 
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using the previous Lemma and recalling that asymptotically the expected com¬ 
plexity of Quicksort is 2n log e (n) it follows that 

lim C n (n • Ujn) = lim f— - + - • (E(AV n _i) + EQA^J - E(A n ))) 

n—>oo n—>■ oo V Tl Tl ) 

= lim (- —- + -■ (2 - \ \ log (U n - 1) 

ih 00 \ n n \ \ n J 


= 1 + 2£ log e e + 2(1 - 0 log e (l - 0 = C(E), V e e [0,1]. 


n 


n ■ U r 


n 


log e (n - U n ) - 2nlog e (n) 


Thus C n (n ■ U n /n) converges to (7(H), (see as well in [59]). Therefore, we 
obtain 


Jzf(Y) = Jgf(y ■ H + y* • (1 - H) + (7(H)). 

But this does not work for the normal. Indeed, if F (and hence Y*) are normals 
with mean zero and variance a 2 , a necessary condition for YH + Y*(l — H) + 
(7(H) to have mean 0 (which would be needed for the equality to hold) would 
be that (7(H) = 0. This equality happens with probability equal to 0, as we can 
easily deduce. Thus, the distribution which the sequence Y n converges is not 
Gaussian. 

2.6.3 Deviations of Quicksort 

It is desirable to derive bounds on the deviation of the random number of 
pairwise comparisons, needed to sort an array of n distinct elements from its 
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expected value, as n gets arbitrarily large. We remind ourselves that the pivot 
is uniformly chosen at random. 

We shall derive bounds on the following probability 


P 




for e > 0, sufficiently small. From Chebyshev’s inequality, we obtain a bound 
for the above probability. Chebyshev’s inequality is as follows [21]: 

Theorem 2.6.10 (Chebyshev’s inequality). Let X be a random variable, with 
E(A" 2 ) < oo. Then, for any real number a > 0 


P(|AT| > a) < 


E(X 2 ) 


IfK(X) = m and Var(AT) = a 2 , then 


P(|2f — m\ > a) < 

CL 

c 

The random variable Y n = " has mean and variance, 

= K (s^o)=i(kj ■ = 1 

Var(K„) = Var (^y) = • Var(C„) 

— 4(n + 1 ) 2 H { n ] - 2 (n + 1 )H n + (7 n + 13 )n 
(2(n + l)H n — 4n) 2 
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For e > 0, we have 


P(|r„- 1| > e) < 



—4 (n + 1 ) 2 H^ ] - 2 (n + 1 )H n + (7 n + 13 )n 
e 2 (2 (n + 1 )H n — 4 n) 2 


It holds that, 


Var(C n ) = —4(n + l) 2 H& - 2(n + 1 )H n + (7n + 13 )n 

< 7n 2 + 13n 

< 20n 2 , 


using that the other terms are negative and n > 1. Further, using that 


lim H ( ^ = ^ 

n—> oo ‘ 

n= 1 




7T 


2 


6 ’ 


we get 


Var(C n ) = —4(n + l) 2 i/, (2) - 2(n + 1)/F n + (7 n + 13)n 

> (7 n + 13)n — 4 (n + 1) 2 tt 2 /6 — 2(n + 1) (log e (n) + 7 + o(l)) 
= (7-2vr 2 /3)n 2 (l + o(l)). 


From above inequalities, we deduce that Var(C' n ) = @(n 2 ). Thus, 


P(|y.-l|> e )<^ 


@(n 2 ' 1 


! (2 (n + l)i7 n — 4 n)‘ 


= 0((elog e (n)) 2 ) 


Rosier [59] derived a sharper bound. He showed that this probability is 
0(n~ k ), for fixed k. McDiarmid and Hayward [50] have further sharpened 
the bound. Their Theorem is 
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Theorem 2.6.11 (McDiarmid and Hayward [50], McDiarmid [51]). 

Let e = e(n) satisfy -—< e < 1. Then as n —> oo, 

log e (n) 

pQy^ _ X| > g) = n - 2e ( lo gi 2) ( n )- lo ge( 1 / e )+°( Io gi 3) ( n )) 

where logl fc+1) (n) := logy (log^fYt)) and log^ n = log e n. 

In the recent paper of McDiarmid [51], Theorem 2.6.11 is revisited using con¬ 
centration arguments. By this result, it can be easily deduced, that Quicksort 
with good pivot choices performs well, with negligible perturbations from its 
expected number of comparisons. As we previously saw, Quicksort can be de¬ 
picted as an ordered binary tree. The root node or node 1 corresponds to the 
input array of length L\ = n to be sorted. 

A pivot with rank U n = {1,2,..., n} is selected uniformly at random and the 
initial array is divided into two subarrays, one with elements less than the 
pivot and a second, with elements greater than the pivot. Then, the node at 
the left corresponds to the subarray of keys that are less than the pivot, with 
length L -2 = U n — 1 and the node at the right corresponds to the subarray of 
keys that are greater than the pivot, with length L 3 = n — U n . 

Recursively, Quicksort runs on these two subarrays and split them in four 
subarrays, until we get trivial subarrays and the initial array is sorted. For 
j = 1,2 ... let Lj be the length of the array to be sorted at j node and M)” be 
the maximum cardinality of the 2 k subarrays, after k recursions of Quicksort. 
It is [50], 

= ma x{L 2 k +i : i — 0,1,..., 2 A — 1}. 
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The following Lemma gives an upper bound for the probability that the maxi¬ 
mum length of 2 k subarrays M]!, will exceed a times the initial set’s length n. 
We easily see that the upper bound is rather small quantity Thus, we deduce 
that the length of the array is rapidly decreasing, as Quicksort runs. 


Lemma 2.6.12 (McDiarmid and Hayward [50], McDiarmid [51]). 
For any 0 < a < 1 and any integer k > log e (it holds 


W 


P (M£ > an) < a 


( 2e • log e (l/a) 

V k 


k 


So far, the analysis explicitly assumed the presence of distinct numbers. How¬ 
ever, in many sorting problems, one can come across with duplicates. As we 
will see next, we have to consider the presence of equal numbers and how this 
affects the algorithm’s performance. 


2.7 Quicksorting arrays with repeated elements 

In this section, we consider the presence of equal keys. In some cases, there 
may be multiple occurrences of some of the keys - so that there is a multiset 
of keys to be sorted. If a key a* appears .s t times, we call s l the multiplicity of 

O'i- 

Thus, we have a multiset of the array {,si • ai, s 2 • a 2 ,..., s n ■ a n }, with s\ + s 2 + 
... + s n — N and n is the number of distinct keys. Without loss of generality, 
we assume that we have to sort a random permutation of the array {si-l,s 2 - 
2,..., s n ■ n}. Obviously, when si — s 2 — ... = s n — 1, then the array contains 


n — N distinct keys. 
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Consider an array consisting of keys with “large” multiplicities. The usual 
Quicksort algorithm is quite likely to perform badly, since keys equal to the 
pivot are not being exchanged. Recall that Hoare’s partitioning routine [30] 
described in the previous Chapter, utilises two pointers that scan for keys 
greater than and less than the pivot, so at the end of partition, keys equal to 
pivot may be on either or both subarrays, leading to unbalanced partitioning 
which we saw is usually undesirable. The algorithm can be further tweaked 
for efficient sorting of duplicates, using a ternary partition scheme [7]. Keys 
less than pivot are on left, keys equal to pivot on middle and on right, keys 
greater than the pivot. 

Hoare’s partitioning scheme can be easily modified, in order for the pointers to 
stop on equal keys with the pivot. In other words, the lower pointer searches 
for a key greater than or equal to the pivot and the upper one for a key smaller 
than or equal to the pivot. Sedgewick [64] considers a partitioning routine, 
where only one of the pointers stops on keys either greater or smaller than the 
pivot. 

The recurrence for the expected number of key comparisons E(C(si, s 2 , ■ ■ ■ , s n )) 
is 

1 n / 

E(C(si, s 2 , • • •, s n )) = N — 1 + — y] Si ( E(C(si,..., Sj-i)) + E(C(sj+i,..., s n )) 

V 4=1 ' 

The analysis has been done in [64]. The solution of the recurrence is 


E(C(si, s 2 , • • ■, S n )) = 2 n. + i j NH n - 3N - n. 


( 2 . 10 ) 
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Note that when n = N, i.e. the keys to be sorted are distinct, Eq. (2.10) yields 
the expected number of comparisons, when the array is partitioned using n — 
1 comparisons. As we saw in the previous Chapter, a different partitioning 
scheme proposed and analysed in [46], [63], utilises N + 1 comparisons for 
partitioning an array of N keys. In this case, an upper bound for large n to the 
expected number of comparisons is [64], 

E(C(si, S 2 , • • •, s n )) = 2N(H]\r + 1) — 2 + O 




Chapter 3 


Sorting by sampling 


The preceding analysis has shown that Quicksort is prone to poor performance, 
when the chosen pivots happen to be (close to) the smallest or greatest keys 
in the array to be sorted. The partitioning yields trivial subarrays, leading to 
quadratic running time taken by the algorithm, and increasing the chances 
of ‘crashing’, i.e. an application of Quicksort may terminate unsuccessfully 
through running out of memory. On the other hand, good choices of pivot yield 
a much more efficient algorithm. The uniform model suggests that in fact any 
key has equal probability to be selected as pivot. In this Chapter, we discuss 
and analyse the general idea of how one can increase the probability that 
the selected pivots will produce (reasonably) balanced partitions, by trying to 
ensure that their ranks are ‘near’ the middle of the array. 

A naive idea towards this, would be finding the median of the keys and use this 
as pivot. However, it is obvious that the finding of median imposes additional 
costs to the algorithm. Therefore, despite the always good choices, this method 
might not be better than choosing the pivot randomly. Instead of finding the 
median of the array to be sorted, it might be more efficient to randomly pick a 
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sample from the array, find its median and use this as pivot. In what it follows, 
we will analyse this idea and its variants. 


3.1 Median of (2k+l) partitioning 


Singleton [68] suggested to randomly select three keys from the array to be 
sorted and use their median as pivot, leading to a better estimate of the median 
of the array and reducing further the chances of worst case occurrence. 

This modification can be generalised as to choosing a sample of (2/c + l) keys at 
every recursive stage, computing their median and using that median as pivot, 
partitioning n > 2k + 1 keys. Arrays that contain at most (2 k + 1) keys are 
sorted by a simpler algorithm, such as insertion sort. The cost of sorting these 
small arrays is linear with respect to n. 

Letting C n { 2 k+i\ denote the number of comparisons required to sort n keys 
when the pivot is the median of a random sample of (2 k + 1) elements, uni¬ 
formly selected from the relevant array, the recurrence for the average number 
of comparisons is given by (see [46]) 


E(C n { 2 fc+i}) — n + 1 + 


n 

2k + 1 


E 

3 =1 


3 ~ 1 
k 


n - J 


E , (Cj_i{2fc+i})- 


Multiplying both sides by 


n 


2k + i r 


2k + ij E (^ n{2fe+1} ) Lk + l ](U 


+ !) + 2 E 


3 = 1 




k 


k J' 
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Multiplying by z n and summing over n, in order to obtain the generating func¬ 
tion for the expected number of comparisons, f(z) = Y^=o E(C n { 2 /e+i})^ n 

£ (» + i) E(c »<“«> )2 ” = 

n=0 v ' 

t ( 2 ; + C * 0 C ; 3 ) 

n =0 v 7 n =0 j=l v 7 v 7 

(3.1) 

It holds 


2k + 1 


E(C'n{2fc+i})- n = (2k+l)\ 2k)E(C n{2k+1 })z r 

' ' n —0 

^,2/c+ly*(2fc-|-l) 

{2k + 1)! ’ 


where f( 2k+1 \z) denotes the (2/c + l)-th order derivative of f(z). For the first 
sum in the right-hand side of Eq. (3.1) 



(n + 1 )z n 



The (2 k + 2)-th order derivative of (z/( 1 — z)) can be easily seen by induction 
that is equal to 

(2k + 2)! 

(1 — z) 2k+3 ' 
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Expanding out the double sum of Eq. (3.1), 

0 - © C>(Co ( , +1 >); + 

(© ( 0 E(c ° |2t+i|)+ (0 Q e(c,,2 ‘ +i, )) z2+ 

(© (0 e(c ° (22+i|)+ Q) (0 e(Ci{2 ‘ +i,)+ © © E(c2(2t+i|) ) 23+ ■ ■ ■ 

= ( © E(Gj(2t+i>) + © 1E ( C ' 1 i 2 ‘+ 1 >)+ • • j ( © z + © * 2 + • • j 

-(e(;)bw '•)(£(:)-) 

z k f { - k \z) z k+1 \ {k) 

~ k\ w 2 ) 

v n=0 7 

_ z 2k+1 f^ k \z) 
k\( 1 — z) k+1 ’ 

The recurrence is transformed to the following differential equation 

z 2 k+ij( 2 k+i)^ _ (2/c + 2)z 2k+1 2z 2k+l f^ k \z) 

(2k + l)\ (1 — z) 2k+3 k\(l — z) k+l ' 

- z\ 2k+1 

Multiplying both sides by ( -j , 

(1 - z) 2k+l f (2k+l \z) 2{k + 1) 2(1 -z) k f (k \z) 

(2k + 1)! “ (l-© 2 + k\ ’ 

which is a Cauchy-Euler differential equation. This type of differential equa¬ 
tions arises naturally to the analysis of searching-sorting algorithms and urn 
models [14]. Substituting x = 1 — z, and putting h(x) = /(I — x), we get 

(—l) 2k+1 x 2k+1 h^ 2k+1 \x) 2(k + l) 2(—1 ) k x k hS k \x) 

(2fc + l)! x 2 k\ 



69 


We use the differential operator 0 for the solution of the differential equation. 
It is defined by 


©(h(x)) xh'(x) 


and by induction 



x k lf k \x) 

k\ 


Applying the operator, our equation becomes 


V 2k+1 (Q)h(x) = 2{k+ 2 1 \ 
oc 

where the indicial polynomial is equal to 

*We) = (- 1 ) 2 ‘ +1 ( 2i .® j) - 2(-i)* 



We proceed to identify the nature of the roots of the polynomial. A Lemma 
follows: 

Lemma 3.1.1. The indicial polynomial V- 2 k+i( 0) has 2k + 1 simple roots, with 
real parts greater than or equal to —2. The real roots are 0,1 ,k — 1 ; — 2 ; 
3 k + 2, if k is odd and the 2 [|J complex roots pi,... , p^j with their conjugates 

Proof. Let a = x + iy being a root of the polynomial. Then 


(3.2) 
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From Eq. (3.2), we deduce that k real roots of the polynomial 7 ? 2fc+i(©) are 
0,1,..., k — 1. For a ^ 0,1,..., k — 1, we have that 


(a-k) (a-k-1) (a — 2k) _ 

(* + !)' (A:+ 2) (2/c +1) 1 J 


Suppose that 1Hc(q;) < —2. Then 


a — k 
k + 1 


VQk - k) 2 + y 2 
k + 1 


\/( — (^ + 2)) 2 + y 2 > k + 2 
k + 1 — k + 1 


a — /c — 1 
k + 2 


(x — k — l) 2 + y 2 
k + 2 


\J{~(k + 3)) 2 + y 2 ^ k + 3 
k + 2 ~ k + 2' 


a-2k \/(x — 2k) 2 + y 2 ^ ^{~ 2 ( k + l)) 2 + y 2 ^ 2k + 2 

2k+ 1 ~ 2k + 1 > 2k + 1 - 2k + 1’ 


so the product of the moduli in Eq. (3.3) is greater than or equal to the 
telescoping product (2k + 2)/(k + 1) = 2, arriving in contradiction. Thus, 
for any root a of the polynomial, it is proved that 9fc(a) > —2. Moreover, 
this argument shows that —2 is the unique root with real part equal to —2. 
To see all roots are simple, assume that there does exist a repeated root a. 
Differentiating the polynomial, 


74 + 1 ( 9 ) 




1 

0^7 


and 




1 

a- j' 


(3.4) 
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Eq. (3.2) and (3.4) imply that for a repeated root must hold 


2k 


£ 


1 


a-j 


k -1 


£ 


1 


a-j 


or 

2k 

T -: = 0 . ( 3 . 5 ) 

a — i 

j=k J 

Consequently, Eq. (3.5) implies that 3m(o;) = 0 and a £ (k, k + 1) U (k + 1, k + 
2) U ... U (2k — 1, 2k). Considering Eq. (3.2), the left-hand side has modulus 
less than 1, since k < a < 2k + 1. However, the right-hand side of Eq. (3.2) 
has modulus greater than 2, leading to contradiction. ■ 


Since the roots are simple, we can factor our polynomial in the form 

(0 - r i)(© - r 2 )... (0 - r 2fc )(@ + 2 )h(x) = 2k + 2 _ 


The solution to the differential equation (0 — a)h(x) = x b is [63], 


( x b 


b — a 


+ cx a , if a ^ b 


h(x) = < 


x b log e (a:) + cx a , if a = b. 


Therefore, applying (2k + 1) times the solution, we obtain 


h(x ) = 


2k+ 2 


log e (aO 


2k+l 


(—2 — r i)(—2 — r 2 )... (—2 — r 2 k) x 2 


+ c i xTi ' 

3 = 1 
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where cj are the constants of integration. Note that 


7 ? 2fc+i(©) — (© + 2)<S 2 fc(©), 


thus, the expression in the denominator <S 2 fc(—2) is 


S 2l (-2)=74 +1 (-2) 


2k 


-2 \ ^ 1 . . 


2 k + u u~ 2 ~ j 


— — {2k + 2){H2k+2 — Hk + 1). 


k 'U~ 2 ~j 


Reverting to the previous notation, 


/(*) = 


1 1 °ge( 1 - z ) + 


H 2k+ 2 ~H k+1 (1 -z) 2 ' ^ 


Using the identity [26], 


1 log/ T 2 -) = E(n« +m - H m ) ( n + m 'ir 




n=0 


m 


and the binomial Theorem, the solution of the differential equation can be 
written in terms of series. 


oo 2/c+l 


/(*) = 


H 2 k+2 — Hft + 1 


E((" + 1 )H n - n)z " + X 5 ] c,(-l)’ 


nil I 


n=0 


n=0 j=l 


n 


z . 
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Extracting the coefficients, the expected number of comparisons of ‘median of 
(2k + 1)’ Quicksort is 

2fc+l / 

E(C n { 2 fe+i}) = — -——((n + 1 )H n — n) + (—l) n 9fc [ Cj 

n 2 k +2 — tlk +1 \ 

The real roots of the polynomial 0, 1,. .., (k — 1) do not contribute to the ex¬ 
pected number of comparisons, since n > 2k + 1 and when k is odd, the root 
(3 k + 2) adds a negligible constant contribution. Further, note that the root 
t 2 fc+i = —2, contributes c 2 k+\ (n + 1), with c 2 k+ 1 £ R. Therefore 

E(Cn{2fc+i}) = 77-77 —(( n + 1 )H n — nj + C2fc+i(n + 1) 

tl2k+2 — H k +1 

Uj 

+ 2^(-l)”Mc (c,-(*)) +0(1). 

The asymptotics of the real parts in the last sum is given in [27]. It is proved 
that 



2(-l)"mc( cA Pj ) ) = 0 ( n -WT+')) 


n 


(SHcO^+iy 


and asymptotically the expected number of key comparisons is 


®(C < n{2fc+l}) — 


1 


H 2 k+2 — Hk+1 

+ o(n), 


n\og e (n) + ( c 2k+ 1 + 


H2k+2 ~ Hk+1 


(7-1)) 


n 


since all the other roots have real parts greater than —2. The solution contains 
the case of ordinary Quicksort, where the pivot is randomly selected and of 
‘median of 3’ Quicksort. The coefficient of the leading term in the latter case 
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is 


1 _ 12 
// 4 - Hi ~ 7 k 

thus the expected number of key comparisons is 

E(Cn{3}) = + l)-^n + 0(n). 

van Emden [19] obtained the asymptotic cost of the expected number of com¬ 
parisons, using entropy arguments. The leading term is equal to 

an log 2 (n ), 


where 


1 

= E W) 


2 


1 


xg(x) log 2 (x) dx 


J o 

Here g{x) denotes the probability density of the median of a random sample 
of (2k + 1) elements and E (H) is the expected information yielded from a 
comparison. The nice and simple asymptotic form for the mean number of 
comparisons is 


E(C'n{2fc+1 } ) ~ lQ ge( 2 ) - n l 0 g 2 ( n ). 

Xi'2k+2 ~ n-k+l 

Note that this result, yields the expected number of comparisons of standard 
Quicksort, for k = 0. In this case, 


a = 2 log e 2, 
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thus 

E (C n ) ~ 1.386n log 2 (n). 

The notion of entropy is of great importance to the analysis presented in this 
thesis. Entropy is a measure of uncertainty regarding the events of a random 
variable. In other words, a higher uncertainty about the outcome of a random 
variable pertains to increased entropy. In the trivial case, where the proba¬ 
bility of occurrence of an event is 1, the entropy is equal to 0, as there is no 
uncertainty. A formal definition follows [65], 

Definition 3.1.2. The Shannon’s entropy H of a discrete random variable X tak¬ 
ing the values x±,x 2 ,... ,x n with probabilities p(x\),p(x 2 ),... ,p(x n ) is defined 

by, 

n 

H{X) = - ^pix^Xag^pixi)). 

i= 1 

The base b of the logarithm will be normally equal to two; in this case we 
measure bits of entropy. We should mention that the notation H(X) does 
not merely denote a function of X; entropy is a function of the probability 
distribution. Generally, in sorting algorithms that utilise comparisons for this 
task, entropy quantifies the amount of information gained from the sorting. 
Consider an unsorted array, with all the n! permutations equally likely. A 
comparison gives 1 bit of information, thus at least log 2 (n!) = fl(nlog 2 (n)) 
comparisons are needed for a complete sort - see in [5], This quantity is 
called information-theoretic lower bound. The range of entropy is given in 
the following Lemma. 
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Lemma 3.1.3. 


0 < H{X) < log b (n). 

Proof. Since log b p(xi) < 0, the left inequality follows immediately. Noting 
that the logarithm is concave function and applying Jensen’s inequality [38], 
which for a random variable X and a concave function /, states that 

/(EPO) >E(/(X)), 


we have 


n 

^2p(Xi) logfe 

i=1 



The next four definitions can be found in [16], [53] and [65]. 

Definition 3.1.4. The joint entropy H(X,Y) of two discrete random variables 
X and Y is defined as 

H(X, Y) = — EE p(x,y) log b p(x,y), 

x€X y£Y 

where p(x, y) is the probability that X takes the value x and Y the value y. 
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Definition 3.1.5. The conditional entropy H(X\Y) of two discrete random vari¬ 
ables X and Y is defined as 

H(X\Y) = -'£'£ p{x,y ) logfa P< Y y • 

xex y eY 

Definition 3.1.6. The information content of a random variable X, with proba¬ 
bility distribution P(X) is 


I(X) = - log 6 P(X). 

From definition 3.1.6, one can easily deduce, that 

H(X) — E(/(A r )). 

In other words, entropy is the expected value of the information. We proceed 
to the definition of mutual information. 

Definition 3.1.7. The mutual information of two discrete random variables X 
and Y is defined as: 

I(X A Y) = 5] y>(z, y) log, -fhf- = H(X) - H(X\Y) 

ttxjtv pix)p(y) 

and quantifies the amount of information provided about X by Y. 

These definitions will be used in a later part of the thesis. 

We have to point out that up until now, our analysis did not take into account 
the added overhead of finding the median at each stage. In the simple case 
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of three elements, the overhead is not significant, but for larger samples this 
might have adversary effects to the efficiency of Quicksort. 

For the selection of the median, we use Hoare’s Find algorithm [29] or Quick- 
select. This simple and intuitive algorithm searches for an element of a given 
rank m in an array of n keys. As in Quicksort, one partitions the array around 
a randomly chosen pivot, which at the end of the partition process is moved to 
its final position, j. If m = j, the pivot is the sought element and the search is 
completed. Otherwise, if m < j, Quickselect is recursively invoked to the left 
subarray of j — 1 keys. Conversely, if m > j, we search in the right subarray 
for the element of rank (m — j). 

Quickselect is ideal in situations where we want to identify order statistics, 
without the need to do a complete sort. The average number of comparisons 
E( C n - m ) required for the retrieval of the m-th order statistic in an array of n 
keys, is given by [46] 

E(C n;m ) = 2 (n + 3 + (n + 1 )H n — (m + 2 )H m — ( n — m + 3)iF n+ i_ m ). 

For the sample of (2/c+l) keys, the rank of the median is k+1. Therefore, 

®(C , 2fc+i ; fc+i) = 2 (2 k + 4 + (2k + 2)H 2 k+i — 2 (k + 3)Hk+i) . 

The cost for finding the median of a random sample of 3 keys is 8/3 com¬ 
parisons at each stage. However, as the sample of keys increases, in order to 
obtain a more accurate estimate of the median, the added overhead imposes a 
bottleneck to the efficiency of the algorithm. 
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The computation of other measures of this variant, such as the expected num¬ 
ber of exchanges and passes can be performed by solving analogous recur¬ 
rences, as the one for the number of comparisons. The average number of 
passes is recursively given by 


E(-P n {2fc+1}) — 1 + 


n 

2k+ 1 


£ 

3 = 1 


J ~ 1 
k 


n 


E(Pj-l{2fc+l})^ 


where the “toll function” is now one recursive call to the algorithm, after the 
chosen pivot is the median of (2k + 1) keys, which yields two subarrays. This 
recurrence can be turned to a differential equation with solution 


! j — j 

lE(T > n,{2fc+l}) ITT TT \ 

£\n-2k+2 — n-k+1) 

noting that for k = 0, the average number of passes is n. In the scheme, where 
arrays containing m or fewer keys are sorted by insertion sort, the average 
costs are reduced. We refer to [27] where this variant is analysed. 


3.2 Remedian Quicksort 

In the previous section, we analysed the modification of Quicksort, where the 
pivot is selected as the median of a sample of (2k + 1) elements. This variant 
offers better protection against the occurrence of trivial partitions. However, 
there are some cases, where the running time of this partitioning scheme can 
go quadratic. Consider the application of ‘median of 3’ Quicksort in an array of 
n numbers, where the keys at positions 1, n and are selected as elements 
of the sample. In case that two keys of this sample happen to be the smallest 
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(or greatest) elements of the array, the chosen pivot will be 2 (or n — 1), lead¬ 
ing to trivial partitioning, making Quicksort everything else, except quick! In 
[63], a permutation of the array (1, — , 15} is given, which leads Quicksort to 
worst-case performance and in [20], an algorithm is presented which forms 
the worst-case permutation. 

In order to remedy this, a bigger sample of (2 k + 1)^ keys is randomly se¬ 
lected, its remedian is found and used as partitioning element of the array to 
be sorted. The remedian of the sample is recursively defined to be the median 
of (2 k + 1) remedians of (2 k + l)' 3 ” 1 elements, where the remedian of (2 k + 1) 
elements is the median, and is shown to be a robust estimator of the median 
[60], [72], 

A particular case of the remedian Quicksort widely used in sorting applications 
is Tukey’s ‘ninther’, where the selected pivot is the median of three medians 
of three samples, each containing three elements [14], [72], In practical im¬ 
plementations, this variant exhibits faster running time [7] with little added 
overhead. Specifically, the computation of the remedian of 9 elements takes 
on average 4 x | comparisons at each call - four times more than finding the 
median of 3 randomly chosen keys. 

Let C n {( 2 fc+i)/n denote the number of comparisons required for the complete 
sorting of an array of n distinct keys, where the chosen pivot at each call is the 
remedian of a random sample of (2 k + 1)^ elements. The recurrence relation 
is much more complicated than the previous ones and the probability p? that 
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the remedian of (2k + 1)P elements is the (j + l)-th element is 


Pi = 


(2 k+l)\ Y P( a i> 


i «2fc+l J 


ai+...+a 2k+1 =j 


j W (2k+l)@ —j—1 > 

<xi,---,a 2 k+i) V(2fc+1)/ 3- ! —1—ai,...,(2fc+l)^ — !—«2fc+i'' 
(2fc+l)/5 \ 

(2fc+l)/ 3 - 1 ,...,(2fc+l)/ 3 - 1 / 


where 


. . . , «2fc+l/ Oil! • • ■ ®2fc+l! 

is the multinomial coefficient and p(a i,..., 02 ^+ 1 ) is defined by 


P(«l> • • • > «2fc+l) = PS 1} (P0 /3 13 + • • • + P£-1 ) ' • • • ' (Pil+J + • • • + P(2fc+!)/3-l)> 

with = 1. See as well [14], where the ‘splitting’ probabilities of 3 d reme¬ 
dian are presented. 

Bentley’s and Mcllroy’s experiments on the ‘remedian of 3 2 ’ Quicksort [7] 
showed that the average number of key comparisons is 1.094n log 2 (n) — 0.74n, 
very close to the information-theoretic lower bound of nlog 2 (n) — 1.44n. This 
Quicksort utilises the ‘ninther’ partitioning for large arrays, then the ‘median of 
3’ is used and as the algorithm proceeds, the partitioning strategy changes to 
the standard uniform pivot selection. The paper written by Durand [18] con¬ 
firmed these experimental results, where the average number of comparisons 
is being given by 


E(C n{3 2 } ) = 1.5697nlog e (n) - 1.0363ra + 1.5697log e (n) - 7.3484 + O 

From a theoretical point of view, this variant yields savings on the expected 
time needed for the sorting, with little additional cost of computing the re¬ 
median. However, the recurrences are quite involved, as the remedian has 
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an inherent recursive definition. A different approach would be to randomly 
choose a larger sample and use its elements as pivots through complete sort¬ 
ing. An obvious advantage of this method, is that the cost of computing the 
median or remedian of a sample at each call of the algorithm is avoided and 
instead all the pivots for the subsequent calls belong in one sample. 


3.3 Samplesort 

Having examined the strategy of selecting the median of a sample as pivot and 
the more complicated ‘remedian’ variant, we proceed to the analysis of Sam¬ 
plesort algorithm invented by Frazer and McKellar [24]. Instead to randomly 
select a sample of keys at each stage, computing the median and using it as 
pivot, a larger sample of 2 k — 1 keys is selected and extracted out of the array. 
It is sorted and its keys are being used as partitioning elements for the sorting 
of the array. 

First the median of the sample is used as pivot, then the lower quartile to 
the lower subarray and the upper one to the subarray of the elements that 
are greater than the median. When the sample is exhausted, the resulting 
subarrays can be recursively sorted by the same procedure or by standard 
Quicksort. 

For convenience, let 2t + 1 = 2 k — 1 and let us denote the total number of 
comparisons of Samplesort applied to n keys by Cn 2 1 b Then, G'i. 2 1 ' is equal 

to the number of comparisons to sort the sample of 2 k — 1 elements, plus the 
number of comparisons to insert its elements to the remainder of the array, 
plus the number of comparisons required to sort the resulting 2 k subarrays, 
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using ordinary Quicksort [24]. For the sorting of the sample, we use Quicksort 
and the average number of comparisons is 

2(2 1 + 2)H 2 t+i — 4(2 1 + 1). (3.6) 


Assume that the sorted sample is x\ < x 2 < ... < x t+ i < ... < x 2t+ 1 - First, 
the median x t +i is inserted to its final position in the array of n — 2 t keys, by 
pairwise comparisons of the n — 2t — 1 keys to x t +i- The cost of partitioning is 
n — 2t — 1 comparisons. Then, the first quartile is inserted to the subarray of 
the elements less than the median x t +\ and the third quartile to the subarray 
of the elements greater than x t +\. The cost of partitioning these two subarrays 
is n — 2t — 1 comparisons, since the sum of their lengths is n — 2t — 1 keys. The 
process is continued until all the elements of the sample are used as pivots and 
this will take 2t +1 partitioning stages. Thus, an approximation to the average 
number of comparisons for the insertion of the sample is [24], [49], 

(n — 2t — 1) log 2 (2t + 1). (3.7) 

After all elements have inserted, there are 2(t + 1) subarrays to be sorted by 
Quicksort. The expected number of comparisons is [24] 


2 (n + l)(H n+ i — H 2 (t+i)) — 2(n — 2 1 — 1). 


(3.8) 
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Putting together Eq. (3.6), (3.7) and (3.8), we have that the expected number 
of comparisons of Samplesort is 


2 (n + l)(H n+ i - H 2 (t+ 1 )) - 2 (n - 2t - 1) + (n - 2t - 1) log 2 (2f + 1) 

+ 2(2f + 2)H'2t+i — 4(2 1 + 1). (3.9) 

For large values of n, the expected number of comparisons is given by the 
following Corollary. 

Corollary 3.3.1 (Frazer and McKellar [24]). The asymptotic expected number 
of comparisons taken by Samplesort for the sorting of an array of n keys, using a 
randomly chosen sample of l keys, is 

E(C'W) = 1.386nlog 2 (n) - 0.386 (n - l) log 2 (l) - 2 n - 0.846/. 

It is worthwhile to note that the programming of Samplesort is simple and 
straightforward, thus making it a suitable candidate to sorting applications. It 
is proven in [24], that the procedure is asymptotically optimal, i.e. as n —» oo, 
the expected number of comparisons approaches the information-theoretic 
bound. The process of randomly drawing the sample out of the array to be 
samplesorted should be carefully selected, to the effect that the elements of the 
sample will produce ‘balanced’ partitions. The main feature of Samplesort is 
that partitioning preserves the order of the sample, thus its keys are exchanged 
with the ones that they have to, so as to the exhaustion of the sample, its 
elements to be spread far apart. 



85 


In the direction of choosing a more ‘centered’ sample, one can proceed by ran¬ 
domly selecting three samples - each containing three keys - and computing 
their medians, at an extra cost of 8 comparisons. For the sorting of larger ar¬ 
rays, the number of samples will obviously be greater. The medians will be 
used as elements of the sample and the remaining elements can be randomly 
chosen from the array. It should be noted, that this is an initial idea, lacking 
the mathematical analysis. 

Albacea [2] derived a modification of Samplesort. This variant starts with 1 
key, that is used as pivot for the partition of 2 keys, so to have a sorted array 
of 3 keys. These keys are used as pivots for the partitioning of 4 keys, so as to 
obtain a sorted array of 7 keys and so on, until the whole array is sorted. It is 
proven [2] that the estimated expected number of key comparisons is 

n|~log 2 (n + 1)] — 2^ log2 ^ n+1 ^ —n + |~log 2 (n + 1)] + 1. 

The derivation of the expected number of comparisons remains an open prob¬ 
lem. 



Chapter 4 


Sorting by multiple pivots 


In this Chapter, different partitioning routines are analysed. These schemes 
utilise many pivots for the partitioning of the array and naturally arise as a 
generalisation of the algorithm. The pivots are chosen uniformly at random 
and the array is partitioned into more than two subarrays. There is an addi¬ 
tional overhead of comparing the pivots before the partitioning, which adds 
very small contributions to the running time. The aim of this modification 
is twofold: first to study if the possibility of worst-case scenario of the algo¬ 
rithm can be reduced further and secondly, to provide a theoretical basis for 
the analysis of the generalisation of the algorithm. 

We show that the average case analyses of these variants can be fully described 
by a general recurrence model, which is transformed to a differential equa¬ 
tion, whose solution provides the expected cost of these variants. Further, we 
demonstrate that the integration constants involved in the solution, can be 
efficiently computed using Vandermonde matrices. 
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4.1 Quicksorting on two pivots 

Along the following lines, we present a variant of Quicksort, where 2 pivots are 
used for the partitioning of the array. Let a random permutation of the keys 
{1,2,..., n) to be sorted, with all the n! permutations equally likely and let 
their locations in the array be numbered from left to right by {1,2,..., n}. The 
keys at locations 1 and n are chosen as pivots and since all the n\ permutations 
are equally likely to be the input, then all the Q) pairs are equiprobable to be 
selected as pivots. At the beginning, the pivots are compared each other and 
are swapped, if they are not in order. If elements i < j are selected as pivots, 
the array is partitioned into three subarrays: one with (i — 1) keys smaller than 
i, a subarray of (j — i — 1) keys between two pivots and the part of (n — j) 
elements greater than j. 

The algorithm then is recursively applied to each of these subarrays. The 
number of comparisons during the first stage is 

Ai .,2 = 1 + ((i — 1) + 2 (j — i — 1) + 2 (n — j)) 

= 2n — i — 2, 

for i = 1,..., n — 1, and j — i + 1,..., n. Note that in the specific partitioning 
scheme, each element is compared once to i and elements greater than i are 
compared to j as well. The average number of comparisons for the partitioning 
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of n distinct keys is 


E(A n ,2) — / \ ^2 ^2 ( 2n * 2 
( ] *=1 j=i+i V 


2 / 5 , o 7 

-77 7+ ~ 2n + r n 

n(n — 1) \o o 


5 n — 7 


Letting C n 2 denote the number of comparisons of dual pivot Quicksort applied 
to an array of n items, the recurrence for the expected number of comparisons 


, 5n - 7 2 

^(C n , 2 ) — -1-7-7T 

3 n[n — 1) 


5] X e(c,.,, 2 )+ e E e (Cj-*-w) +EE e < c »-ja 

2=1 j = i -\-1 2=1 7 = 2+1 2=1 7 = 2+1 


Note that the three double sums above are equal. Therefore, the recurrence 


becomes 


^77 — 7 ft 

E(C„, 2 ) = —— + --- E(« - mCi- 

3 nyn — 1) 


Letting a n = E((7 n 2 ), we have 


5 n — 7 6 \ ^ 

a n = —--h —r —~ > ^(n - n > 2. 


3 n(n — 1) ' 

v 7 2=1 



89 


It holds that CLq — Ci\ — 0. Multiplying both sides by ( ), we obtain 


n 


n 


0>n. 


n—1 


n\ t 5n — 7 6 

7j — + — ( -7T > (n - l 

3 ii(n — 1) z —' 


n(n — 1)(5 n — 7) 


n—1 


6 


3^(n - i)di- 1. 




This recurrence will be solved by the difference method. We have 

A F(n) := F{n + 1) — F(n) and for higher orders 

A k F(n) := A k ~ l F{n + 1) - A fc “ 1 F(n). 

Applying the difference operator 


. n \ f n + 1\ /n\ 5n 2 — 3 n 

^ I ~ I &n — I 9 / V 2 / ^ 71 = 2 ^ ^ / y ^ 

/ ^ / i =o 

A 2 f jM a n = A ^ a n+ i - A Q a n = 5n + 1 + 3 a n . 


By definition, 


A 2 f 2 Vn = A r 7 2 1 J a n+ i — A f ^ ) a n 


™ + 2\ Ar+1\ /V 

2 ) a n+2 21 ^ ) a n +1 + I 2 ) a r 


and the recurrence becomes 


(n + l)(n + 2)a n+ 2 — 2 n(n + l)a n+1 + n(n — 1 )a n = 2(5n + 1 + 3a n ) 

==>■ (n + 1) ((n + 2)a n+2 - (n - 2)a n+ i) - (n + 2) ((n + l)a n+1 - (n - 3)a n ) 
= 2(5n + 1). 
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Dividing by (n + 1 )(n + 2), we obtain the telescoping recurrence 

(n + 2 )a n+2 - (n - 2)a n+1 _ (n + l)a n +i - (n - 3 )a n 2(5 n + 1) 

n + 2 n + 1 (n + 1) (n + 2) ’ 


which yields 


(n + 2)a„, +2 - (w - 2)a n+ i 
n + 2 




i=o 


5 j + 1 

(j + l)(i + 2) 


18 

n + 2 


+ 10H n+ i - 18. 


The recurrence is equivalent to 


na n — (■n — 4)a n _i = 18 + l()n// n _i — 18n. 


Multiplying by — 
scoping one [63], 


l)(n — 2 )(n 
24 


3) 


this recurrence is transformed to a tele- 



( 4 ) a «-' + 


18(n 


l)(n — 2 )(n 
24 


3) 




Unwinding, we have 


n 


&n, 


is ± (j - i)(j ~ 2)(j - 3> + io ± a ) 


1=1 


1=1 


l=i 


J 

4 


24 


(4.1) 
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The second sum of Eq. (4.1) is 



The expected number of comparisons, when two pivots are chosen is 

a n = 2 (n + 1 )H n — An. 

This is exactly the same as the expected number of comparisons for ordinary 
Quicksort. 

Next, the expected number of key exchanges will be computed. The exchanges 
during partitioning are performed as follows. Set two pointers l 4— 2, u •<— n —1 
and store temporarily the pivots in another array of size two, so the cells at 
locations 1 and n are empty, leaving two “holes”. 

After the pivots are sorted by one comparison, the key at position 2 is compared 
to the first pivot (i.e. the smaller of the two pivots); if it is less than the pivot, 
it is put into the left hole, which now is moved one position to the right and 
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l is increased by one. If it (i.e. the key at position 2) is greater than the first 
pivot, it is compared to the second pivot (i.e. the greater of the two pivots). If 
it is less than the second pivot, l is increased by one, otherwise / stops. 

Now, the u pointer starts its downward scan. If an examined key is greater 
than both pivots, it is put into the right hole, which is moved one position to 
the left and u is decreased by one. If a key is less than the second pivot and 
greater than the first, then u is decreased by one. In case that a key is less 
than the first pivot, then u stops its scan and the key that is greater to the 
second pivot, where / has stopped is put to the right hole, which is moved one 
position to the left and the key where u has stopped is put to the left hole, 
which is moved one position to the right. Then, l is increased by one, u is 
decreased by one and l resumes its scan. 

When pointers are crossed, the first pivot is put to the left hole, the second 
pivot is put to the right hole and partition is completed, since keys less than 
the first pivot are on its left, keys between two pivots on the middle and keys 
greater than the second pivot are on its right subarray. Note that the auxiliary 
space required for the storing of pivots is 0(1), since at the end of the partition 
routine, the pivots are moved back to the array and two other new pivots can 
be stored, as the algorithm operates on a given subarray. We refer to [63] for 
further details of this scheme. 


The average number of swaps during the first stage is 



n— 1 n 


i = 1 j = i -\-1 


1 n— 1 
0 - 1 



'n —1 n —1 

(n — i)i — ^^(n — i) 
. 2=1 1=1 
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since (i — 1) keys are less than pivot i. Thus, the average contribution is 

2 (n 3 — n n(n — 1) ^ n — 2 

n(n - 1) V 6 2 ) ~ 3 

For the (n — j) keys greater than j the average value is the same, because the 
sums are equal. Adding the two final “exchanges” to get the pivots in place, the 
average number of exchanges during the partitioning routine is 
Letting S n , 2 denote the number of exchanges of dual pivot Quicksort, the re¬ 
currence for the mean number of exchanges in course of the algorithm is 




E(S' n , 2 ) = 


2 (n + 1) 2 


3 ’ n(n — 1) 

n —1 n n —1 n 

x i E E +E E E (v<-1.2) 


, 2—1 j=i -\-1 

2 (n + 1) 6 


2—1 


n— 1 n 

E E e < s ”-j.^ 

2=1 j =2+1 


n—1 




71(71 — 1 ) ^ 

v 7 i=i 


Multiplying both sides by 




n(n 


l)(n + 1) 
3 


n— 1 

+ 3^(n-i)E(A,_i, 2 ). 
2 — 1 


This recurrence is solved in [63]: here we present a solution using generating 
functions. Letting b n = E(5' nj2 ) and g{z) = b n z n be the generating func¬ 
tion of the average number of exchanges, the recurrence is transformed to the 
following differential equation: 


d 2 g(z) 
2 d z 2 


z 2 d 3 
3 dz 3 



v n+l 


00 n 

aYsYXn-i)^. 

n =1 i=l 
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The double sum is equal to 


Y YX n ~ i)bi-\z n = b 0 z 2 + (2 b 0 + b^z 3 + (36 0 + 2&i + b 2 )z A + ... 

n =1 i= 1 

= z~(bo + + •••) + 2 z 3 (bo + 61Z + b 2 z 2 + ...) + ... 

= (V + 2z 3 + 3z 4 + .. .)g(z) 



and our differential equation becomes 


£ 2 d 2 g(z) 
2 dz 2 


2 z 2 

(1 - *) 4 


+ 3 g(z) 


2 


Changing variables v — 1 — 2 , we have / |A : ( v) — i -1 i' '/ ' 1 (1 — v). Thus, 


(1 — v) 2 d 2 f(v) 2(1 — a) 2 , „„ . fl-v 
— -+ 3 /<” 


The differential equation can be simplified by multiplying both sides by 


2 


v 2 d 2 /(v) 
2 dv 2 


4 + 3 /<t)- 

v 


(4.2) 


An elementary approach to solving this differential equation, is to assume that 
the solution is of the form x m [9]. Substituting the “trial solution” to Eq. (4.2), 
the characteristic or indicial polynomial is 


V 2 (m) = m(m — 1) — 6, 
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with roots rri\ = 3 and m 2 = —2. Thus, the solution to the corresponding 
homogeneous equation is c\v 3 + c 2 v~ 2 , with cj,c 2 6 R. A particular solution of 
Eq. (4.2), which can be found e.g. using the method in [61], is 


log e (v)v 2 


By the initial conditions /(1) = —/'(1) = 0, the solution is 


f(v) 



20 log e (u) + 4 
25t4 


In the next section, we will examine the generalised version of this differential 
equation. Reverting to variable z and discarding terms for n < 3, we see that, 
expanding out the fraction term as a series, 


9{z) = 


OO 

E 

n=0 


— ( (n + 1 )H n 


n 


25 (n+1 ) 


Finally, the mean number of swaps of dual pivot Quicksort is 



(n + 1 )H n 


24 n + 4 
25 


which is nearly 2.4 times greater than the expected number of exchanges of 
standard Quicksort. 

The recurrence for the number of partitioning stages P n)2 is much simpler; 


P'ri.2 — 1 + Pi-1,2 + Pj-i- 1,2 + Pn-j. 


3 , 2 ' 
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By the same reasoning, as in the derivation of the expected number of ex¬ 
changes, the solution is 


E(Pri,2) — + 1) - 


4.1.1 The variance of the number of key comparisons 


It is desirable to compute the variance of the number of key comparisons of 
dual pivot Quicksort, as this measure provides a grip of the deviation of the 
random number of comparisons from its expected value. By the recursive 
relation, we have 


n— 1 n 

P(Cn,2 = t) = yyy y P(^n,2 + Cj-1,2 + Cj-i-1,2 + Cn-j',2 = f), 
2=1 j=i -\-1 


noting that the resulting subarrays are independently sorted, the above is 


n —1 n 


n 


i =1 j=i+l l,m 


EEE P(CU ,2 — /)P(C'j-j-i )2 — m)P(C ' n _ Ji 2 — t — m — l — 2 n + i + 2) 


OO 

Letting f n (z) = ^P(C n , 2 = t)z f be the ordinary probability generating func- 

t =o 

tion for the number of comparisons needed to sort n keys, we obtain 


fn(z) 


n— 1 n 

2=1 ji=2+l 


n 

2 


(4.3) 
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It holds that /„( 1) = 1 and /'(l) = 2(n+l)iT n —4n. The second order derivative 
of Eq. (4.3) evaluated at z = 1 is recursively given by 



The fourth and fifth sum turn out to be equal and by simple manipulation of 
indices, the sums involving products of expected values are equal. The double 
sum of the product of the mean number of comparisons can be simplified as 
follows, using Corollary 2.3.5: 
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Further, using 


n — i + 1 


+ (n - i). 


Zi( n ~i +1 )= Ef(i-l) + l)f"?’ +1 ) 


= J2 (i - !) (" 2 ‘J + £ (” 2 j H> -' H 

n— 1 n—1 

- l)(n - + ^(n - 


The four sums can be evaluated using Corollary 3 in [70], 

After some computations in Maple, that can be found in Appendix A, the 
recurrence is 


/n(l) — 2(n + l)(n + 2) (if 2 — H — H n n 2 + —n + 6^ + 

731 13 6 ^ 

+ + ¥ + ^ (n “ ?)/t - l(1) ‘ 


Subtracting Q /"(!) from ^ ^ /"+i(l), we have 


A (2) = 4n ( n + + } ) “ ^p( 84 ™ 2 + !98n + 42) 

n 

+ 3^/'^ + ¥79n 2 + 231n + 14), 


using the identity [63] 


tt ‘2 _ rr(2) _ tt 2 _ rr(2) ,_ 

■*■*71+1 -^*n+l -*■*n ^ ^ _|_ ]_ * 
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Also, it holds that 

A 2 Of) m = 12(n + l)(n + 2 )(ff„ 2 - ) _ ff „(20n 2 + 32n - 12) 

+ 17n 2 + 37n + 3/"(l). 

The left-hand side of the previous equation is the same as 

and the recurrence becomes 

(n + 1 )(n + 2)/" +2 (l) - 2n(n + l)/" + i(l) + n(n - 1)/"(1) 

= 2 f 12(n + l)(n + 2)(H 2 - //f 2) ) - 77 n (20n 2 + 32n - 12) + 17n 2 + 37 n + 3/"(l) j 

Dividing by (n + l)(n + 2), we obtain the telescoping recurrence 

(n + 2)/; + . ; (l) - (n - 2)/" +1 (l) 
n + 2 

, ( n + l)/n + l(l) ~ ( n ~ 3)/n(l) 

n + 1 

H n (20n 2 + 32n - 12) 17n 2 + 37n \ 

(n + 1) (n + 2) + (n + 1) (n + 2) y ’ 

with solution 

(n + 2)/" +2 (l) - (n - 2)/" +1 (l) = (24n 2 + lOOn + 104)(fl^, - A,2i) 

- 77„+i(88n 2 + 292n + 224) + 122n 2 + 346ra + 224, 
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which is equivalent to 

<(1) - (ft - 4)/"_j(l) = (24n 2 + 4 n)(Hl_, - H^\) 

- H n _i(88n 2 - 6On - 8) + 122n 2 - 142n + 20. 

. . , c . L1 -j , (n — l)(n — 2)(n — 3) . 

Again as before, multiplying both sides by -——-, the recur¬ 

rence telescopes with solution 

/"( 1) = 4(n + 1 ) 2 (tf 2 +1 - H^) - 4 H n+1 (n + l)(4n + 3) + 23n 2 + 33n + 12. 

Using the well known fact that 

Var(Ga) = /:(!)+ /(,(!)-(/(.(l)) 2 , 


the variance of the number of key comparisons of dual pivot Quicksort is 

7n 2 - 4(n + 1) 2 ^ 2) - 2(n + 1 )H n + 13n. (4.4) 

Note that the variance of dual pivot Quicksort is identical with the variance of 
ordinary Quicksort. In the next subsection, we provide the theoretical expla¬ 
nation of this fact. 

4.1.2 Distribution of the number of key comparisons 

Our results have shown that dual pivot Quicksort has the same expected num¬ 
ber of comparisons, and the same variance, as in the case of ‘one-pivot’ Quick¬ 
sort. Thus, it is natural to ask if the two random variables have the same 
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distribution. We now show this, after an argument sketched by Prof. Colin 
McDiarmid [52], 

Suppose that an array of n distinct keys x\, x 2 ,..., x n is to be sorted by Quick¬ 
sort and let, as usual, C n be the random number of comparisons required for 
the sorting. Obviously, C\ — 0, C 2 — 1 and for n > 3 we pick uniformly at 
random an ordered pair of distinct indices (/, J) in [n] = {1,2,..., n} and we 
use xi as the first pivot. Given that I — i, the pivot x, partitions the array of n 
keys to the subarray of (i — 1) keys less than x, and to the subarray of (n — i) 
keys greater than x t by (n — 1) comparisons. 

Given that I — i, if xj < x i} then xj is a uniformly at random chosen pivot 
from the subarray of (i — 1) elements less than x t . In this case, for I = i and 
J = j, the subarray of (i — 1) keys is partitioned to the subarray of (j — 1) 
keys less than xj and to the subarray of (i — j — 1) keys greater than Xj by 
( i — 2) comparisons. Note that ( i — 2) keys are compared to both pivots in two 
partitioning stages. Therefore, the following recurrence holds: 


n— 1 n 


P(Q„, — t)— / x ^ ^ P((n - 1) + (i ~ 2) + i + C^li — t ) 


k3) 


j =i *=l+i 


n— 1 n 


n / n _ x') P (( n + 1 ~ 3 ) + C j-1 + C i-3 -1 + C n-i ~ *) , 

K ’ 1 = 1 *= 1+1 


where Cn \ Cn' 1 and CrP are independent copies of C n - that is, are ran¬ 
dom variables with the same distribution as C n , independent of it and each 
other. 


If xj > Xi, then xj is a uniformly at random selected pivot from the subarray 
of (n — i ) keys greater than x^. Given that I = i and J = j, the subarray of 
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(n — i) keys is partitioned to the subarray of (j — i — 1) keys less than x 3 and 
to the subarray of (n — j) keys greater than x j by (■n — i — 1) comparisons. The 
recurrence relation is 


n— 1 n 

P(C n = t) = ^(( n — 1) + ( n — — 1) + C-J 1 + Cj]j _i + C„lj = t ) 

2—1 ^‘=2+1 


n—1 n 


E E p « 2 ™ - - 2)+<£’.+cfE,+= t ), 

' 2 i=l j=i+l 


where as in the previous recurrence, Cn\ C'n ' 1 and Cn'* are independent copies 
of C n . Observe that 


n— 1 n n— 1 n—1 n 

= _ N 2n = ( 2n _ * - 2 )’ 

j=l i=j+l j =1 i=l j=i+l 

thus for any two pivots selected uniformly at random, the recurrences are the 
same. 

Recall that for the random number of comparisons C n> 2 of dual pivot Quicksort, 
it holds that C\ 2 = 0, C 2)2 = 1 and for n > 3 we choose uniformly at random 
an ordered pair of distinct indices (/, J) in [n] = {1, 2,..., n}. The pivots 
xi and xj are sorted by one comparison and we assume that its outcome is 
xi < xj. Given that I = i and J = j, the array is partitioned to the subarray 
of (i — 1) keys less than x i} the subarray of (j — i — 1) keys between two pivots 
and the subarray of (■n — j ) keys greater than Xj . Since keys greater than Xi are 
compared with the other pivot as well, the recurrence for the random number 
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of comparisons is 


Ti—l n 


P(C„, 2 = t ) = £ £ P(1 + (*-!) + 2 (j - i - 1) + 2 (n - j) 


i=l j=i -\-1 


+ cfh + <?£-!, 2 + OS*. = *) 


■( 2 ) 
7i—l n 


(3) _ 


- n(n _ ^ S P (( 2n - * - 2 ) + ^-1,2 + ^i-i-1,2 + C'i-j,2 ~ *) , 

' ’ i =1 j=i+l 

where C'^, C^ \ and C^ \ are independent copies of C n 2 . Note that when x r > 
xj, the recurrence is the same. Thus, since dual pivot Quicksort and ordinary 
Quicksort satisfy the same recurrence and have the same initial conditions for 
n = 1,2, we deduce that the random variables C' n 2 and C n are identically 
distributed. 


4.2 Multikey partitioning 

A natural extension of having two pivots would be to have some other number 
k of pivots. Here, we study the idea of randomly picking k pivots i 1; i 2 ,..., ik 
and partitioning the array simultaneously according to these. 

Again, let a random permutation of the array {1,2,..., n} be given to be sorted 
using this variant, with all the n! permutations equally likely to be the input. 
The k rightmost keys are chosen as pivots, are compared to each other and 
exchanged, if they are out of order. The sorting of the pivots can be effi¬ 
ciently implemented by insertion sort. Since all n\ permutations of the keys 
are equally likely to be the input, this amounts to the fact that any k -subset of 
keys has equal probability to be selected. 



104 


The remaining (n — k) keys are compared to the pivots and the array is parti¬ 
tioned to (k + 1) subarrays. The partitioning can be performed as follows. We 
compare the leftmost key to a randomly chosen pivot; if it is smaller than this 
pivot, it is compared with another smaller pivot (if one exists). Otherwise it is 
compared with a larger pivot (to the right) and after a series of comparisons, 
is inserted to its place between any two pivots, or to the left of the smallest 
pivot or to the right of the greatest pivot. We continue in the same fashion, 
until all keys are examined. 

In [28], each of the (n — k ) keys is compared to the pivots by binary search, 
so a key is compared first to the median of the sorted array of the pivots. If 
it is less, is compared with the first quartile, otherwise is compared with the 
third quartile and after a series of comparisons is inserted to its position. In 
worst case, it takes 0(log 2 (/c)) comparisons for the insertion of a key. Then, 
multipivot Quicksort is recursively applied to each of the resulting segments 
that contains at least (k + 1) keys and arrays with less than (k + 1) keys are 
sorted by insertion sort in 0(n) time. 

This is equivalent to [k + l)-ary search trees, which is a generalisation of binary 
trees. Indeed, if n > k + 1, the k pivots are stored in the root node of the tree 
in increasing order and the remaining (n — k) keys are placed in the resulting 
(k + 1) subtrees of the root. In case that n = 0, the tree is empty and if 
n < k, the tree has a single node, which stores the keys in order. Under the 
assumption of uniformity, this is a (k + l)-random tree of n nodes - see the 
article of Chern et al. [14] and Mahmoud’s book [49] for the correspondence 
between trees and variants of Quicksort. 
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Let f(n, k) denote the expected cost of the algorithm applied to an array of n 
keys. We deliberately allow some flexibility in the form of cost; a typical ex¬ 
ample might be the number of comparisons. The expected cost of this variant 
is recursively given by 

f(n, k) = T(n, k ) 

+ T^Y (/(*! “ k ) + /(*2 -h~l,k) + ... + f(n - 4 , k)^j , 

—<*fc 

where i[ < i' 2 < ... < i' k are the pivots in increasing order, T(n, k) = a{k)n + 

6(/c) is the average value of a “toll function” r(n, k) during the first recursive 
call and f(i\ — 1, k) denotes the average cost for sorting the subarray of (i\ — 1) 
elements less than i\ by multipivot Quicksort on k pivots. 

Though this looks a complex k -index summation, the recursion can be sim¬ 
plified, by noting that the pivots are randomly selected and the sums are 
equal, 
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n 


Multiplying both sides by , the recurrence relation becomes 

k 


n A f(n, k ) = V(n, k) + (k + 1) y ^ ) f{%\ - 1, k). 


*l=i 


For notational convenience, let f(n, k ) = a n and consider the generating func- 

OO 

tion h(x) = ^ a n x n ; 


n=0 


E 

?i=0 


n 


a n x n = 


n=0 


oo / n 


n =0 \i'=l 


£(")r(n,*)z” + (fc + l) EIECE 1 


, . a,;' _ 1 X 

k- 1 1 1 1 


The recurrence is transformed to a k -th order differential equation 


h^ k \x)x k 

k\ 


n =0 
oo 


Y / ( h )T(n,k)x n + h{x)(k + l)J2( k _ 1 


n =0 


n — 1 


x' 


1 — x 


Y, ( a{k)n + b(k))x n + (k + l)h(x) ^ 

x k (a(k)(x + k) + b( k)( 1-x)) | | / x 


(1 -x) 


1 — x 


since it can be easily seen by induction that the A;-th order derivative of 


y^(a(k)n + b(k))x r ‘ 


n =0 


a(/c)x + b(k)( 1 — x) 
(1 — x ) 2 


is 


/c!(a(/c)(a; + k) + b(k)( 1 — x)) 


_ rr>\k-\~2 


(! ~x) 
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Multiplying by 


x 


x 


-k 

, the differential equation is simplified to 


h^ k \x){ 1 — x) k 
k\ 


a(k)(x + k) + b(k)( 1 — x) 

(1 — x ) 2 


+ (k + l)h(x). 


This differential equation is an equidimensional Cauchy-Euler equation, as the 
one encountered in the previous Chapter. Changing variables x — 1 — z, it is 
h(x) = g( 1 — x). Applying the differential operator 0, where Og(z) = zg'(z), 
the differential equation becomes 


((-i) fc e ( e -1 ) ... (© - k +1 ) - (k + \)\) g {z) 


k\[a{k){ 1 — z + k) + b(k)z) 


and the indicial polynomial Vk{®) is equal to 


V k (&) = (-l) k Qk- ( k + 1)!. 


Using the notation from [26], 0- = 0(0 — 1)... (0 — k + 1) with k > 0, 
denotes the falling factorial. Again, we need a Lemma regarding the roots of 
the indicial polynomial. 

Lemma 4.2.1. The indicial polynomial 7^(0) has k simple roots with real parts 
in the interval [—2, k + 1]. The real roots are —2; (k + 1), if k is even and the 
2 L^J complex roots on,..., a^fe^j with their conjugates ai ,..., . 



108 


Proof. Let a = x + iy be a root of the polynomial. It holds 


a (a — 1)... (a — k + 1) = (— l) k (k + 1)! 

=>• a(a — 1 )... (a — k + 1 ) = (— 2 )(— 3 ). .. {—{k + 1 )) 

aa — 1 a — k + 1 
^ ^2 -3 ~(k+ 1 ) 


(4.5) 


Suppose that 9fc(a) < —2. Then 



V% 2 + y 2 V(-2 ) 2 + y 2 > 

2 2 ~ 

a/(^ - I) 2 + y 2 \/_ (—3) 2 + y 2 

3 3 


> 1 


a — (k — 1) 
-(* + 1 ) 


\J {x — (k — l)) 2 + y 2 \J (~(k + l)) 2 + y 


k + 1 


k + 1 


> 1 . 


Considering the moduli in Eq. (4.5), we see that the left-hand side is a product 
of numbers which are all greater than 1, and so the overall product is greater 
than 1, but the right-hand side is equal to 1, leading to contradiction. There¬ 
fore every root has real part greater than or equal to —2. Further, looking over 
the same argument, we see that the only way we can have the real part being 
equal to —2 is if the imaginary part is equal to zero. 
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By the Fundamental Theorem of Algebra a polynomial of degree n has n com¬ 
plex roots with multiplicities. Note that —2 is always a simple root since, 


r k (- 2) = (-l)‘(-2)(-3)... (-(k + 1)) - (* + 1)! 

= (-l) a (fc + 1)! - (k + 1)! = 0, 


k— 1 .. 

2) = (k + 1)! X -y— = -(fc + l)!(/4+i - 1) < 0. 

j-0 1 

Suppose that a is a repeated root. Since 91c(a;) > —2, we can write a = 
(a: — 2) + iy with x > 0. By the above comments we have that 


k -1 

E 

i=o 

1 


a - 3 


k -1 

V — 

^ (x - 2 — j) + iy 


k -1 


(x — 2 — j) - iy 
j^ 0 {x-2- j) 2 + y 2 


0 

0 

0 


In particular, Jm(o;) = 0. But that imaginary part is equal to ( c _ 9 _^ )2+)/2 

which is clearly only equal to zero if y — 0 i.e. the root is real. We will thus 
have obtained our contradiction if we can show that there are no real roots 
other than —2 and (for k even) k + 1. 

To do this, suppose that we did have a real root a > —2. We note first that 
a > 0: because if not, then since —2 is also a root, there is a root of V' between 
—2 and 0 by Rolle’s Theorem. But since V' k (Ji) = 0 implies that ^ = 0 
and when /3 < 0 this number is clearly negative. Thus any real root a is 
positive. It is also < k + 1 as if it were greater than k + 1 we would have 
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a (a — 1) ... (a — k + 1) > (k + 1)!. Further, note that k + 1 is a root if and 
only if k is even and there cannot be a root in a e [fc, k + 1) as the product 

a (a — 1)... (a — k + 1) would be < (k + 1)!. 

Suppose then that a G (j, j + 1) for some 0 < j <k-i. Then the product 
of the non-negative numbers in the sequence a, a — 1,.a — /c + lisat most 

a (a — 1)... (a — j) < (j + 1)... 2 • 1 = (j + 1)!. Thus, to get a being a root, we 
have to have that 


(—l) A (o; — j — 1) ... (a — k + 1) > yyy — (.j + 2)(j + 3) ... (k + 1). 

However the largest in modulus of a ■=< j — 1,..., a — k +1 is a -- k +1 > j +1 — k 
and so their product is less than (k — 1 — j)\. Consequently our inequalities 
together imply 


(fc + 1)! 

a +1 )!(fc -1 - j)i 


< 



< 1 


and this is a contradiction, completing the proof. 


The differential equation can be written as 


Sk- 1 (©)(© + 2)g(z) 


k\ (a(k)( 1 — z + k) + b{k)z) 


Letting r k = —2 and the remaining (k — 1) simple roots be ri, r 2 ,..., rk- 1 , we 
have 


k\(a(k)( 1 — z + k) + b(k)z) 


(0 - n)... (0 - r fc_i)(0 + 2 )g(z) 
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For the solution of our differential equation, let two functions gi(z) + g 2 (z) = 
g(z). Then 

N n \, , . , .x a(k)(k + 1)! (b(k) —a(k))k\ 

(Q-r 1 )...(Q-r k _ 1 )(Q + 2)(g 1 {z) + g 2 {z)) = -~+ V V ~ 


and by the property of linearity of differential operator 


(0 - n)... (0 - r fc „ 1 )(0 + 2)g 1 (z) 


(0 - n)... (0 - r fc_i_)(0 + 2 )g 2 (z) 


a(k)(k + 1)! 

{b(k)-a{k))k\ 

z 


In the same manner as in the analysis of‘median of (2/c+l)’ Quicksort, applying 
k times the solution, we obtain 


9i{z) 

92{z) 


a(k)(k + 1)! 


l0g e (2) 


(—2 — r i_)(—2 — r 2 )... (—2 — r k -i) z 2 


i= 1 


CjZ 


k\ 


( b{k ) — a(/c)) 


(-l-ri)(-l-r 2 )...l 




i =1 


where q and d* are constants of integration. In order to evaluate i(—2), 
note that 

5*-i(-2) = V t l- 2), 


thus 


,(^2)^1). 


Moreover, 


Pfc(-l) = -**!• 
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Combining both solutions, 


a(k) log e (z) 1 (a(k)- b(k)) „ _ r . 

yy*) TLT 1 -9 1. „ "T" > 


R fc+ 1 -1 


k 


(4.6) 


2=1 


where s* = c t + d t . The constants of integration can be found solving the 
following system of equations 


g(l)= 9 '(l) = ...= 9 <*- 1 >(l)= 0 . 


In terms of series; 


h(x) = 


a(k ) 


oo k 


H, 


k+1 


3T E(( n + - ™))* n + E E s *(- 1 ) r 


n=o 


n=0 i=l 


n 


X 


a(k ) — 6(/c) 
fc 


E 

n =0 


X . 


(4.7) 


The third sum of Eq. (4.7) adds to the solution a constant negligible contribu¬ 
tion. Also, the root (k +1), when k is even, contributes a constant and the root 
r k = —2, adds s fc (n + 1), with E R. Extracting the coefficients, the expected 
cost of multipivot Quicksort is 


Ch<n. 


a(k ) 


LVJ 


Hk+i — 1 


((n + 1 )H n — rd) + s k {n + 1) + 2 ^ ^ (—l) ,! 91c 


2=1 


0 ( 1 ). 


The asymptotics of the last sum can be found by the well-known Stirling’s 
formula, that states [26] 

_ / 77, \ n 

n\ ~ v27rn J ■ 
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Expressing the binomial coefficient in terms of T functions, we have 


(- 1 ) 



—+ n — 1 
n 


r(n - at) 

n\T(-ai) ' 


The relation 6.1.26 in [1] reads for x,y e R. 


|r(a; + ij/)| < |r(x)|. 


thus 


|T(n — «j)| < |T(n — 9 c tc(o'i)) |. 


Using Stirling’s formula, 


T(n — 91c (otj)) 


Tl\ 


2,71 ( n — (91c(«j) + l) 


n — + l) 


n— (lKc(ai)+l) 


V27to 


n 


n — (91c (atj) + l) 


n- (ote(ai)+l) 


'n - (9te(c*i) + 1) ^ ^ (a0+1 ^ f 1 _ 91 c(o,) + 1 


n 


C7n-( OTe ( ai ) +1 ), 


where C = e ( We (";)+i) j s an unimportant constant. 


Therefore the term is bounded by 


2(-l) n 9te^“^ = 0(n" (iHc( “ i)+1) ) 
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and asymptotically, the expected cost is 


a(k ) 

Hk+l ~ 


ynlog e (n) + 



a(k) 

Hk+i ~ 1 


(7 - 1)^ 


n + o(n), 


since all the other roots have real parts greater than —2. 

Knowing the coefficients, any mean cost of the generalisation of the algorithm 
can be directly computed, using this solution, which assumes a simple form. 
These coefficients are related to the number of pivots used during the parti¬ 
tioning scheme. In [28], the average number of comparisons of the first stage 
is given by: 

a(s)n + 0(1), 

where s denotes the number of partitions, when s — 1 pivots are used or equiv¬ 
alently the maximum number of descendants of a node of an s-ary tree. The 
coefficient a(s)n is equal to 

o _ 2 r io S2( s )l 

a(s) = pog 2 (s)l +-. 


Thus, the average number of comparisons of Quicksort on k pivots is 


/ flog 2 (fc + l)] + 1 


2 rio S2 (fc+ i )i 

k +1 


Hk+i — 1 


(n + 1 )H n + 0(n). 
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4.2.1 Derivation of integration constants using Vandermonde 
matrices 


The constants of integration can be found using Vandermonde matrices. Dif¬ 
ferentiating m times Eq. (4.6), 



i(k) (-1 ) m+1 m!((m + 1) log e (,s) - ((m + 1 )H m - m)) 


H, 


k+1 


y m -\-2 


+ (-l + 


kz m+1 


2=1 


(-l) m m! ( a{k)({m + 1) log e (^) - ((m + 1 )H m - m)) 


y m +1 


z(H k+ 1 - 1) 


ci(k) — b(k )) ^ ^ 


m 7 *— 777 , 

srr—^ 1 . 


2=1 


The result can be easily proven by induction or by Leibniz’s product rule. Using 
the initial conditions, namely that g and its first (k — 1) derivatives are 0 when 
evaluated at z = 1, we obtain 




m+1 rn! 


i =1 


a(k)([m + l)Tf m - m)) a(/c) - 6(/c) 


-fffc+i — 1 


k 


for m — 0,1,..., (k — 1). In matrix form, the linear system is 
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Here we use the well-known identity x n = EL, [1], where {”} are the 

Stirling numbers of the second kind, to transform the coefficient matrix into a 
Vandermonde matrix. The determinant of this Vandermonde matrix is equal 
to 

n {r -i - r ') ^ o, 

l<i< 7 <n 

as the roots are all simple. Considering the expected number of passes of 
multipivot Quicksort, it holds that a(k) = 0 and b{k) — 1, for k — 1,2,... . The 
system is, 
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Turning the coefficient matrix into a Vandermonde one and using the identity 
^" =0 (—= (—l) n , (see subsection 24.1.4 in [1]), we obtain 


r 1 1 .. 

1 


Sl 


1 

k 

1 

r\ r 2 

-2 


S 2 


~k 

k-l k- 1 

1 '2 

1 

7 


Sfc 




In [31] and [73] the inverse of Vandermonde matrix is given, in terms of 
product of an upper and lower triangular matrices. Letting A^ 1 denote the 
inverse, it is equal to 


A ^ 1 


1 

i 

1 



r\— r *2 

(ri—r 2 )(n— rz) 



0 

1 

1 



r2— ri 

(r 2 -r 1 )(r 2 -r 3 ) 



0 

0 

1 



(r 3 —ri)(r 3 —r 2 ) 



0 

0 

0 









1 


~r i 


rir 2 


—t ir 2 r 3 


0 

1 

-{n + r 2 ) 
rir 2 + rir 3 + r 2 r 3 


0 

0 

1 

-(n + r 2 + r 3 ) 


It is clear that the lower triangular matrix, post-multiplied by the vector 
(l/k, —1 /k,... (—1 ) k ~ l /k) T , will give us 



3 

, ~(r 1 + l)/k, (n + l)(r 2 + l )/k, - Y[(ri+l)/k, 

i =1 


,,(-l)‘- i p[(r i +l)/fc) T 


i —1 
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Thus, the solution is 


n (o+ x ) 


^ = (- 1 ) 


i¥=i 

k-l l <3< k 


k n 

1 <j<k 


Ti ~ V; 


and the expected number of partitioning stages of multipivot Quicksort on k 
pivots is 


11(0 + !) 

-(n + 1) + o(n). 

k JJf-2 - o) 

j=i 


Note that 


n<r, + D 

J =1 _ 

fc-1 

n <- 2 - »v) 

i=l 


( — 1 ) fc_1 


kk\ 

(k+mH^-i)’ 


therefore the mean number of partitioning stages is 


n+1 

(k + l)(H k+1 -l) + ° {n) - 

We remark that a generalised version of this result can be found in [28]. In 
[35], it was shown that the constants of integration can be computed for ar¬ 
bitrary values of the coefficients a(k) and b{k) by the same method, as in the 
derivation of the integration constants in the simple case of a(k) = 0 and 
b(k) = 1. 
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At the end of this section, it should be noted that the worst-case probability is 
not eliminated, but is less likely to occur. In an unfortunate situation, where 
the k smallest or greatest keys are selected as pivots, partitioning will yield 
trivial subarrays and one containing the remaining elements. If the chosen 
pivots happen to be close to the quantiles of the array, this yields an optimal 
partitioning of the array. In the next section, we examine ways of a more 
efficient selection of pivots. 


4.3 Multipivot-median partitioning 

The preceding analysis of multipivot Quicksort, where k pivots are uniformly 
selected at random has showed that the worst-case scenario is less likely from 
the standard ‘one-pivot’ model. Is any other way, where we can reduce further 
the probability of such scenario? We have seen that choosing the median from 
a random sample of the array to be sorted, yields savings to the running time 
of the algorithm. Since we have examined the analysis of multiple pivots, then 
we can select these pivots as the quantiles of a bigger random sample. 

Thus, one can randomly choose a larger sample of k{t + 1) — 1 keys, find the 
(t + l)-st, 2 (t + l)-th, . .., (k — l)(t + l)-th smallest keys and use these ( k — 1) 
statistics as pivots. Note that for k = 2, this variant contains the median 
of 2t + 1 Quicksort as a special case and for t = 0, we have the multipivot 
algorithm, whose mathematical analysis was presented in the previous section. 
This ‘generalised Quicksort’ was introduced by Hennequin [28]. Let T(n^,t}) 
be the average of a “toll function” during the first pass and f(n{ k ,t}) the total 
expected cost of this variant, when applied to an array of n keys. The following 
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since the pivots i\,..., 4-i ate selected to be the (t + l)-st, ..., (k — 1)(£ + l)-th 
smallest keys of the sample and each of the k resulting subarrays (— 1), (i 2 — 
ii — 1),... j,(n — ik~i) contain t elements of the sample. 

As before, the general recurrence of average cost is translated to a differential 
equation with indicial polynomial [28], 

T W .(e> = (-i)^ ,) - 1 (*<,+* _,)-M-i)'(!)■ 

A Lemma follows concerning the whereabouts of the roots of this polyno¬ 
mial: 

Lemma 4.3.1. The indicial polynomial Vkp+ i)-i(@) has k{t + 1) — 1 simple 
roots, with real parts greater than or equal to —2. The real roots are the integers 
0,1— 1), —2; k{t + 1) + t, when t is odd and k is even or when t is even 
and k is odd and the 2 complex roots Ai,..., A y (k-nt+k- 2 j with their 

conjugates Ai,..., A ^ (k~i)t+k ~2 j. 
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Proof. It can be easily deduced that the integers 0,1, ..., (t — 1) and —2 are 
roots of the polynomial. Now, when t is odd and k is even, then k{t + 1) — 1 is 
odd and so a root a of the polynomial will satisfy 


U + ViM“) (4 - 8) 

and by simple manipulations we can now verify that k(t + 1) +1 is also a root. 
Similarly, if t is even and k is odd, we have that k{t + 1) — 1 is even and Eq. 
(4.8) is valid, again making k(t + l)+ta root. 

It will be proved by contradiction that all roots have real parts greater than 
or equal to —2. Note that the argument is similar with the proofs of Lemmas 
3.1.1 and 4.2.1. For any root r = x + iy, with x,y E R holds 






(4.9) 


For r / 0. 1,_ :(l — 1), Eq. (4.9) can be written as 


/ T xfe(t+n—i (r — t)(r — t — k(t + l) + 2) 

( j (t + l)...(*(t + l)-l) ( } 


(4.10) 
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Assume that 9te(r) < —2, then 


r — t 
t + 1 


y/{x - t) 2 + y 2 
t + 1 


\J (-0 + 2 )) 2 + y 2 ^t + 2 

7TT _ t +1 


r — t — 1 
t + 2 


v/Qk - t - l) 2 + y 2 
t + 2 


\/ (—(^ + 3)) 2 + y 2 ^ t + 3 
7T2 “ t + 2 


r — (fc(t + 1) — 2) 
k(t + 1) - 1 


^(x ~ k(t + l)+ 2) 2 +y 2 ^ yj(-k(t + l)) 2 + y 2 
k(t + 1 ) - 1 > k(t + 1 ) - 1 

> + 1) 

— &(£ + 1) — 1 


Considering the product of moduli, we see that the left-hand side of Eq. (4.10) 
is greater than or equal to (in modulus) the telescoping product k(t + 1)/(t + 
1) = k, which gives a contradiction. Further, this argument shows that —2 is 
the unique root with the least real part. 


We now show the roots are simple. Assuming, for a contradiction, that r is a 
repeated root, we get: 


- 1 ) 


fc(t+i)-i 


fc(t+l)-2 / A t-1 i 


k(t + l)-l; “ r-j 


t J ' r — j 
7 i =o J 


(4.11) 


Eq. (4.9) and (4.11) imply that 


k(t+ 1)-2 t- 1 

E — = E— 

f ^ y n f ^ y i 

3=0 J 3=0 J 


or 


k{t+ 1)-2 

E 


j=t 


1 

r-j 


0. 


(4.12) 
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From Eq. (4.12), we deduce that 3m(r) = 0 and re (t, t + 1) U {t + 1, t + 2) U 
. .. U ( k(t + 1) — 3, k(t + 1) — 2). However, the modulus in the left-hand side of 
Eq. (4.9) is smaller than 1, while the right-hand side is greater than k, proving 
that all roots are simple. ■ 

By the Lemma, the polynomial can be written in terms of simple factors and 
the differential equation can be solved using the same way, as in other variants 
of Quicksort, previously analysed. The average cost of ‘generalised Quicksort’ 
is 


a(k, t) 

Hk(t+ 1 ) — H t+ 1 


((n+ 1 )H n 


- n ) + 0(n), 


when the “toll function” is linear and its average is a(k, t)n + 0(1). 

Our analyses of the average cost of the algorithm and its variants has showed 
that any Quicksort needs on average Qnlog e (n) + 0{n) key comparisons for 
the complete sorting of a file consisting of n distinct keys. The constant C can 
be made very close to the information-theoretic bound, as we saw. In many 
sorting applications the ‘median of 3’ is being used, with savings on the average 
time and little overhead for the computation of median. For large arrays, one 
can use the ‘remedian of 3 2 ’ Quicksort. 



Chapter 5 


Partial order of keys 

5.1 Introduction 

Here, we investigate the analysis of Quicksort under the assumption of prior 
information of the order of keys. Specifically, we assume that there is a partial 
order on the keys. The rough idea is to see how much having partial infor¬ 
mation compatible with the true order allows us to speed up the process of 
finding the true order. 

Let us illustrate the idea first with a simple example. Suppose that there are 
d levels with k keys at each level, so that n = kd. Anything in a higher level 
is known to be above everything in a lower level. Computing the ratio of the 
expected complexities, we have 

E(C n ) d(2(k + 1 )H k - 4 k) 

E(Q) “ 2(n + l)H n — An ’ 

where C* and C n denote the number of comparisons of Quicksort with uni¬ 
form pivot selection and in case of partial order, respectively. We consider the 
following cases: 
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1. When d is fixed number, then as n tends to infinity, 


lim ^ = oo. 

n—>-oo d 


Thus, 


lim ~ lim d j21pM = llm ^ = lim = !. 

n-xx> E(C*) n-Hx> 2nlog e (n) «->oo log e (n) n ->oo log e (n) 

2. d — k — \fn. 

li m _ lim Vn(2(y/n + 1 )H s/Tl - Ayjn) ^ ^ 2nlog e (y/n) _ 1 

n—>oo E(C7*) ih oo 2(n + l)H n — An n-+oo 2nlog e (n) 2 

We see that Quicksort is on average twice as fast, when sorting a partially 
ordered array 

3. k — - where c is a constant. Then, 

lim E(C n ) _ cn{2{\ + l)H 1/c - |) _ w((2 + 2c)H l/c - 4) _ d 
n-»-ooE(C'*) 2{n+l)H n — An 2nlog e (n) log e (n)’ 


(2 + 2c)Hi/ c — 4 

where d =- - -. 

2 

We should think a little about variability too. It is unsurprising that having 
the additional information about levels reduces variability of the number of 
comparisons, let us get a preliminary result. If we have the level structure, 
then Var(Q*) is the sum of the variances of sorting each of the d independent 
levels. Each of these variances, since there are k keys in each level, is just 
Vax(C'fc). For simplicity, we assume k —> oo as n —> oo and do asymptotics. We 
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then have (for m either n or k) 


Var (C m ) ~ 



• m 


2 


and thus we get 


Var (C n ) (7 — 2tt 2 /3)n 2 
Var (C*) (7 — 2-k 2 /?>)k 2 d 

so the variance of the version with the presorting is reduced by a factor of 
about d. These suggest there is interest in studying this situation, we now do 
so in more detail. 


5.2 Partially ordered sets 

An approach of having additional information is the partial order of the keys. 
We shall employ this assumption along the following lines. First, we present a 
definition [56], 

Definition 5.2.1. Let a finite set P equipped with a binary relation ‘<’ which 
has the following properties. (Here, x, y and z are elements of P). 

(i) x < x, Vx E P. (That is, < is reflexive) 

(ii) If x <y and y < x, then x = y. (< is antisymmetric) 

(iii) If x <y and y < z, then x < z. (< is transitive) 

Then the pair (P, <) is called Partially Ordered Set. 

Henceforth, in this thesis we abbreviate ‘partially ordered set’ to ‘poset’. We 
also present two key definitions [56]. 
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Definition 5.2.2. Let (P, <) be a poset. We say that two elements x and y of this 
poset are comparable if x < y or y < x. Otherwise they are incomparable. 

Definition 5.2.3. Let (P, <) be a poset. 

(i) A minimal element of (P, <), is an element with the property that no other 
element is smaller than it. A maximal element of (P, <), is an element with the 
property that no other element is greater than it. 

(ii) A chain in P is a set T of elements, where every pair of elements of T are 
comparable. The number of elements of P in the longest chain in P is called the 
height of P and denoted by h(P). 

(Hi) An antichain in P is a set U of elements, no two of which are comparable. 
The number of elements of P in the order of the largest antichain is called the 
width of P, and is denoted by w(P). 

(iv) A total order in P is a partial order where every pair of elements are compa¬ 
rable. 

For example, the set of subsets of X = {1,2} has an antichain of order 2, 
namely {1} and {2}. A chain of length 3 in it, is 0 < {1} < {1,2}. This partial 
order is not a total order as {1} and {2} are not comparable. Usually, if we 
have a partial order on a set, there will be several ways of extending it to a 
total order on that set. Often, we will use -< rather than < to denote the partial 
order. Further, we present the following definition, that we will come across 
later. 

Definition 5.2.4. Let a poset (P, <J. Its comparability graph G(P ) is the graph 
with the poset’s vertex set, such that the elements are adjacent if and only if they 
are comparable in (P.\ <). Its incomparability graph G(P ) is the graph, such that 
the elements are adjacent if and only if they are incomparable in (P, <). 
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Another example of a partial order which usually is not a total order is the 
collection of subsets of a fixed set X, with the partial order < being inclusion, 
normally denoted as C. It is easy to check that, for any A C X, we have that 
A < A since any set is a subset of itself: if A < B < A then we indeed have 
that A = B, giving asymmetry: and finally, if A C B C C then of course ACC 
and so < will be transitive. This is not a total order if A" has order at least 2, as 
{xi} C X and {x 2 } are not comparable for x\ f x 2 members of X. However, 
when we have a partial order on a set P there will be at least one total order 
on P extending it, and in fact usually there will be several such: 

Definition 5.2.5. A linear extension of a partial order (P, -<) is a total order < 
on the set P such that whenever x -< y in the partial order, then we have x < y 
in the total order too. The number of linear extensions of a poset P is denoted by 

e(P). 

In other words, a linear extension of a partial order is a total order on the same 
set which is compatible with the partial order. This is of course of great rele¬ 
vance to us, as the situation we are in is that we are given partial information 
on the true order of the set of elements and want to know how many more 
pairwise comparisons we have to do to work out the true order on it: that is, 
we are trying to identify which of the numerous linear extensions of the partial 
order is the true order on it, with as few comparisons as possible. 

The number of linear extensions of a poset can vary substantially according to 
the structure of the poset. For example, trivially, if the partial order happens 
already to be a total order there is only one extension, namely itself. Equally 
trivially, if the partial order contains no comparisons - i.e. it provides no 
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information whatsoever - then all n\ possible orderings of the n elements of P 
are linear extensions. 

Here is a generic lower bound on the number of pairwise comparisons we need 
to make in order to find the true order of our data, given a partial ordering P 
of it. 

Theorem 5.2.6. Given a partial order (P, -<) which is partial information about 
the true total order on the underlying set P, it takes at least |~log 2 (e(P))] pah~wise 
comparisons to find the total order. 

Proof. Recall that linear extension or total order of a poset (P, -<) is a total 
order compatible with the partial one. All elements are comparable, forming 
a unique chain. Now let x, y be incomparable, distinct members of poset. In 
A total orders, we will have that x < y and in B linear extensions, x > y. It 
holds that A + B = e(P) so, 


max(A, B) > ^ . 

Thus, after making one comparison, there are at least —p- candidates. Simi¬ 
larly, of those number of linear extensions, choosing again two elements, we 
have to consider at least Ap- total orders. Thus, after r comparisons, we will 
examine at least c 4p linear extensions. We want to identify the unique total 
order among the number e(P) of all possible linear extensions. Therefore, this 
fraction becomes equal to unity when r — [log 2 (e(P))] comparisons. ■ 

Remark 5.2.7. The quantity log 2 (e(P)) is the information-theoretic lower bound, 
that we first came across in section 3.1. 
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Kislitsyn [41] and independently Fredman [25], showed that often this lower 
bound is close to the truth. 

Theorem 5.2.8 (Fredman [25], Kislitsyn [41]). Sorting an array of n keys, 
which obeying a partial order P, can be achieved in worst case by log 2 (e(P)) +2 n 
comparisons. 

Corollary 5.2.9. If P is a poset for which log 2 (e(P)) grows faster than n, we 
have that the number of comparisons of finding the true total order, in worst 
case, is log 2 (e(P)) (l + o(l)). 

Proof. We have that this quantity is bounded below by log 2 (e(P)) and above 
by log 2 (e(P)) + 2 n which is log 2 (e(P))(l + o(l)) by the assumption in the 
statement of the corollary ■ 

Kahn and Kim [40] gave an algorithm for actually doing the finding of the 
true total order, which uses at most 54.45 log 2 (e(P)) comparisons. This has 
recently been reproved by Cardinal et al. [13] whose proof manages to avoid 
Kahn and Kim’s use of the ellipsoid method, a technique which though it is in 
theory polynomial-time, is difficult to do in practice. What all this makes clear 
is that, in considering how much information we can deduce from a random 
partial order, we will need to know about the logarithm of the number of linear 
extensions the partial order typically has. 

Often a useful notion in studying posets is the theory of levels, [10], [42]. 
Definition 5.2.10. If (P, <) is a poset, we define 


Li = {x G P : fly G P, y < x A y f x} 
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the set of minimal elements of our poset to be the first level of our poset. The next 
level is the set of minimal elements in P \ L x : each of these will have (at least 
one) element of Li below it. We then continue by induction, defining Li to be the 
level of minimal elements of P\ (iTflfLf). 

Note that every level of a poset is an antichain: for two minimal elements in a 
poset cannot be comparable with each other. Moreover, every time you go up 
in a chain, you go up to a higher level. Thus the height of the poset will be the 
number of levels. 

Definition 5.2.11. The linear sum of two posets (Pi, -<i) and (P 2 , -< 2 ) is a poset 
with vertex set the disjoint union of Pi and P 2 and with x -< y if and only if: 

(i) if x and y are in P\ and x -<y y: 

(u) if x and y are in P 2 and x -< 2 y: 

(Hi) if x E P\ and y e P 2 , then automatically we have x -< y. 

A useful definition for us will be the following [56]. 

Definition 5.2.12. 

(i) Suppose P is a partially ordered set, and L is a particular total order on the 
same set which agrees with P on every pair of elements which are comparable in 
P. (In other words, L is one of the linear extensions of our partial order). Then a 
setup is a pair x, y of elements which are incomparable in P but are consecutive 
in L. 

(ii) The number of setups S(P, L ) which must be made comparable to obtain the 
linear extension L is denoted s(P, L). 

(iii) The setup number of P is the minimum, over all the linear extensions L of 
P, of s(P, L) and it is denoted s(P). 
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So the point is that, if information is given in the partial order P and using 
pairwise comparisons to obtain the rest of the order, s(P) is a lower bound on 
the number of comparisons which have to be done to find the true order. 

The setup number is also known as the jump number. The following Lemma 
is used for the derivation of a simple lower bound on the setup number of a 
poset. 

Lemma 5.2.13 (Dilworth [17]). The minimum number of chains into which a 
poset P can be partitioned is the width w(P). 

Proof. This is standard and there are several proofs, we refer to [17]. ■ 

Theorem 5.2.14. For any poset P, s(P) > w(P ) — 1. 

Proof. The best case is to partition the poset into w(P) chains C), C' 2 ,..., C W ( P ) 
and hope that there is some ordering of these chains (which without loss of 
generality is the order given) such that every element in is less than the 
minimum element of C i+ 1 for each 1 < i < w — 1. Because if this happens, then 
the total order is just the direct sum of the C) and we only have to do w(P ) — 1 
comparisons of the maximum element of C, with the minimum element of 
C'+i- ■ 

What we will do for the next while is consider various ways in which we could 
have a partial order given to us before we start using Quicksort to determine 
the complete order. So we are imagining that a previous researcher had carried 
out some of the comparisons and we want to know how many more compar¬ 
isons we have to carry out to determine the total order. We will consider cases 
where the poset is randomly generated. 
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5.3 Uniform random partial orders 

Definition 5.3.1. A uniform random partial order is a partial order selected 
uniformly at random from all the partial orders on S — {1,2,..., n}. 

This means that all partial orders on S are equally likely to be chosen. The ba¬ 
sic structural result on such posets is the following, rather surprising, one. 

Theorem 5.3.2 (Kleitman and Rothschild [42], Alternative proof by Brightwell, 
Promel and Steger [11]). Suppose that <is a uniform random partial order on 
{1,2,, n}. Then, whp. (henceforth, ‘whp.’ stands for ‘with high probability’, 
which denotes the fact that as n —> oo, the probability of an event approaches 1) 
there are three levels: the bottom level has approximately n /4 elements in it, the 
middle layer approximately n/2 elements and the top layer about n /4 elements 
in it. 

Remark 5.3.3. This feature of a uniform random partial order - that it has 
height only 3 - is surprising to most mathematicians when they hear it, and 
perhaps suggests that “in nature” posets do not occur uniformly at random - 
some posets are favoured over others. 

Here is the key information on the number of linear extensions. 

Theorem 5.3.4 (Brightwell [10]). Given any function u(n) tending to infinity 
with n (think of it as doing so extremely slowly) the number of linear extensions 
of a random partial order chosen uniformly at random is whp., between 

(n/2)!((n/4)!) 2 ( ^ 

-y—- — and (n/2)!((n/4)!) uj(n). 

oj(n) 
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Proof. We refer to Brightwell [10], page 66. ■ 

Corollary 5.3.5. When P is selected uniformly at random, we have whp. 


log e (2) 


nl°g e (n) - 0(nj) < log 2 (e(P)) 


< 


log e (2) 


nlog e (n) + 0(n)). 


In other words, in this situation, the entropy lower bound on the number of 
comparisons required is essentially the right answer. 


Proof. We have by the previous result, taking u (n) to go to infinity very slowly, 
in particular more slowly than log e (n), that 


lo g 2 ( j - lo S 2 (e(P)) < log 2 ^(n/2)!((n/4)!) 2 o;(n; 

1 lo g e 1 < l°g 2 (e(P)) < lQ g e ( (»/2)!(W 4 ) ! )w(n) 


!og e (2) 

1 


l°g e (2) 

< log 2 (e(P)) 


w(n) J - v " ~ log e (2) 

((n/2) log e (n/2) + 2(n/4) log e (n/4) - (n/2) - 2(n/4) + 0(log e (n)). 


< 


l°g e (2) 
1 


((n/2) log e (n/2) + 2(n/4) log e (n/4) - (n/2) - 2(n/4) + 0(log e (n)). 


log e (2) 


(nlog e (n) - 0(n)) < log 2 (e(P)) < 


loge(2) 


nloge(n) + O(n)), 


as required. 


Using this Corollary, a key result follows regarding algorithm’s time complex¬ 
ity: 
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Corollary 5.3.6. Given a uniform partial order on a set of n keys, the time taken 
to sort them by pairwise comparisons is approximately 


1 

21og e (2) 


« 0.72135 


times the number of comparisons required by Quicksort to sort them without the 
partial information. 


Proof. The expected number of comparisons required by Quicksort for the 
sorting of n keys is 2 n log e (n) (l + o(l)) , and since the variance of this is asymp¬ 
totically (7 — 27r 2 /3)n 2 (l + o(l)), we have by Chebyshev’s inequality, letting C n 
be the number of comparisons 


P(|(C' n -E(C n ))| > nlog e (log e (n))) < 


(7-2vr 2 /3)n 2 (l + o(l)) 


n 2 (log e (log e (n))) 2 
\C n - 2nlog e (n)(l + o(l))| > nlog e (log e (n))) ->• 0, as n -» oo. 


Thus the probability of the complementary event, namely that C n is within 
n log e log e (n) of its mean will tend to 1. Therefore, whp. the number of com¬ 
parisons is in 

(2nlog e (n)(l + o(l)) - nlog e (log e (n)), 2nlog e (n)(l + o(l)) + n log e (log e (n))). 

Thus, whp. it takes about 


2 n log e (n) (l + o(l)) comparisons. 
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On the other hand, we have just seen that with a uniform partial order of keys, 
it takes about 

-—— • n log. (n) comparisons. 

lo g e ( 2 ) 

This number is indeed 1/(2 log e (2)) times the number Quicksort needs and the 
numerical value of this fraction is as stated. ■ 

It is of interest to compare this with the naive lower bound on setup number, 
which performs rather poorly here. 

Corollary 5.3.7. The setup number of a uniform random poset is at least 

Proof. This is immediate from the lower bound s(P) > w(P) — 1 and the fact 
that, by the Theorem of Kleitman and Rothschild [42], we clearly have that 
w(P ) is greater than n/2 + 0(1). Thus in this case the simplest setup number 
lower bound is not a very good one, as we have seen that the true answer is 
0(nlog(n)) in this model. ■ 

This in turn implies that if we were to use Quicksort, even in the optimal cases, 
we would have to compare n/2 + 0(1) pairs of keys. The expected time to do 
this would be asymptotically 2{n/2) log(n/2) ~ n log(n). 
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5.4 Random bipartite orders 

In this section, we inspect the number of linear extensions for bipartite orders. 
Let present a definition [10]. 

Definition 5.4.1. Let X and Y be two disjoint sets, each one having cardinality 
equal to n. A random bipartite order A P (X, Y) is the poset X U Y with both X 
and Y antichains, and for each pair (x, y) G X x Y there is a relation x -< y with 
probability p and no relation (i.e. they are incomparable) with probability 1 —p, 
independently of all other pairs. 

We now think about the number of linear extensions. A linear extension of 
such an order will have to, amongst other things, put the set X in order - there 
are nl ways to do this - and there are similarly n! ways to order Y. However 
there will also be some choices to make elsewhere, because while we have 
some relations x < y for (x, y) G X x Y in the partial order, we will also have 
some incomparable pairs. More precisely, the probability that a total order on 
X, a total order on Y and a decision rule a for each pair (x, y) G X x Y on 
whether x -< y or y -< x, is a total order compatible with the partial order - i.e. 
a linear extension of the partial order - is (1 — p)^ a \ where £(a) is the number 
of reversals in a, that is the number of pairs (x, y) G X x Y such that x > y in 
the total order. 

Thus the expected number of linear extensions is (n!) 2 Ed — p)^ a \ The func- 

Oi 

tion multiplying (n!) 2 here is discussed at length in [10]: it is defined 

OO 

v(p) ■■= n( x “ ( X -pt)- 

i=l 
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Though some more work needs to be done to check this, it turns out that whp. 
the number of linear extensions is close to this mean value. We quote the 
result from Brightwell [10]. 

Theorem 5.4.2 (Brightwell [10]). The number of linear extensions e(P) of a 
random bipartite partial order A p (X,Y) with \X\ = |Y| = n and probability p 

p(n ) • n 1 /' 

such that lim -— = oo, satisfies whp., 

(log 2 (n)) 1/7 

e(P) = < " !)2 '^' (l + 0(1)) - 

In particular this applies when 0 < p < 1 is a constant. More precisely, we 
have that whp. 

(»!f ■ ■ (l - tbMM) < e(P) < («!)=. v( P r ' (l + . 

Thus, for p constant, letting C = log 2 {v(p)) an d noting that both logarithms 
approach 1 as n —>■ oo, we have that whp. 

2log 2 (n!) -(7+0(1) < log 2 (e(P)) < 21og 2 (n!) -C + o(l). 

(Of course the two o(l) terms are different). Thus the order of magnitude of 

log 2 (e(P)) is whp. 

oi „ , n 2 • log e (n!) 2nlog e (n)(l + o(l)) 
log e (2) " log e (2) 


We need to be careful about comparing this example with Quicksort: we 
must remember that the total number of keys being sorted in this example 
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is n + n — 2 n. Therefore, the expected time would be An log e (2n) (l + o(l)) = 
4nlog e (n)(l + o(l)). Thus the factor by which we are quicker here is again 
2 log e (2). In other words, we get the same speed-up as for the uniform and bi¬ 
partite cases. We now move forward to the analysis of random A;-dimensional 
orders in Quicksort. 

5.5 Random k-dimensional orders 

Here, the application of Quicksort in a random k-dimensional order is consid¬ 
ered. For this purpose, a definition follows: 

Definition 5.5.1. A random k-dimensional partial order on the set P = {1,2,..., n} 
is defined as follows. We select k total orders on {1,2,..., n} uniformly at ran¬ 
dom from all n\ total orders on that set, say we chose <i, < 2 ,..., <k■ We then 
define the partial order -< by 

x Ay x <i y for all 1 < % < k. 

Of course it is highly likely that some of the total orders will be inconsistent 
with each other, and so we will only get a partial order. We now have to change 
perspective: we assume that the partial order which results from these k total 
orders is given to us as partial information about the order on the set P, and 
that we have to use pairwise comparisons to find the true total order on P. 
Therefore, the total order we are looking for, and the k total orders we used to 
define the partial order, may have little to do with each other. 
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We aim to estimate time complexity of finding the true order when the partial 
information given is a random /c-dimensional order. Again, we use the Theo¬ 
rem 5.2.8 coupled with the information lower bound. Thus we need to know 
about the number of linear extensions of a random /c-dimensional order. An 
important result follows, 

Theorem 5.5.2 (Brightwell [12]). The number e{P) of linear extensions of a 
random k-dimensional partial order P with \P\ — n satisfies whp. 

(e-V- 1 ^)” < e(P) < (2kn 1 ~ 1 ' k ) n . 

Consequently we have that log 2 (e(P)) is bounded below by n ■ (1 — 1/k) ■ 
log 2 (n)(l + o(l)) and similarly is bounded above by n ■ (l°g 2 (2 k) + (1 - 1/k) ■ 
log 2 (n))(l + o(l)). This is of course larger in order of magnitude (for fixed 
k, say) than 2 n so the information lower bound is tight. Now, we deduce 
that 

Corollary 5.5.3. The time complexity of finding the true total order on a set 
given a random k-dimensional partial order on it, where k is a constant, is whp. 
asymptotically equivalent to 

rclog e (n) L _ 1\ 

l°ge(2) V k ) 
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Proof. For the logarithm of the number of the linear extensions as n tends to 
infinity holds, 


n • ( 1 - • log 2 (n)(l + 0(1)) < log 2 (e(P)) 

< n • (log 2 (2 k) + (l - ' lo S 2 ( n )) 

* ■ (* ■" s)' lofc(,!) 


< log 2 (e(P)) < 


n 


108 .( 2 ) 


log e (n)). 


Then, we obtain 




which completes the proof. ■ 

Therefore, the speed up relating to Quicksort with no prior information is on 
average 


n • log e (n) 


1 - 


2n-log e (n)log e (2) k 

1 / i 

2 log e (2) \ k 

1 


0.72135 • 1 


k 


Remark 5.5.4. Note that the factor 


2 l°g e (2) 


0.72135. 
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was previously encountered in uniform and bipartite random orders. It is worth 
to point out that for the multiplier 


1 


1 


k' 


when k is arbitrarily large, the speed up is as of the case of uniform random 
orders. Whereas, having few k-dimensional orders, Quicksort runs much faster. 


5.6 Random interval orders 

In this section, we examine the case where poset forms an interval order. A 
definition follows [23], [71]. 

Definition 5.6.1. A poset (P, <) is called an interval order if there exists a func¬ 
tion I such that each element x e P is mapped to a closed interval I(x) = 
[a x , b x \ C R. Then \/x, y E P, it holds that x -< y if and only ifb x < a y . 

In other words, there exists a mapping x \-E I(x) [a x , b x ] for every element of 
(P, <), having the property that any two elements of the poset are comparable 
if and only if their corresponding intervals do not intersect. Otherwise, they 
are incomparable. Hence, the size of the largest chain of the poset is the 
maximum number of pairwise non-intersecting intervals. Conversely, the size 
of the largest antichain is the maximum number of intersecting intervals. We 
present the definition of random interval order. 

Definition 5.6.2. A random internal order is one where we generate 2 n inde¬ 
pendent numbers Xi ,..., X n , Y u ..., Y n from the uniform distribution on [0,1] 
and form n closed intervals I,, for 1 < j < n, where I :j = [X,-, Yf if Xj < Yj 
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and [Yj,Xj] otherwise. (The event that X t = Yj has probability zero so can be 
ignored). Then we define a partial order by saying that h I I if and only if the 
maximum element of f is less than the minimum element of I.j. 

Remark 5.6.3. In fact any continuous probability distribution can be chosen to 
the analysis of random interval orders. 

Again, we want to estimate how many linear extensions there are of these. 
This time, there does not appear to be an immediate bound for the number of 
linear extensions in the literature. However one can obtain the relevant bound 
showing that log 2 (e(P)) is whp. at least cn log e (n) for some c > 0, which will 
of course be enough to show that the 2 n term in the Fredman [25] - Kislitsyn 
[41] bound log 2 (e(P)) + 2 n is small compared with the term log 2 (e(P)). The 
main results that we need in this direction are the following two Theorems, 
regarding the size of the largest antichain and chain of a random interval order 
respectively. For their proofs, we refer to [39]. 

Theorem 5.6.4 (Justicz et al. [39]). Let A n denote the size of a largest set of 
pairwise intersecting intervals in a family ofn random intervals. Then there exists 
a function f(n ) = o(n), such that whp. we have 

2 ~ f ( n ) < A n < 2 + /( n )- 

Theorem 5.6.5 (Justicz et al. [39]). Let Y n denote the maximum number of 
pairwise disjoint inteixals in a family ofn random intervals. Then 



in probability. 
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The following Corollary gives a lower bound for the number of comparisons 
required to sort a random interval order: 

Corollary 5.6.6. The number of comparisons for sorting n keys, given a random 
interval order is whp. at least cri log e (n) ( 1 + o(l)), where one can take c = 
1/2 log e (2) ft* 0.72135. 


Proof. Theorem 5.6.4 shows that the largest antichain of the random interval 
order is whp. at least r — [(1 — e)n/ 2] for any e > 0. This is because a family 
of intersecting intervals forms an antichain. Thus we need to sort all these r 
incomparable elements of the partial order in a total order extending it, and 
there are at least r! ways of doing this. Using Stirling’s formula, we obtain 

=► l°g 2 (e( p )) > log 2 (r!) > ^ 2 e ^ lo g 2 2^™ ) + ^ log2 (^ “ e ) n7r )’ 


Then we have that log 2 (e(P)) has order of magnitude nlog 2 (n): in particular 
2 n = o(log 2 (e(P)) and so the complexity is log 2 (e(P)) (l + o(l)). Further, we 
have that log 2 (e(P)) is at least, by the above, n{ 1 — e) log 2 (n)(l + o(l))/2 which 
is equal to n log e (n) ^ °^ and so we can take the constant c to be at least 


2 log,(2) 


1/2 l°g e (2). 


We now present a much stronger result that gives sharp bounds on the num¬ 
ber of linear extensions of a random interval order, following the insightful 
suggestions of Prof. Colin McDiarmid [52], 
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Theorem 5.6.7. The number of comparisons for sorting n keys, given a random 
interval order is whp. for 0 < e < 1, between: 

(l-e + o(l))nlog 2 (ra) < log 2 (e(P)) < (l + e + o(l))nlog 2 (n). 

Proof. Let 0 < a < b < 1 and consider the interval (a, b). Let I(i,j) be the 
interval (^4, where i and j are positive integers with i odd and i < 2 j . 
Further, let j(a, b) be the least j such that jj e (a, b ) for some positive integer 

i. There is a unique such i since if there were at least two odd i, then there 
is at least one even k between them and considering k/2 j , we can replace it 
by (k/2)/2 j ~ 1 giving a smaller value of j and contradicting the definition of j. 
Since i is unique, we may call it i(a, b ). 

We denote the interval l(i(a,b),j(a,b)) by J(a, 6 ). Recall that i is odd with 
i < 2 j and observe that if b — a > 2~ j , then j(a, b) < j. For a given such i and 

j, let A{i, j) be the set of all intervals (a, b), such that J(a, b) = The sets 

A(i,j) are antichains: for if we had two intervals (ai, bf) and (a 2 , bf) in A(i,j) 
with (ai, 6 i) being less than ( 02 , 62 ) (of course this is equivalent to b\ < a 2 ) 
then saying that 


61 ) — J(a 2 , 6 2 ) — 

would imply that j(a 1 ,b 1 ) = j(a 2 ,b 2 ) and i(ai,bi) = i(a 2 ,b 2 ). But given that 
j(ai,6i) = j(a 2 , b 2 ), there is clearly some if2 j with i odd in (ai,6i) which is 
less than any such in (a 2 , b 2 ) and the result follows. Indeed, if the midpoint of 
I(i,j) (i.e. i/2 j ) is less than the midpoint of /(*',/) (i.e. i'/2 j ') then no interval 
in A(i',f) can precede any interval in A(i,j) in the interval order. 


fi- 1 i + 1\ 
V 2i ’ 2J ) 
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Let X and Y be independent and uniformly distributed random variables 
which denote the endpoints of a random interval in (0,1). We have that 

P(|X — Y\ < e/4) < e/2, 


since 


p(|*-y|< j) = r ( y -\< x<y + ^ 

= / F(Y- e -<X<Y+ e -\Y = y)f Y (y)dy 

J 0 

r i 

= / F(Y-^<X<Y + ^\Y = y)dy. 

J 0 

The last equation follows because the probability density of Y, fy(y) is 1 on 
[0,1] and 0 elsewhere. Given that Y = y, the probability that X is in the 
interval (y — e/4, y + e/4) is (as it is uniformly distributed) at most the length 
of the interval (y + e/4) — (y — e/4) = e/2. 

Thus the random number J\f of intervals with length at most e/4 is stochasti¬ 
cally dominated by a binomial random variable B(n , e/2) with n independent 
trials and success probability e/2. Therefore, 

P(A/” > en) < P (13 (n, e/2) > en) 


and this probability tends to 0, as n —> oo. This is a consequence of Chernoff’s 
inequality, in the following form: if A" is a binomially distributed variable with 
n independent trials and success probability p, then for 5 > 0 



p+s f i-p y-p~ s \ 

\l-p-5j ) 


¥(X > n(p + 8 )) < 
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Chernoff’s inequality appears in various places: we refer to [8] and to [54]. 
Since the number being raised to the power n is < 1, this will indeed tend to 0 
(in fact will do so rapidly). The result in our case follows plugging in P = e/2 
and 6 = e/2 and shows that the number of intervals with length at most e/4 is 
whp. less than en. 

Let j 0 be a positive integer and m be the number of intervals I(i, j) with j < Jo- 
Let X be a set of intervals and X' be the set of intervals in X with length > 2~ jo . 
Let n(i,j ) be the number of intervals of X' in A(i,j). Then, the number of 
linear extensions e(P ) satisfies, since we have to put each of the antichains in 
order and there are r! ways to order an antichain of r elements, 


e (P) > > J J 




n(i,j) 


and by convexity 

l°g 2 (e(P)) > J ) log 2 > |Xj log 2 —. 

Choosing j 0 sufficiently large that 2 _J0 < e/4 , the number of intervals of 
length at most 2”- 70 is less than the number of intervals of length < e/4 which 
by Chernoff’s inequality is whp. < en: thus almost all intervals have length at 
least 2 -J0 and so T has order at least (1 — e)n. This will give us that whp. 

log 2 (e(P)) > (1 - e + o(l))nlog 2 (n). 

This completes the proof, as this gives the lower bound and the upper bound 
is just a consequence of the fact that there are at most n\ linear extensions of a 
partially ordered set with n keys, and then we use Stirling’s formula again. ■ 
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In the next Chapter we proceed to the analysis of partial orders where the 
information-theory lower bound is not uj(n). Central to the subsequent analy¬ 
sis are random graphs. 



Chapter 6 


Linear extensions of random graph 
orders 


Recall that in the previous Chapter, we examined the number of linear exten¬ 
sions of various partial orders, where the information-theoretic lower bound 
dominated the linear term. In this Chapter, we examine the case where both 
terms are asymptotically equivalent. We obtain bounds of the expected height 
of a random graph and we derive a new bound on the number of linear exten¬ 
sions of a random graph order. 


6.1 Random graph orders 

In this section, we consider random graphs. A definition follows [37]: 

Definition 6.1.1. The Erdos-Renyi random graph G(n,p ) has labelled vertex 
set { 1 , 2 ,..., n} and for each pair of vertices, the probability of an edge arising 
between them is p, independently of all the other pairs of vertices. 

By the definition, here and throughout this thesis, a random graph G(n,p ) 
denotes a simple graph, without loops or multiple edges. An independent set 
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of a graph is a subset of the vertex set, such that there is no edge connecting 
any two vertices. On the other hand, a clique of a random graph G(n,p ) is 
a subset of its vertex set, with the property that an edge is arising between 
every two vertices. Note that p may very well depend on n. We will usually be 
interested in the behaviour as n —> oo. 

Definition 6.1.2. The random graph order P(n,p), with partial order relation 
-<, is a partially ordered set with underlying set the vertices of G(n,p) and we 
initially say that i -< j if and only if i < j and the edge i ~ j is present in the 
random graph G(n,p). We then take the transitive closure of this relation to get 
a partial order. 

In the same manner, we write P(Z,p) for the infinite partially ordered set ob¬ 
tained by taking vertex set Z, the set of integer numbers, saying initially i -< j 
if and only if i < j in the usual total order on Z and the edge i ~ j is present 
(which it is with probability p independent of all other edges), and then taking 
the transitive closure. 

The point is that for a partial order we of course require transitivity by the 
axioms for a partially ordered set. But of course this is not guaranteed in a 
random graph. We could, for example, have the edges 1 ~ 2 and 2 ~ 3 in 
the random graph, but no edge between 1 and 3. Then of course we would 
have put in the edge 1 ~ 3 as since 1 -< 2 and 2 -< 3 we must have 1 -< 3 by 
transitivity. Note that there are efficient algorithms for finding the transitive 
closure of a relation. What we aim to do next, following the analysis in the 
last section, is to consider how many comparisons will be needed to finalise 
the order of a set of keys when a random graph partial order on the set is given 
already. For our needs, we present the definition of graph entropy, introduced 
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by Korner [47]. This definition is from Simonyi’s survey on graph entropy 
[67]. 

Definition 6.1.3. Let G = G(n,p ) be a random graph. Let X be a random 
variable taking its values on the vertices of G and Y taking its values on the 
stable (independent) sets of G. Suppose further that their joint distribution is 
such that X e Y with probability 1. Also, the marginal distribution of X on 
V(G ) is identical to the given distribution P. Then, the graph entropy H(G, P) 
of the random graph G is 


H(G,P) = min/(X AF), 
where I(X A Y) is as in Definition 3.1.7. 

As in the last Chapter, we need to know about the number of linear extensions. 
The following Theorem from [4] will give us what we need. 

Theorem 6.1.4 (Alon et al. [4]). Let 0 < p < 1 be fixed and consider e(P), 
where the partial order is from P(n,p ). Then we have that there are p(p) > 0 
and <j 2 {p) > 0 such that 

i°g .(.(q)-W-») A 

a{p) ■ y/n 

Corollary 6.1.5. There is a constant c(p) such that whp. we have 

log 2 (e(P)) = c(p)n + 0{n}^ 2 u(n)) j: 

where u(n) is any function tending to infinity with n (we usually think of it as 
doing so very slowly). 
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Proof. For random variable N(0,1), the probability that it is between —u(n) 
and u(n) is 1 — o(l) as n —» oo. Thus for large enough n 


log e (e(P)) -fj.jp) -n 
a(p) ■ yfn 


G (— uj(n),cu(n)). 


Changing the base of the logarithm, we get: 


log 2 (e(P)) 

log 2 (e) 


- M ■ n 


G (— u(n),co(n)) 


a(p) ■ ^fn 

log 2 (e(P)) e (log 2 (e) • ((i(p)n - w(n)a(p)y/n) ,log 2 (e) • (fi(p)n + u(n)a(p)y/n) 


and this gives the claim, with c{p) = n(p) ■ log 2 (e) > 0 . 

Corollary 6.1.6. For a random partial order P(n,p ) with 0 < p < 1 constant, 


log 2 (e(P)) = log 2 (e)//(p)n(l + o(l)). 


Proof. The proof follows directly from the previous Corollary ■ 

In other words, the logarithm of the number of linear extensions is linear in 
n. Thus when we use Theorem 5.2.8, we see that both terms in it log 2 (e(P)) 
and 2 n are linear. So we do not get the exact asymptotics as in the previous 
Chapter. Recall that in that Chapter, we always had log 2 (e(P)) = u(n) so out¬ 
weighed the linear term. (In many cases, log 2 (e(P)) was of order of magnitude 
n log e (n) so won comfortably). However here it is not clear what multiple of n 
will be the time complexity of finding the total order by pairwise comparisons. 
Note that it will be at most some constant multiple of n, so this will be quicker 
than just using about 2 n log e (n) comparisons in Quicksort - in other words, the 
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partial order here does substantially speed up the process of finding the true 
order. 

To deal with this question in more detail, we need to know about the structure 
of a random graph order. We concentrate to begin with on the case where p is 
a constant. In this case, we shall see that basically the partial order consists of 
a linear sum of smaller partial orders. 


6.2 Additive parameters and decomposing posets 

Definition 6.2.1. A vertex v in any partial order is said to be a post if and only 
if every other vertex of the partial order is comparable with it. 

If v is a post in a partial order (P, -<), then we can write the partial order as a 
linear sum of two subposets (Pi, -<) and (P 2 , -<), namely 

Pi = {x E P : x ^ v} P 2 = {x E P : x >- v}. 

Clearly two elements in the same Pi which were comparable before still are, 
and any element x in Pi is smaller than any element y in P 2 since x -< v -< y 
by assumption. 

Definition 6.2.2. A parameter f of partial orders is said to be additive if and 
only if whenever P is the linear sum of two subposets Pi and P 2 we have /(P) = 

/(Pi) + /(P 2 ). 

An example of an additive parameter is the height, as if we have a longest 
chain in P x and a longest chain in P 2 , we can concatenate them to form a 
chain in the linear sum, so /(P) > /(Pi) + /(P 2 ): and in the other direction, 
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given a longest chain in P, we restrict to the subsets and get chains in Pi and 
P 2 , so f(P) < /(Pi) + /(P 2 ), and so they are equal. Another example is the 
number of further comparisons needed to sort a partially ordered set which is 
a linear sum of two subposets. An important Lemma follows 

Lemma 6.2.3. The number of further comparisons c(P) of Quicksort which need 
to be applied to aposet P to obtain the true ordering, with knowledge of where the 
posts are and where elements are relative to the posts, is an additive parameter. 

Proof. Write P = Pi © P 2 . If we sort the whole linear sum of P 1 and P 2 with 
k comparisons, we have sorted Pi and P 2 as well. Thus, taking k = c(P), we 
see that c(P) > c(Pi) + c(P 2 ). (We have sorted Pi and P 2 using k comparisons: 
we might have done it with fewer). Conversely, if we have sorted P x and P 2 
with a total of k comparisons, then we have sorted the whole of P because, by 
definition, everything in Pi is above everything in P 2 , thus c(P) < c(Pi)+c(P 2 ). 
The result follows. ■ 

Now here is a key result from Brightwell [10], 

Theorem 6.2.4 (Brightwell [10]). With probability 1, the set of posts in P(Z,p) 
for 0 < p < 1 constant, is infinite. 

Basically, each bit between posts will be small. Indeed Brightwell’s survey [10] 
also shows that for sufficiently large k, given our p (which remember is con¬ 
stant) there is a constant dependent of p, c(p) > 1 such that the probability 
that none of {2k, 4 k,..., 2k 2 } are posts is less than or equal to c~ k . We can 
now start showing how to use this idea to break down various invariants of 
P(n,p) into small units. We need some notation about posts. Let their posi- 
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tions be 


• • • U—\, Uq, U \... 


where U 0 is the first post at or to the right of 0. Then we say that P, is the 
poset induced on the interval ( Uj , U ]+ \\. These posets are called the factors of 
the partial order. The next Theorem from [10] presents an important result, 
regarding convergence in distribution of additive parameters. 

Theorem 6.2.5 (Brightwell [10]). Let pbe a constant with 0 < p < 1. Let f be 
an additive parameter of partial orders which is not proportional to \P\. Let Y 
and Z be the random variables /(P 0 ) and f(P-i) respectively. Further, suppose 
that the moments E(F r ) and E (Z r ) are finite for all r e N. Then there exist 
constants p = pip) > 0 and a = o(p) > 0, such that E(/(P(n,p)))/n —>• p and 
Var(/(P(n,p)))/n —> a 2 . Furthermore 


f(P(n,p)) - Pip) -n 

a(p) ■ y/n 



N(0,1), 


with convergence of all moments. 

The following Corollary is a consequence of Theorem 6.2.5. 

Corollary 6.2.6. Given 0 < p < 1 constant, the height of P(n,p), which is an 
additive parameter is whp. equal to p(p) ■ n ■ (l + o(l)). 


Proof. This follows easily from the previous Theorem. 
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6.3 Average height of random graph orders 


What we really need to do now is to obtain bounds for the average number 
of linear extensions. Albert and Frieze [3], derived estimates for the aver¬ 
age height of a random graph order, which is an additive parameter as we 
previously saw. The idea is the following. A random graph is sequentially con¬ 
structed; at each step a new vertex j is added and the probability of an edge 
from it to any previously existing vertex i with i < j is 1/2. They consider 
both an underestimate of the height and an overestimate. We present a gener¬ 
alisation of this construction with a constant probability p e (0,1) of an edge 
arising. 

Theorem 6.3.1. The underestimate f(p ) and overestimate h(p ) increments of the 
average height are given by: 


f(P ) 

h(p) 



Proof. Let l k be the length of the longest chain in a random graph order of 
size k e N and d k be the number of top endpoints of longest chains. Consider 
the addition of vertex (k + 1). The event that the new vertex (k + 1) is the 
new, unique, endpoint of a longest chain is the event that one or more of the 
edges from (k + 1) to the d k endpoints actually arises. In this event, what will 
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happen is that the length of the longest chain will increase by 1 and the number 
of endpoints of longest chains will drop to 1, with probability 1 — (1 — p) dk . The 
complementary event is that the number of endpoints will be increased by one. 

It is at the next step that we make a pessimistic assumption. The pessimistic 
assumption is that, in these cases where (k + 1) does not become the unique 
endpoint of a longest chain, the number of endpoints increases by 1 with prob¬ 
ability only p. In fact, though there will certainly be at least one vertex one 
level below all the d k upper endpoints, in most cases there will be more than 
that - say r of them - so if (k+1 ) is joined to any of them the number of longest 
chains will increase, and the probability of this happening will be 1 — (1 — p) r . 
In this event, the number of endpoints will increase by at least 1, but the length 
of the longest chain will remain unchanged. 

Let the random variables rj and 0 denote the underestimates of the length of 
the longest chain and the number of endpoints of the longest chain(s) respec¬ 
tively. These variables obey the following recurrence: 


( 


iyik+i) Ok+i) \ 


(Vk + 1 , 1 ) 
(Vky @k + 1 ) 

(Vk, o k ) 


with probability 1 — (1 — p)° k 
with probability p(l — p)° k 
with probability (1 — p)(l — p) 6k 


Also, let p and (j> be the overestimates of the length of the longest chain and 
the number of endpoints of the longest chain(s) respectively. In this case, the 
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recurrence relation is: 

{ (/j fe + 1,1) with probability 1 — (1 — p)^ k 
4>k + 1) with probability (1 — p)^ k 

We consider only the second component. This is clearly a positive recurrent, 
irreducible, aperiodic, Markov process with state space the positive integers. 
Thus a stationary distribution does exist with limiting probabilities 

Pj = lim P (9 n = j ) 

n—> oo 

and when 9 k+i = j + 1, then this could come about from 9 k being j (with 
probability p(l — p) j ) or from 9 k being j + 1 (with probability (1 — p)b‘+ 1 )+ 1 = 
(1 — p) j+2 )- The solution giving the stationary distribution pj = lim ¥(9 n = j) 

n—>• oo 

for the first case is 


Pj+i = P{ 1 - p) J Pj + (1 - p)(i - p) 3+1 Pj +i 


Pj +1 


P(1 ~ P) j _ 
1 — (1 — p)l +2 


/b • i= n 


Z=1 


p(l — p)* 

1 — (1 — p) i+2 


■Pi- 


The value of pi can be found by the following equation 


oo /j —1 

n 

3= 1 \i=l 


P(l-P)' 

1 — (1 — p) i+2 J 


1 . 
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The expected height increment generally for p e (0,1) is 


OO / ,j —1 


p( i — p ) 1 






Substituting p l3 which is a function of p, the average height increment is 


OO / ,j — 1 


/(p) = 


OO /j—1 

e n 


p( i — py 




i+2 






i+2 


which is equal to 


1 



p(i-pY \ 
i - (i -pY +2 J 


(i -pY 


OO / j—1 

an 

3= 1 \i=l 


i - (i -pY+ 2 J 


For the second (overestimate) case, the stationary distribution obeys the fol¬ 
lowing recursive relation 


Pi+i = (! “ PYPj 

7(7 I 1) 

Pi+i = (! - pY ■ (! - pY _1 ••■■■(!- P)Pi = (! - P)^~Pi- 


Thus, the value of p\ can be retrieved by the following equation 


OO 


Pi = Pi 

3 =1 


~P) 

3 =1 


JU- 1) 
2 


1 
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and the average height increment is 


OO / 

Pi^U 1 - i 1 - (! ~P) J ) 

3 =1 


which further simplified yields the simple expression 


E /,, - 3 ( 3 - 1 ) j(j+i) \ 

(Jl -p) 2 - (1 -P) 2 J — Pi 
3 =1 


Therefore, the average height increment in the overestimate case is 


1 

h (p) = Si -• 

3= 1 


The argument is complete. ■ 


We can rewrite the inverse of the overestimate as 


^(l-p) jU 2 > = (1 - p)° + (1 - p) 1 + (1 - p) 1+ ' 2 + (1 - p) 1+2+3 + . . . 

3 =1 

= 1 + (1 —p) + (1 ~p)( 1 -p) 2 + (1 -p)(l -P) 2 (l -P) 3 + • • • 

00 (j-1 \ 

= e np-p)' ■ 

i=i \i=i / 


The average height increment can be computed in terms of 6 elliptic func¬ 


tions. This class of functions arises in many different areas of Mathematics 


and a comprehensive account of their analyses can be found at [1], [74]. The 


overestimate height increment assumes the simple form of 
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where the 0- 2 function is defined by, (see section 16.27 in Abramowitz and 
Stegun [1]), 

OO 

9 2 {z, q) = 2 q x ' A ( l n{n+1) cos((2 n + 1 )z ), 

71=0 

with q = 1 — p. 

The previous equation can be expressed as a product, instead of summation. 
It holds [74], 

OO 

6 2 {z, q ) = 2c? 1/4 G cos(^) n (! + 2 q 2n cos(2 z) + q An ) , 

n=l 

oo 

where G = lb 1 —q 2n ). Therefore, the expected overestimate height increment 

n =1 

becomes 


1 


n<i - (i - ?)”) n (i+2<i - pt +(i - p) 2 ”) 

n =1 n=l 

l 

n((i-(i- P r)(i + (i- P )”) 2 ) 

n =i ^ 2 

The following Maple graph plots the underestimate and the overestimate 
functions. Note that for p > 0.7, both f(p) and the overestimate h{p) are 
in fact very close to p. 
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p 

Figure 6.1: Plot of underestimate (blue) and overestimate (red) functions of 
the expected height increment of a random graph order. 


In the next section, we give sharp bounds for the average height. 


6.4 Bounds on the height of random graph or¬ 
ders 

We have managed to obtain functions f(p), which is the underestimate incre¬ 
ment of the expected height and h(p), which is the overestimate of the ex¬ 
pected height of a random graph partial order. We first give a lower bound on 
the underestimate function, which is easier to work with (more tractable). 

Lemma 6.4.1. For all p e (0,1), f(p) > p. 
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Proof. The probability of an edge arising between a newly added vertex [k+ 1) 
and any of the k existing vertices in a random graph is p, independently of 
the other edges. Thus, with probability p the length of a greedy chain will 
increase by one and with probability (1 — p), the length will remain unchanged. 
Thus, the expected increase in the height of a chain is p and the claim follows 
immediately. ■ 

Here is an upper bound on the overestimate function. 

Lemma 6.4.2. It holds that 


h(p) <p + 


(P- l ) 2 

2 — p 


Proof. To obtain an upper bound on 


hip) 


1 

1 + (1 -p) + (1 - p) 3 + ... 


it is enough to bound below the denominator 1 + (1 — p) + (1 — p) 3 + ... . 
A crude lower bound is 2 — p (as all other terms are positive). Thus 


1 

2 — p' 


h(p) < 
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This in turn is 


1 

2 — p 


= 1 + (p - 1) + (p - l) 2 + ... 


l-(p-l) 

p + (p - l) 2 + (p - l) 3 + ... 
(P- l ) 2 


p + 
p + 


1 - (P - 1) 
(P ~ l ) 2 

2 — p 


Note that as p tends from the left to 1 and this fact explains the convergence of 
the underestimate and overestimate, seen previously in Figure 6.1. Collecting 
more terms, we will be able to obtain a sharper bound in the following manner. 
It holds that 


1 + (1 - p) + (1 - p) + ... > 2 - p + (1 - py + ... + (1 - p) , 


jYl( Tfl — 1 ) 

for a finite number K = -, where m e N, so the bound has the 


form 


1 

m 

3 =3 

Its sharper than the previous one but again infinitely many terms are dis¬ 
carded. In order to obtain a bound using all terms, note that j(j — l)/2 < j 2 
and 


1 + (1 - p) + (1 - pf + ... > (1 - p) + (1 - p) A + (1 - pf + ... . 
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The overestimate can be bounded above by another theta function. Therefore, 
a sharper bound is 


,I(P) * «,( 0 , 1 -„) + !• 

The 0 3 function is defined [1], 

OO 

6*3 (z, q) = 1 + 2 ^ g” 2 cos( 2 nz). 

n= 1 

Thus, we obtained suitable bounds, such that 


p<f(p)< h(p ) < 


e 3 (o,i- P ) + i 


In the following section, we present bounds on the number of linear exten¬ 
sions. 


6.5 Expected number of linear extensions 

Here we obtain bounds on the average number of linear extensions of a ran¬ 
dom partial order. We start by quoting a useful Theorem from [4], 
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Theorem 6.5.1 (Alon et al. [4]). Let a random variable Y be geometrically 
distributed, i.e. P(Y = k) = pq k ~ l , for k = 1,2, .... Then, 


E(io ge (y» 


1 •' /'O') 

1/p — k(p) 


k= 1 


< (! - Hp)) lo S e 


1 — k(p) 


1 


< log e - = log e (E(Y)), 


where k(p) is defined as, 

OO 

kp)= 

k =1 

with q—l—p. 

The following Theorem regarding the expected number of linear extensions is 
the main contribution in this Chapter. 


Theorem 6.5.2. Let p e (0,1). The average increment of the natural logarithm 

E(log (e(P))) 

of the number of linear extensions pip) =--— - - is whp. bounded below 


n 


by: 


Pip) > 


log 2 h(p) log e (2) 


Proof. We consider the longest chain C. We know that whp. E(|C|) < h(p)n 
by Theorem 6.3.1. Now, by Lemma 5 of Cardinal et al. [13], which states that 
\C\ > 2 ~ H ^n, we deduce that 


hip ) > E(2-^ (p) ) t 
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where H(P) denotes the entropy of the incomparability graph. By Lemma 4 of 
that paper, we also have that 

nH(P) < 2 log 2 (e(P)). 


Therefore 


h(p) >E (e(P)~ 2/n ). 

Taking on both sides logarithms and applying Jensen’s inequality [38], we get 
that whp. 


log, hip) > log, E (e(P) 2/ ”) 

> E (log 2 (e(J>)- 2/ ")) 
= -|E(log 2 (e(P)) 

= -^bP p) ' 


Thus, n(p) is whp. bounded below by 

^ > log 2 %) log e (2) 

and the proof of the Theorem is complete. ■ 

As we can see from the following graph, the bound isn’t tight compared with 
the one of Alon et al. [4]. However, its derivation provides an insight into 
the relation of the height of a random graph order with the number of linear 


extensions. 
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Figure 6.2: Plot of Alon et al. [4] bound (red) and of bound of Theorem 6.5.2 
(green) on fi(jp). 




Chapter 7 


Conclusions and future research di¬ 
rections 


At the last Chapter of this thesis, the conclusions of the research and possi¬ 
ble future directions are discussed. In the first section, we consider the sort¬ 
ing of partially ordered sets and in the second section, the fast merging of 
chains. 


7.1 Sorting partially ordered arrays 

Central to the analysis of the time complexity of sorting partially ordered sets, 
was the number of linear extensions, as a measure of the ‘presortedness’ of 
the array. Recall that the quantity log 2 (n!) is the lower bound of compar¬ 
isons needed to sort an array of n keys, with no prior information. In all 
cases of partially ordered sets considered in Chapter 5, the constant 1/2 log 2 (e) 
appeared to the asymptotic number of comparisons. A future direction to re¬ 
search might be the sharpening of these results. For example, one might ask, 
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what is the average number of key exchanges or the computation of exact 
expected costs. 

Generalising Albert-Frieze argument [3] and using entropy arguments, a new 
result was the lower bound on the number of linear extensions of a random 
graph order. However, we have seen that it does not directly compete the 
bounds of Alon et al. [4], thus there is space for further improvement of this 
bound or to derivation of new sharper ones and this might be a suitable topic 
for further research. 

A different bound on the number of linear extensions of a random interval 
order can be derived as follows. By Theorem 5.6.5, one can deduce that whp., 
the size of the longest chain C, is 


|C1 


2 y/n 




Using the result \C\ > 2 H ^n from Cardinal et al, as in Theorem 6.5.2, we 
have 


2 > 2 _ h{P) 
\fnn 


Taking logarithms, a lower bound for the entropy of the incomparability ran¬ 
dom interval graph is, 
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By the following inequality [13] 

nH(P) < 21og 2 (e(P)) , 
log 2 (e(P)) is whp. bounded below by 

77 / 

log 2 (e(P)) > — log 2 (n) + 0(n). 

Thus, whp. 

77 / 

l°g 2 (e(P)) > -log 2 (n)(l + o(l)) 

= 4l^(2) l0g » (n) ( 1+0(1) ) 

« 0.36n log e (n) (l + o(l)). 

Remark 7.1.1. The speed up factor of the derived bound of the number of linear 
extensions of a random internal order is l/41og e (2) « 0.36 - exactly half of the 
constant stated in Corollaiy 5.6.6. This fact shows that the bound in this Chapter 
performs rather poorly, comparing with the results in section 5.6. Note that, with 
further information about the entropy, improvement might be possible. 

7.2 Merging chains using Shellsort 


In the last section, we consider the merging of chains. This problem is anal¬ 
ysed in the paper of Hwang and Lin [32], where an efficient algorithm is 
presented. 
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Here, we discuss some preliminary ideas, that might be worthwhile for further 
study. Specifically, we propose the application of Shellsort for the merging of 
linearly ordered sets. Shellsort was invented by Donald Shell [66] in 1959 
and is based on insertion sort. The algorithm runs from left to right, by com¬ 
paring elements at a given gap or increment d e N and exchanging them, 
if they are in reverse order, so in the array {cti,a 2 ,.,. , a n } the d subarrays 
{aj, dj + d, aj + 2 di ■ • • }, for j = 1,2,... ,d are separately sorted. At the second 
pass, Shellsort runs on smaller increment, until after a number of passes, the 
increment becomes d — 1. This final insertion sort completes the sorting of the 
array. 

The sequence of the increments is crucial for the running time of the algorithm, 
as the pivot selection is important to Quicksort. Shell [66] proposed the se¬ 
quence Li J , Lf J > • • ■ ’ w hich leads to quadratic time. Pratt [57] suggested a 
sequence of the form 2 a 3 6 < n, where a, b e N, which yields Q (n log^n)) time. 
Incerpi and Sedgewick [36] have shown that there do exist log 0 (n) increments, 
for which the running time of the algorithm is 0{n V 1 °s 2 m with a = 2 e ' 8 
and e > 0. Despite the extensive analysis of Shellsort, there are many open 
problems, as whether the algorithm can achieve on the average O (n log 2 (n)) 
run-time. 

In our problem, Shellsort can be fruitfully applied to merging chains, which 
can be done quite fast, using the knowledge of the partial order. Consider two 
chains C\ and C' 2 , 


Ci = {ai < a 2 < ■ ■ ■ < a n j 
C'2 = {b\ <b 2 < ... < b m }. 
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Starting from a\, with initial increment d = rna x{n,m}, the algorithm sepa¬ 
rately sorts d subarrays. At the second pass, the algorithm iterates from a 2 
with increment d — 1. The final comparison of the element occupying the lo¬ 
cation d with its adjacent key in the position d + 1, terminates the merging 
process. 

We illustrate this informal idea, with a simple example. Suppose that we want 
to merge the chains C\ = {5, 7,9,11,12} and C 2 = {4,6,10}. We start with the 
array {5,7,9,11,12,4,6,10} and initial increment max{5,3} = 5. Then, the 
following subarrays are independently sorted: {5,4}, {7,6}, {9,10}, so at the 
end of the first iteration, the array becomes {4,6,9,11,12,5, 7,10}. Starting 
from 6, with increment equal to 4, the algorithm proceeds to the subarrays 
{6,5}, {9, 7}, {11,10}, so we obtain {4,5, 7,10,12,6,9,11}. In the same man¬ 
ner, starting from 7, with gap equal to 3, the subarrays {7,6}, {10,9}, {12,11} 
are sorted, giving {4,5,6,9,11,7,10,12}. Then, the subarrays {9,7,12} and 
{11,10} are sorted, yielding {4,5,6,7,10,9,11,12}. The final comparison to 
{10,9} returns the merged chain. This algorithm took 13 comparisons for the 
merging of 8 keys. A different increment sequence might speed up the pro¬ 
cess, noting that there are some redundant comparisons, e.g. the comparison 
of {9, 7} in the second pass. Its appealing feature is that it merges ‘in-place’, 
without the need of auxiliary memory. 

Central to the argument is the number of inversions of a permutation of n 
distinct keys {ai, a 2 ,..., a n }. An inversion is a pair of elements, such that 
for i < j, ^ > aj. Obviously, an upper bound for the number of inversions 
is 1 + 2 + ... + (n — 1) = . On the other hand, a sorted array has 0 

inversions. In other words, the number of inversions determine the ‘amount’ 
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of work needed for the complete sorting. Generally, when one has to merge k 
chains, with cardinalities ..., m k and rrij = n, an upper bound to 

the number of inversions is 



Note that this bound corresponds to the case, where the chains are presented 
in completely wrong order, e.g. the elements of a chain Cj are greater from the 
elements of chains, which lie to its right. In practice, this case occurs rarely, so 
the bound can be greatly improved. 

The application of Mergesort, as proposed by Cardinal et al. [13] for the merg¬ 
ing of chains completes the sorting in (1 + e) log 2 (e(P)) +0(n) time - see their 
Theorem 3. Cardinal et al. remark in their paper that their Mergesort algo¬ 
rithm is better than an earlier one of Kahn and Kim [40], provided log 2 (e(P)) 
is super-linear. As we have seen in this Chapter, log 2 (e(P)) is linear, thus the 
application of Shellsort might constitute an alternative choice for the merging 
of chains. 



Appendix A 


Maple Computations 


Here is the Maple worksheet for the computation of the variance of the num¬ 
ber of key comparisons of dual pivot Quicksort. 

> restart; 

> f_{n}(z):= 1/binomial(n, 2)sum(sum(z''2n-i-2f_{i-l}(z)f_{j-i-l}(z) 
f_{n-j}(z), j=i+l..n), i=l..n-l); 

> diff(f_{n}(z), z$2) ; 

> subs(z = 1, diff(f_{n}(z), z$2)); 


n— 1 n 


n—1 n 


/ " (1) = ( E E < 2 " - > - 2) 2 -EE ( 2 « - i - 2 ) 

' ' ' i =1 j=i +1 i= l j=i+l 

n—1 n n—1 n 

+ 2 EE (2n - i - 2)E(C i _i, 2 ) + 2 EE (2 n - i - 2)E(Cj-_i_i j2 ) 

2 = 1 j=i+1 i= 1 j=i+1 

n—1 n n—1 n 

+ 2 EE (2n - i - 2)E(C n _ J - 2 ) + 2 EE E(C i _ li2 )E(C , J -_ i _ li2 ) 

2 = 1 j=i -\-1 z=l j=i-\-1 

n—1 n n—1 n 

+ 2 EE E(C , i_i,2)E(C' n _j i 2) + 2 EE E(C J -_ i _ li2 )E(C n _ i)2 ) 

2=1 j=i-\- 1 2=1 J=2+l 

n—1 n n—1 n n—1 n \ 

+ fi-ii 1 ) + E E fj-t-ii 1 ) + fn-A 1 )J • (A- 1 ) 

2=1 J=2+l 2=1 J=2+l 2=1 j =2-f-l ' 


175 



176 


The following Maple commands compute the sums of Eq. (A.l) 

> sum(sum((2n-i-2)~2, j=i+l..n), i=l..n-l); 

> sum(sum((2n-i-2), j=i+l..n), i=l..n-l); 

> 2(sum(sum((2n-i-2)(2iharmonic(i-l)-4(i-l)), j = i+l..n), i = l..n-l)); 

> simplify(°/„) ; 

> 2(sum(sum((2n-i-2)((2(j-i))harmonic(j-i-l)-4(j-i-l)), 
j=i+l..n), i=l..n-1)); 

> sum(4 n harmonic(k)(k+i+1)-4nharmonic(k)i-2iharmonic(k)(k+i+1) 
+2harmonic(k) i~2-4harmonic(k)(k+i+1)+4iharmonic(k), k=0..n-i-1); 

> 2(sum(3n-3i+2nharmonic(n-i)(n-i)~2 -iharmonic(n-i)(n-i) 
+2nharmonic(n-i)(n-i)-(3(n-i))n+(3/2(n-i))i-iharmonic(n-i)(n-i)~2 
-2harmonic(n-i) (n-i)~2 

+(1/2)i(n-i)~2-2harmonic(n-i)(n-i)-n (n-i)~2+(n-i)~2 , i=l..n-1)); 

> 2(sum(4nharmonic(j)(n-j)~2-5n~2harmonic(j)(n-j)-2nharmonic(j) 

+(n-j)harmonic(j)n+(2(n-j))harmonic(j)-(n-j)~2 

harmonic(j)+2n~3 harmonic(j)-(n-j)~3 harmonic(j), j=l..n-l)); 

> sum((2(k+l))harmonic(k)-4k, k = 0..n-i-1); 

> 2(sum((2iharmonic(i-l)-4(i-l))(2binomial(n-i+1, 2) 
harmonic(n-i)+(n-i-5(n-i)~2)(1/2)), i=l..n-l)); 

> sum((8(n-k))binomial(k+1, 2)harmonic(k), k=l..n-l); 

> sum(8binomial(k+1, 2)harmonic(k), k=l..n-l); 
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Using these results, the second-order derivative evaluated at 1 is given recur¬ 
sively by: 


a 1) = 2 (n + l)(n + ) - H n (h„ 2 + fn + e) + 

n—1 


209 

~36” 


-n 


731 13 6 

~36 U + If + n(n - 1) ^ 

v 7 1=1 


(A.2) 


The Maple input for the solution of the recurrence is as follows 

> f’’_n(l)=2 (harmonic(n))~2 n~2-17/3 harmonic(n) n~2+209/36 n~2 

-2 harmonic(n,2) n~2+731/36 n-6 harmonic(n,2) n-47/3 harmonic(n) n 
+6 (harmonic(n))~2 n+4 (harmonic(n))~2-4 harmonic(n,2)-6 harmonic(n) 
+13/6+6/(n(n-l))sum((n-i)f’’_i-l(1), i=l..n-1); 

> binomial(n, 2)f’’_{n}(l); 

> binomial(n+1, 2)f’’_{n+l}(l)- binomial(n, 2)f’’_{n}(l); 

> binomial(n+2, 2)f”_{n+2}(l)-2binomial(n+1, 2)f”_{n+l}(l) 

+ binomial(n, 2)f’’_{n}(l); 

The output of these commands gives the second-order difference 


A' 2 a 1) = 12(n + !)(« + 2)(Hl - Hf>) - H n ( 20n 2 + 32n - 12) 


+ 17n 2 + 37n + 3/"(l) 


(A. 3) 


and after standard manipulations of Eq. (A.3), (recall subsection 4.1.1), the 
solution of the recurrence is 


/"(1) = 4 (n + 1 f{H 2 n+l - H^) - 4 H n+l (n + l)(4n + 3) + 23n 2 + 33 n + 12. 
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